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Abstract—The rate of oxidation of the 110 faces of zinc sulphide has been studied by the use 
VOL. of a thermobalance. Experiments were carried out at temperatures between 680°C and 940°C 
and with cxygen partial pressures between 0-014 and 0-50 atm. 


956/57 Below 830°C the oxidation rate varies as an exponential function of the temperature and is 
largely independent of the thickness of the zinc oxide layer ; it is thus “ chemically controlled.” 
Above 830°C there is evidence of the onset of a diffusion resistance. At temperatures above 
900°C, the thickness of the oxide layer becomes very uneven, due, it is believed, to a phenomenon 
of thermal instability. At still higher temperatures an entirely different mechanism sets in, 
namely the dissociative vaporisation of zinc sulphide, followed by its oxidation in the vapour 


phase. 


Résumé—Les auteurs étudient a l'aide d'une balance thermique, les vitesses d’oxydation des 
faces 110 du sulfure de zinc. Ces essais ont été faits entre 680 et 940°C avec une pression partielle 
d’oxygéne de 0,014 et 0,50 atm. 


Au-dessous de 830°C la vitesse d’oxydation est une fonction exponentielle de la température 
et elle est indépendante de |’épaisseur de la couche d’oxyde de zinc ; elle est alors chimiquement 
contrélée. A des températures supéricures 4 830°C il se produit une résistance due a la diffusion. 
Au-dessus de 900°C la couche d’oxyde de zine devient irréguliére : ceci est dd & un phénoméne 
d'instabilité thermique. Enfin & des températures plus élevées, un mécanisme entiérement 
différent entre en jeu : la vaporisation et la dissociation du sulfure de zinc suivies d’une oxydation 
en phase vapeur. 


1. INTRODUCTION and at present the engineer could not hope to 


Reactions between solids and gases are an 
important class which includes the combustion of 
solid fuels, together with many examples from 
the metallurgical industry, such as the roasting 
and smelting of ores. The study of such processes 
falls within the field of the chemical engineer, 
especially in so far as they involve mass and heat 
transfer taking place simultaneously with the 
chemical reaction. However, there is a lack of 
knowledge concerning the rate limiting factors, 


carry out a successful design for such a process 
in the same a priori manner as he would use for 
a distillation column or a vessel for a homogeneous 
reaction. 

In this paper we shall describe some work on a 
reaction of the above type, the oxidation of zinc 
sulphide, which is of technical importance in zinc 
production. The work is confined to a study of 
the reaction rate, the diffusion of the gases 
through the layer of reaction product, and the 
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conditions of thermal stability. For the purposes 
of a fundamental study, the reaction in question 
has the merit of simplicity, for over a considerable 
range of temperature it proceeds in accordance 
with the equation 


ZnS + 50s — ZnO + SO, 


to the exclusion of other reactions, such as the 
formation of zinc sulphate. The high exother- 
micity of the reaction, on the other hand, gives 
rise to some difficulty in the control of tempera- 
ture. Earlier studies on the same reaction have 
been described by Ocawa [1] and by One, 
WapswortH and Fassexy [2]. 


2. EXPERIMENTAL 


(a) Method. The experimental method consisted in 
following the weight change of a single crystal, a few 
millimetres in size, suspended in a furnace from a sensitive 
balance. The crystal was surrounded by a stream of air, 
or other oxygen-nitrogen mixture, and conditions were 
chosen to avoid mass transfer limitation within the 
gaseous phase. On the other hand, as reaction proceeded, 
the core of zinc sulphide progressively diminished and was 
replaced by an adhering layer of zinc oxide, more porous 
than the sulphide. There was, therefore, the possibility 
of a mass transfer limitation due to the need for oxygen 
to pass inwards, and for sulphur dioxide to pass outwards, 
through this layer. 


(b) Preparation of specimens. It was considered im- 
practical to grow artificial crystals of zinc sulphide of a 
sufficient size and use has therefore been made of the 
best available natural sphalerite, as obtained from the 
Picos de Europa mine at Santander, through the good 
offices of the University of Madrid. This material had a 
yellowish colour, but was transparent in small thick- 
nesses. The main impurities were 0-1-0-3 %, of Fe, 0-1-0-3% 
of Cd and about 1 % of silica.* 

A disadvantage in the use of the naturally occurring 
crystal is that it tends to decrepitate when placed in 
the furnace, probably due to small inclusions of water or 
of free sulphur. To obviate this difficulty, we have used 
for our rate measurements only those crystals which had 
successfully undergone a previous heat treatment. 


* Independent experiments carried out by Mr. J. F. 
PEARSON at these laboratories have shown that a very 
impure black zinc sulphide oxidises at the same rate as a 
fairly pure yellow variety. It appears, therefore, that the 
presence of the commonly occurring impurities in varying 
amounts does not significantly affect the observed reaction 
rate. 


The crystals, a few millimetres in size, as obtained from 
this treatment, usually manifested two or three of the 
natural cleavage planes (110 faces), but were otherwise 
irregular in shape. Before using for the rate measurements, 
these crystals were therefore brought to a regular shape 
by grinding and polishing on a suitable machine having a 
double dividing head with horizontal and vertical rotating 
axes. By using any two of the pre-existing 110 faces on 
the specimen for purposes of orientation, it was thus possi- 
ble to reduce the whole surface to a set of known faces 
(usually chosen as all 110) to within an accuracy of about 
+1°. 

(c) The thermobalance. <A large vertical tube furnace 
was mounted on slides and provided with counterweights 
so that it could be moved up and down to allow loading 
of the crystal without upsetting the balance. The latter 
- an Oertling Model 52M with optical projection scale — was 
supported on a strong shelf, two feet above the furnace, 
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Fic. 1. Schematic layout of thermo-balance. 
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A silica thread hung vertically downwards from one of the 
balance pans, through a hole in the shelf, and supported 
at its lower end, within the furnace, a small platinum tray 
for carrying the crystal, which was supported on a platinum 
wire frame above the tray to allow O, to reach its under- 
side. Details of the arrangement are shown in Figs. 1 
and 2. 

The furnace contained two concentric silica tubes, and 
the current through the primary winding on the outer 
tube was permanently “on.” The secondary winding on 
the upper part of the inner tube, in the vicinity of the 
crystal, was used for control purposes and was activated 
by a thermo-junction within the furnace, together with 
a Cambridge Direct Deflectional Indicating Controller. In 
the reaction zone the temperature remained steady to 
within + 0-5°C and the vertical variation of temperature 
was no more than 5°C per inch; since specimens were 
of the order } in. in size, and thermo-junctions could be 
suspended as close as 4 in. to the specimen, the errors 
due to temperature gradients were probably small. 
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Fic. 2. Detail of reaction zone. 


Radiation shields were provided and thermo-couples 
were periodically recalibrated. 


(d) Control of the gas flow. A prepared nitrogen-oxygen 
mixture, contained in a gas holder, was passed through 
drying tubes, and a flow meter, and thence into the lower 
part of the inner tube of the furnace, where it was heated 
by the ceramic packing shown in Fig. 2. The flow of gas 
(maintained constant at 400 c.c./min) gave rise to a certain 
amount of “ lift’? on the thermo-balance, but this was 
allowed for in its calibration. 


(e) Experimental procedure. The various dimensions of 
the crystal were first determined by use of a travelling 
microscope and its weight was measured on an ordinary 
balance. Meanwhile the furnace had been adjusted to the 
required temperature with an approximate oxygen-nitro- 
gen mixture passing through. The crystal was mounted 
above the thermobalance pan and, at ¢ = 0 for the rate 
measurement, the furnace was lifted on its counterweights 
until it surrounded the crystal. Thereafter the reaction 
rate was followed over a period of hours on the optical 
scale of the thermobalance. 

At the completion of the run, the total weight change, 
as measured on the thermobalance, was checked by trans- 
ferring the crystal, now at least partly zinc oxide, to an 
analytical balance. The crystal was finally cut across in 
order to check the equality of the oxide layer on all 
faces. This was necessary as several specimens were found 
to have contained undetected cracks, which were penetrated 
by the oxygen, thus invalidating the particular experi- 
ment ; it served also to confirm the equal ease of diffusion 
of oxygen up to the various crystal faces. 


(f) Temperature correction. Due to the exothermic 
nature of the reaction, the temperature of the crystal 
was usually somewhat higher than was measured by the 
thermo-junction in its vicinity. The magnitude of this 
effect was investigated by placing a_ thermo-junction 
between a pair of crystals, the leads from the junction 
being wound severa] times round the pair and the whole 
assembly being hung in the furnace in the position normally 
occupied by the platinum tray. It was found that the 
temperature within the pair of crystals exceeded the 
external temperature by as much as 10° at 800°C and 
20° at 900°C, in fair agreement with values calculated 
theoretically (Part II). The temperatures used in con- 
nection with equations (1) and (2) below are the actual 
temperatures within the crystal, as estimated by means 
of this correction. 


8. EXPERIMENTAL RESULTS 


(a) Chemical control. For the rate measure- 
ments of Tables 1 and 2 below, the crystal faces 
were all 110 and thus had known angles between 
them. From a knowledge of the original dimen- 
sions of each individual crystal, it was therefore 
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possible to calculate what would be its weight 
change for a progressively increasing inward 
movement, assumed equal on each face, of the 
ZnO-ZnS interface. This calculated weight 
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Fic. 3. The effect of temperature on the rate of film 


growth in air. 


change could then be plotted against the move- 
ment of the interface, or “* penetration "’ as it will 


be called. From the thermo-balance readings 
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Fic. 4. The effect Uf oxygen concentration on the rate of 
film growth at 780°C. 


each run also gave a plot of the observed weight 
change against time. Thus, from a comparison 
of the two plots, the penetration was obtained as 
a function of the time. 
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. The effect of oxygen concentration on the rate of 
film growth at 850°C. 


These results are shown in Figs. 8, 4 and 5. It 
will be seen that in almost all runs the rate of 
penetration diminished somewhat during the 
early stages of reaction, but subsequently 
approached an almost constant rate, the thickness 
of the oxide layer increasing linearly with the 
time.* The cause of the initial decrease of rate 
is not known, but it appears to be a fairly common 
feature of such processes and a number of possible 
causes have been put forward by GULBRANSEN 
and Anprew [3]. On the other hand, the linearity 
which follows the initial period is consistent with 
the assumption that the observed rate is ‘ chem- 
ically controlled,”’ i.e. the rate of movement of 
the ZnO-ZnS interface depends on the processes 
taking place at this interface and not on the 
diffusion through the overlying layer of zinc 
oxide. 

What will be called the “ steady reaction rate,” 
as measured in centimetres penetration per second, 
has been calculated from the linear portion of 
these curves and the results are shown in Tables 1 


* Onc, Wapswortu and Fasse. [2] report linearity 
over the entire range. However, an inspection of their 
Fig. 3 indicates an initial fall-off in rate in their runs at 
771°C and 799°C. 
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The first of these refers to variable tem- 
perature at constant gas composition (air) and 
the other refers to variable partial 


and 2. 
oxygen 


Table 1. 


Oxidation in air at atmospheric pressure 





7 T 


Calculated 
reaction reaction reaction 
temperature rate rate 
°C | (em/sec) x 108 | (em/sec) x 10% 


Corrected Steady 
Furnace | 
| 


temperature 


0-14 
0-23 
0-41 
0-69 
1-13 1-05 
1-54 1-88 
3-3 3-3 
40 3-4 
3-6 3-7 
42 40 
46 46 
6-0 (9-3) 
8-8 (14-2) 
11-8 (29) 
10-0* (29) 
11-3* (41) 
12-5* (49) 


' 


0-15 
0-24 
0-45 
0-54 





* These results were not obtained by use of the thermo- 
balance, but by direct measurement of the oxide thickness. 


Table 2. Ovidation in oxygen-nitrogen mixtures 


| | 
| Steady 
Oxygen | reaction 
Furnace partial rate 
temp. (ltemperatu re pressure (cm /sec) 


Cc ( atm. 





Calculated 
reaction 
rate 
(cm /sec) 
x 108 
0-41 
0-60 
0-73 


Corrected 
reaction 


0-014 

0-030 
0-040 
0-114 1-25 
0-21 1-88 
0-35 ‘ 2-4 
0-50 , 3-0 





0-014 
0-029 
0-064 


(1-7) 
(2-7) 
(4-2) 
(9-7) 
(12-5) 


0-36 








pressure at two approximately constant tempera- 
ture levels. 

Using the data of Table 1, the logarithm of the 
steady reaction rate is shown in Fig. 6 plotted 
against the reciprocal of the absolute temperature. 
Over the range 680°C to 830°C, the relationship 
is a linear one and the rate can thus be represented 
by the relation : 


E aT (1) 
E = 50 kceal/g mole. However, at tem- 
peratures above 830°C the measured rate is seen 
from Fig. 6 to become less than corresponds to 
this equation and this is believed to be due to 
the onset of an appreciable rate limitation arising 
from diffusion through the oxide layer. 


rate = const X e~ 


with 











reaction rate, crn/sec) 


(R 








+ x 104 (7 = ternp. °K) 


Fic. 6. The effect of temperature on the oxidation rate of 
sphalerite in air, © Thermobalance results ; A Microscope 
measurements. 


Turning now to the effect of the oxygen partial 
pressure, it will be seen that the two parts of 
Table 2 do not refer to temperatures which are 
quite constant. Equation (1) above has been used 
to reduce the rates to temperatures of 780°C and 
850°C and these were then found to be roughly 
proportional to the square root of the oxygen 
partial pressure (Figs. 7 and 8). The combined 
effect of temperature and partial pressure may 
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Fic. 7. The effect of oxygen concentration on the oxidation 
rate of zine sulphide at 780°C. 
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Fic. 8. The effect of oxygen concentration on the oxidation 
rate of zinc sulphide at 850°C. 


thus be represented approximately 
relation* 


by the 


rate = Ae~*/*Tpi 


(2) 
where p is the partial pressure in atmospheres. E 
has the value 50 kcal/g mole; A has the value 
8-0 x 10* when the rate is expressed as centi- 
metres penetration per second, or has the value 


* An alternative expression for the dependence of the 
rate on the oxygen partial pressure will be considered in 
the discussion. 


3-2 x 10° when the rate is expressed as gm. moles 
of ZnS per cm® per second. The rate as calculated 
by the use of equation (2) is shown in the last 
column of Tables 1 and 2; the figures enclosed 
in brackets refer to temperatures above 830°C 
where equation (2) is no longer applicable due to 
the onset of diffusion control. 

The above results apply to zinc sulphide crystals 
having entirely 110 faces. In order to investigate 
the possibility of anisotropy in the roasting rate, 
a series of additional faces were ground on to the 
crystal at 5° intervals between the 110 faces, 
which met at an external angle of 60°. Thus the 
most oblique of the new planes was at 30° to the 
110 faces. After roasting for one hour at 780°C, 
the crystal was cut through a plane normal to the 
previously established faces ; the thickness of the 
oxide layer was found to be greater the larger the 
angle to the 110 face, reaching a maximum of 35% 
greater for the face at 30°. 

Similar results were obtained for more pro- 
longed periods of heating and it appears therefore 
that there is a significant variation jn roasting 
rate from face to face, the cleavage planes (110) 
reacting the most slowly. A 100 plane was cut 
on another crystal and was found to develop an 
oxide layer 10°, thicker than on a 110 plane. 

(b) Onset of diffusion control. At all tempera- 
tures below 830°C, the penetration was very 
uniform over the faces, the ZnO—ZnS interface 
was sharply defined and the corners of the 
sulphide core remained sharp, as shown by the 
tull lines in Fig. 9. On the other hand, at tem- 
peratures above 830°C, the corners began to show 
a rounding as indicated by the dotted line in 
Fig. 9. This is clearly what would be expected if 
the diffusion of oxygen or sulphur dioxide 
through the oxide layer was becoming significant 
as a rate controlling factor. The fact that diffusion 
becomes intrinsically slower than reaction above 





ZnO 


; 2nS 


Fic. 9. Showing the sharp corners resulting from oxidation 
below 830°C, and the rounded corners obtained at higher 
temperatures (dotted line). 
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a certain temperature is due, of coures, to its 
smaller temperature coefficient and this effect has 
previously been demonstrated very clearly in the 
ease of carbon combustion by Tu, Davis and 
Horret [4). 

(c) Uneven penetration. At temperatures over 
900°C the penetration became very uneven over 
different parts of the crystal. By making sectional 
examinations at various stages of roasting, it 
became clear that, whereas the reaction was 
initially uniform, small deviations soon developed 
and magnified themselves during the progress of 
reaction. The characteristic appearance of a 
section, showing also the laminar structure of 
the oxide layer, is shown in Fig. 10. 


Fic. 10. Non-uniform growth pattern of oxide films on 


ZnS above 900°C. 


The uneven penetration was at first thought to 
be due to cracks or other faults in the crystal, but 
this does not explain why the effects become 
marked only above 900°C, It seems more probable 
that it is due to thermal instability : the thickness 
of the oxide layer impedes the escape of the heat of 
reaction which thereby causes a local rise in tem- 
perature and therefore a still further increase in 
reaction rate. This theory will be discussed in 
more detail in Part II. 

The occurrence of uneven penetration made it 
impossible to obtain a satisfactory measure of 
the reaction rate at temperatures much above 
900°C, by use of the thermobalance. Accordingly 
the last three results recorded in Table 1 were 
obtained by direct microscopic measurement of 
the thickness of the oxide layer, over those parts 
of the crystal where penetration was uniform, at 
early stages in roasting. 

(d) Vaporization. When crystals were roasted 
at temperatures above 950°C there was a con- 
siderable deposition of woolly zinc oxide on the 
platinum and silica supports of the roasting tray. 


At these temperatures it is known (5) that zinc 
sulphide has an appreciable dissociation pressure, 
whilst zine oxide itself is still almost involatile. 
It seemed therefore that an evaporation of zinc 
sulphide was taking place, followed by a vapour 
phase oxidation and a deposition of zinc oxide on 
available surfaces. 

The structure of the oxide layer surrounding 
the crystal was also entirely different to that 
which was obtained at lower temperatures. Below 
900°C the oxide layer was laminar, fragile and 
porous (the zinc oxide having a smaller molar 
volume than the zinc sulphide). Above 950°C the 
layer was dense, nonlaminar and quite hard, and 
under electron diffraction was found to have a 
completely random orientation of its microcrystal- 
line structure.* These observations suggested 
that the porous structure which would otherwise 
have been formed was being filled in by the 
deposition of oxide by the vapour phase reaction. 

As noted above, penetration was far from 
uniform at higher temperatures and, over parts 
of the crystal, the thickness of oxide layer was 
sometimes less than at lower temperatures, in the 
same time. This was probably due to its impervious 
structure preventing the inflow of oxygen. Such 
layers commonly showed blisters and eruptions, 
due to the dissociation pressure of the zinc 
sulphide causing them to be lifted off the under- 
lying solid sulphide. 


4. DIFFUSION AND REACTION 
ON A PLANE SURFACE 


We shall put forward an approximate theory of 
simultaneous reaction at the ZnO—ZnS interface 
and of diffusion of oxygen through the overlying 
layer of zinc sulphide. The theory will be based 
on a number of approximations, the main purpose 
being to illustrate the transition from chemical 
control at low temperature to diffusion control 
at higher temperature, with an intervening region 
where both processes are significant. 

For simplicity we shall consider the case of 
plane surfaces and neglect the edge effects. The 
treatment will thus apply approximately to the 

* At lower temperatures the oxide formed on zinc 
sulphide shows preferred orientation [9], [10] and [11]. 
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roasting of a crystal when the penetration is still 
comparatively small. A comparable theory can 
be readily developed for the case of a spherical 
particle, allowing for the progressive diminution 
of the interfacial area. 

Again for simplicity we shall apply Fick’s law 
to the inwards diffusion of oxygen through the 
porous zine oxide, despite the existence of a flow 
of sulphur dioxide in the reverse direction.* The 
concentration gradient will also be assumed to 
be linear ;+ the number oxygen 
diffusing in unit time up to unit area of the 
interface is therefore 


— c)/2, (3) 


of moles of 


D (eg 


where D is the diffusion coefficient and a is the 
oxide layer thickness at time ¢, ¢ and ¢y are the 
oxygen concentrations at the ZnO-ZnS interface 
and at the outer surface of the crystal respectively 
(the latter being equal to the concentration in 
the bulk gas phase, in the absence of any appre- 
ciable mass transfer resistance in this phase). 

The reaction rate at the interface will be taken 
as being approximately proportional to the square 
root of c, as indicated by the experimental 
results. 

Thus the number of moles of oxygen reacting 
in unit time at unit area of the interface is 


ke}, (4) 


where k is a velocity constant. Equating (4) and 
(5) at the steady state, 


ke? — D (eg — ¢)/2 
and hence 
c=¢(1+r 


r4/l+2 r), 
where r is the dimensionless quantity 
r = 2k? /2D%, (6) 


The penetration rate, dx/dt, is obtained by multi- 
plying (4) by V, the molar volume of zinc sulphide 
and by a factor 4, from the stoichiometry of the 
reaction. Thus, from (4) and (5), together with 
these factors, 


* For a more elaborate treatment of such cases, see 
Jost [6). 

+ For a discussion on this point, see Boorn [7] and 
DancKWERrts [8]. 


- = $kVei (lL +r—ry/1+2/r)* (7) 


This expression reduces to simpler forms under 
two limiting conditions : 


(1) If r <1, which occurs at temperatures at 
which & is small, 


dx ; 
7 = tkVed (8) 


This corresponds to “chemical control ” 
since the rate is independent of the thickness, 


or of the value of the diffusion coefficient. 

Ifr 

where k is relatively large, 
DVey 


os . * (9) 


- 1, which occurs at higher temperatures 


which corresponds to “ diffusion control.” 


Integrating (9) 


2* = : D Veet, 


(10) 


and the thickness is thus proportional to the 
square root of the time. 


Below 830°C, the limiting case represented by 
eq. (8) is in satisfactory accord with the present 
results. On the other hand, the “ parabolic "’ law 
(10) has not been encountered, probably because 
a high enough temperature could not be reached 
before the onset of the phenomera of unever 
penetration and vaporisation. For the sam 
reason we have not made any detailed application 
of the more general equation (7), which is applic 
able in the intermediate range of temperature 
where chemical reaction and diffusion are both 
significant as controlling factors. 


5. Discussion 

As mentioned previously, earlier studies on the 
oxidation rate of naturally occurring crystals of 
zine sulphide have been made by Ocawa [1] and 
by Onc, Wapswortn and Fasse.t [2]. At 800°C 
Ocawa obtained a penetration rate in air of 
4-2 x 10° cm" sec as compared to 2-6 x 10° em 
sec’ by interpolation from Fig. 6 of the present 
work, 
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Much more extensive results were obtained by 
Onc, Wapswortn and Fasseit [2]. Under 
similar conditions the rates described by these 
workers are between 25°, and 50%, of those 
obtained in the present work, and thus lie in the 
opposite direction to the figure quoted from 
Ocawa [1]. 

As described above, we have expressed the 
dependence of the “ chemically controlled ”’ rate 
on the oxygen partial pressure in terms of a half 
power. However an inspection of the calculated 


‘ 


results in the upper part of Table 2 shows that 
the representation of the experimental data by 
means of this relationship is by no means good, 
especially at the extremes of partial pressure. The 
results of Onc, Wapswortn and Fasse. [2] - 
published since our own work was completed — 
cover a much wider range of partial pressures. The 
general characteristic of their results at tempera- 
tures below 846°C is that the rate increases more 
and more slowly with increasing partial pressure 
and eventually becomes independent of it. This 
behaviour they interpret as due to the rate 
limiting step being dependent on the concentration 
of adsorbed oxygen molecules on the zine sulphide 
surface, this adsorption being determined by a 
LANGMUIR type isotherm : 


__K,(0,] 
1 + K,(0,) 


(11) 


where @ is the fraction of the surface covered and 
[O,] is the concentration of oxygen in the gas 
phase. This function accounts for their measured 
rates very satisfactorily, and it can also be used 


to express our own results quite as well, if not 
better than the half power relationship.* Never- 
theless the latter is an adequate approximation 
in the region of oxygen partial pressures which are 


* However, at an oxygen concentration at which the 
rates measured by the American workers have become 
independent of the concentration, our own rates are still 
increasing. This indicates « different value of K,. 


of interest in the technical process of zine pro- 
duction and is a simpler relationship for practical 
use. 

As discussed above, it is believed that the 
diffusion of oxygen through the bulk oxide become 
significant as a rate limiting factor at temperatures 
in excess of 830°C, At still higher temperatures, 
however, it seems that an entirely different 
mechanism of reaction sets in, and it may event- 
ually become faster than the mechanism discussed 
above. This is the dissociative evaporation of zine 
sulphide, followed by its oxidation in the vapour 
phase. The rate of this process may be expected 
to depend on temperature according to an 
exponential factor e~**/"", where the activation 
energy E, is at least equal to the latent heat of 
evaporation. According to VesELovski [5] and 
the unpublished work of Dr. Ricuarpst, this 
latent heat is considerably larger* than the value, 
50 kcal/mole, obtained in Section 3 above for the 
activation energy in the low temperature region 
of “ chemical control.’’ This evaporation mech- 
anism may therefore be expected to overtake the 
mechanism discussed above at a sufficiently high 
temperature ; at even higher temperatures this 
evaporative process itself may run into a second 
region of diffusion control, connected with the 
mixing of the reactants in the gaseous phase. It is 
these processes which are probably the most 
significant in the hottest regions of the bed, as 
used in the industrial roasting method. In parti- 
cular, the deposition of zinc oxide from a vapour 
phase reaction may contribute to the mechanism 
of sintering. 

Financial assistance, which enabled this work 
to be done, was given by C.S.LR.O. and the 
Imperial Smelting Corporation, and is gratefully 
acknowledged. 


+ Dr. A. W. Ricnarps of the Imperial Smelting Corpor- 
ation,Avonmouth: 
* VeseLovski's figure is 65 keal per mole. 
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Studies on gas-solid reactions—II 
Causes of thermal instability 


K. J. Cannon* and K. G. Densicut 
The Chemical Engineering Department, Cambridge 


( Received 18 December 1956) 


Abstract—The paper describes two distinctive forms of thermal instability which can make 
their appearance in gas-solid reactions. The first is due to an unstable type of crossing of the 
curves which represent the rate of heat generation by reaction within the solid and the rate 
of heat loss from its surface. In this connection the significance of the concept of ignition tem- 
perature is discussed. The second type of instability is one which affects localised regions of the 
reacting solid ; under conditions where the temperature of the reaction interface becomes higher 
the larger the thickness of the layer of reaction product, any small initial variations in the thickness 
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of this layer will rapidly augment themselves. A theory of this effect is put forward. 


Résumé— Les auteurs décrivent deux formes différentes d'instabilité thermique susceptibles de se 


produire dans des réactions gaz-solides. 


La premiére serait die & un type instable d'intersection des courbes : vitesse de génération 
de chaleur a l’intérieur du solide et vitesse de diperdition de chaleur par sa surface. Les auteurs 
discutent, & ce sujet, du concept de la température d ignition. 


La seconde forme d'instabilité affecterait des régions bien localisées du solide en réaction ; 
si la température de l‘interface de réaction s‘éléve, l'épaisseur de la couche du produit de réaction 
devient plus grande ; de petites différences initiales dans l'épaisseur de cette couche s’amplifient 
rapidement d'elles-mémes. Les auteurs approfondissent cette derniére théorie. 


1. INTRODUCTION 


In Part I [1] we have described experimental 
work on the oxidation rate of zine sulphide and 
during the course of this work it was noticed that 
there are interesting phenomena of thermal 
instability. Since these phenomena may be met 
with throughout the field of gas-solid reactions, 
including combustion, it has seemed worth while 
presenting our conclusions in a separate paper. 

It appears that thermal instability may occur 
in two distinctive forms according to whether it 
is due to (1) an unstable crossing of the curves of 
heat generation and heat loss, (2) the effect of 
solid reaction product in impeding the loss of 
heat. The former affects the solid particle as a 
whole; the latter results in uneven growth of 
the layer of reaction product on different parts 
of the surface. 


2. UnsTaBLe CRrossINnGs 
Phenomena of this kind have been described in 
the literature, notably by van HEERDEN [2] and 
van Loon [3], but are not yet well known. Con- 
sider a reaction of the type 


ZnS + 50s — ZnO + SO, 


As discussed in Section 4 of Part I, there is a 
transition from chemical control to diffusion 
control with rising temperature. Thus at low 
enough temperatures the overall rate is deter- 
mined by the chemical kinetics and increases 
rapidly and exponentially with temperature ; at 
higher temperatures the rate limiting factor 
becomes the diffusion of the gases through the 
solid reaction product and this process has a 
relatively low temperature coefficient. The 


* Present address : Division of Industrial Chemistry, C.S.1.R.O., Fisherman’s Ben, Melbourne. 
+ Present address : The Department of Chemical Technology, King’s Buildings, West Mains Road, Edinburgh. 
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overall rate, as a function of temperature, is 
thus as shown by the sigmoid curve of Fig. 1 
(which has been drawn from equation (7) of 
Part I, using the known temperature dependence 
of k and D). 


approximately independent of temperature, the 


In so far as the heat of reaction is 


rate of heat generation in the reacting solid will 
be of the same form and the sigmoid curve will 
be taken henceforth to refer to the heat generation. 
Consider now the rate of heat loss from the 
solid. If convective transfer to a gas stream is 
the most significant mechanism, the heat loss will 
be approximately proportional to T — T,, where 
T and T, are the temperatures of solid and gas 
respectively. On Fig. 1, the heat loss will thus 
be represented by straight lines* 4B, A'B’, etc., 
Lo, 
ees | 
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3. 1. Rate of heat generation and loss from a reacting 
particle, 


whose gradient is determined by the heat transfer 
coefficient and whose intercepts A, 4’, 
the horizontal axis represent the particular values 
of T, in the system in question. 

The points C and C’ clearly represent stable 


etc. on 


states of low reaction rate in the chemically 
controlled region; they are stable in the sense 
that any momentary rise or fall of temperature 


* If loss of heat by radiation is also a significant factor, 
the straight lines will be replaced by curves which are 
convex towards the temperature axis. 
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Denna 


will cause the heat loss to become greater or less 
respectively than the heat generation and will 
thus cause the system to return to the points C 
or C’. Similarly points FE’, E’’ and E’” represent 
stable states of high reaction rate in the diffusion 
controlled region. 

Point D’, on the other hand, is an unstable 
crossing point; although heat loss is equal to 
heat generation, any small downward variation 
of temperature will cause the system to fall to 
point C’ and any small upward variation will 
cause the system to rise to point EF’, Thus the 
region C’ to E’ cannot be realized under the 
prevailing conditions. 

Consider the introduction of the solid into the 
gas when the latter is at a temperature 4’; if 
the solid is initially cold its temperature will rise 
to the lowest point at which there is a stable 
On the other hand, if the 


solid is initially hot and at a temperature exceeding 


crossing, namely C”’. 


the value corresponding to D’, the temperature 
will rise to EB’, where there is a high reaction rate 
in the diffusion controlled region. Thus D’ may 
(The 


use of this terminology implies that the solid is 


be referred to as an ignition temperature. 
not regarded as being “ ignited” in the region 
of points C and C’ where the reaction rate is 
very low.) 

The introduction of the cold solid into the gas 
when the latter is at a slightly higher temperature 
A’ will result in the solid reaching the high 
temperature corresponding to E’’’. It is evident 
that the line A’ B"’, which is tangential to the 
heat generation curve at D’’, represents a critical 
condition ; D” is the minimum ignition tempera- 
ture of the solid and the gas temperature A” is 
the minimum gas temperature for ignition of the 
cold solid. 

The tangent, (which is not shown), at point F 
also represents a critical condition of a different 
type; the temperature at F is the lowest tem- 
perature, under the prevailing conditions, at 
which there is rapid reaction and is referred to 
by van Loow [3], in his discussion of carbon 
combustion, as the minimum combustion tem- 
perature. 

The actual values of these critical temperatures 
are determined not only by the nature of the 
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particular reactants but also by the partial 
pressure of the reactive gas and by the value of 
the heat transfer coefficient under the prevailing 
conditions, since it is these which determine the 
relative positions and slopes of the curves. The 
notion of ignition temperature is thus by no 
means absolute, but is relative to a particular set 
of conditions. 

It may be noted also that if the heat transfer 
coefficient is large, or if the reaction is such that 
the rate of heat release is always low, the unstable 
crossings may not be obtained. 


8. Tur Errect or THE LAYER 
oF REACTION PrRopuUcT 


As described in Part I, at temperatures above 
900°C small initial variations in the thickness of 
the oxide layer tend to magnify themselves, giving 
rise to a sort of autocatalytic growth of the layer 
on different parts of the surface. This type of 
instability evidently occurs when the local 
reaction rate at the interface is greater the thicker 
is the overlying layer of reaction product. For 
when this condition is satisfied, any regions of the 
particle where the layer is accidentally slightly 
thicker than elsewhere will tend to develop 
relatively still thicker layers. It is therefore of 
interest to examine the conditions under which 
the reaction rate at the interface increases with 
increasing penetration. 

Consider a partially oxidized spherical particle 
of zine sulphide whose radius is ry. Let r; be the 
radius up to the sulphide-oxide interface, i.e. 
ro — 7, is the thickness of the oxide layer. Let 
T, and T, be the temperatures at the outer 
surface and at the interface respectively and let 
Q be the rate of heat generation per unit area of 
the interface due to reaction. 

The total rate of heat generation is 47 r7Q; at 
the stationary state at which the temperatures 
are time independent, this is equal to the rate at 
which heat passes outwards through any spherical 
shell within the zine oxide layer. Thus 


4n r?Q = «= K 4nr” “ ’ 


r 


where T is the temperature at a radius r and K 


is the thermal conductivity of the layer. 
grating between the limits 


we obtain 


(1) 


The spherical particle will be assumed to be 
suspended in a moving gas stream at a tem- 
perature T, and to be “ seeing ” furnace walls at 
a ternperature T,. If the particle may be regarded 
as being “ grey,” the rate of heat loss will be 


4n _¢ {oe (7 — T/) +h (Ty — T,)}; 


where o and « are the Stefan-Boltzmann constant 
and the emissivity respectively and h is the heat 
transfer coefficient. At the stationary state this 
expression must equal 477? Q, the rate of heat 
generation at the interface. Thus 


=oe(Ty — T/) +h(T, — T,) (2) 


If T,> T, —T, and T, = Ty; (as under the 
experimental conditions described in Part I), 
equation (2) may be approximated by 


= (4oe T +h) (T. — Ty) (8) 


Eliminating T, between (1) and (8) we obtain 


. doe T;* + hj 


which gives the value of 7; — T;, the amount by 
which the temperature at the reacting interface 
exceeds the controlled temperature of the environ- 
ment of the particle. 

In this equation the rate of heat generation, Q, 
is itself a function of 7;. It will be assumed for 
simplicity that the chemical kinetics at the inter- 
face remain the controlling factor over the whole 
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temperature range.* (In a more elaborate treat- 
ment it would be preferable to use the equation 
for a spherical particle, analogous to equation (7) 
of Part I, which allows for both reaction and 
diffusion.) With this assumption Q may be 


written 


Q = Be-®/R7,, (5) 


where B is a constant which includes the heat of 


reaction. Substituting (5) in (4) we obtain 


T,—T,= Be~®/*7; "oP -p)+a(l pl, (6) 


J 


where p, the fractional penetration, is defined by 
(7) 
(8) 


—_ r:) To 


Pp = (To 
a = 1/(40eT/?? + h). 


and 


By differentiation of (6), 


Be~®/*7, ("21 —2p) — 2a(1—p) 


ae 
RT? 


(9) 
1 (7; — T;) 

An instability of the type under discussion 
will occur if d7,/dp is positive: for in this case 
the temperature at the interface is higher the 
greater is the penetration. Areas of the solid 
where the oxide layer is already slightly thicker 
than elsewhere will have a higher interface 
temperature and a higher reaction rate, and any 
initial disparities in thickness will thus tend to 
magnify themselves. 

The denominator of (9) is normally positive : 
the condition of instability is therefore 


(10) 


This condition is the more readily satisfied (a) the 
smaller is the degree of penetration, (b) the smaller 


* In the experimental work described in Part I, diffusion 
through the zinc oxide layer began to have an appreciable 
influence at temperatures above 830°C. Nevertheless, the 
chemical reaction remained the main controlling factor to 
above 900°C, It may be noted that under conditions where 
diffusion becomes the controlling factor the instability 
described in this section cannot occur, as may be seen by 
replacing equation (5) above by an equation analogous to 
(9) in Part I. Diffusion control causes the reaction rate to 
diminish with increasing penetration. 


is the thermal conductivity K, (c) the larger is 
the radius r, of the particle, (d) the larger is the 
furnace temperature, the emissivity and the heat 
transfer coefficient, which determine the value of 
the quantity a, defined by equation (8). 

The physical situation may perhaps be under- 
stood clearly by considering the two 
opposing factors : 


most 


(1) with increasing penetration of a spherical 
particle r, diminishes and this tends to cause 
a reduction in the total heat generation 
4nr7Q, within the crystal ; 


with increasing penetration the thickness of 
the oxide layer, r, — 7,, increases and this 
tends to reduce the loss of heat from the 
interface by thermal conduction. 


If the second factor more than offsets the first, 
the interface temperature will increase with 
increasing penetration and thus cause instability. 

The effect of greater heat loss from the surface 
of the particle in favouring instability (item (d) 
above) appears paradoxical, but is due to the fact 
that 7,, the surface temperature, tends to 
diminish with increasing penetration, as may be 
seen from equation (8); however the larger is 
the surface heat loss the smaller is the relative 
magnitude of this decrease and therefore the 
more readily does 7; increase with increasing 
penetration. 

Equation (6) cannot be solved directly for 
T — T,;, as a function of T; and p, since 7’; also 
occurs in the exponential term. However, by 
trial and error solution, a set of values of 7; — T, 
have been obtained which satisfy the equation. 
For this purpose it was necessary to use estimated 
values of the emissivity and thermal conductivity 
of the porous layer of zinc oxide, and also of the 
heat transfer coefficient h.* It was found that a 
reasonable choice of these quantities leads to 
values of T; — T, which are of the same order 
of magnitude as the experimentally measured 
“ self-heating,” as described in Section 2(f) of 
Part I. 


* The values of the parameters E and B, which also 
appear in equation (6), have been taken from the results 
of Part I, together with the known heat of reaction. 
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The family of curves obtained in this way is radius ry = 0-25cm. Due to the uncertainty in 
shown in Fig. 2. which pertain to a particle of the assumed values of the emissivity etc., these 
curves have actually been drawn using a value 
a - for a, as defined by equation (8), which brings the 

calculated self-heating effect into exact agreement 
with the experimentally observed values, i.e. 
about 10° at 800° and upwards of 20° at 900°C. 
It is a useful check on the theory that the same 
value of a, viz. 3-7 x 10°, when applied in the 
expression (10), predicts an onset of instability 
at about 900°, as observed experimentally. This 
point of instability is also made apparent in 
Fig. 2 as the temperature at which there com- 
mences an upward slope of the 7; — T;, curves 
at low degrees of penetration. 

It may be remarked that a self-heating effect 
(or self-cooling, if the process is endothermic) in 

















VOL. gas-solid reactions has been observed previously 
esa res [4], [5], [6], but the theory of instability to which 


it may give rise seems not to have been described 
+0 hitherto. 

p (Fractional penetration) Financial assistance by C.S.1.R.O. and the 

Imperial Smelting Corporation, which enabled 











Fic. 2. Variation of reaction interface temperature with 


penetration. this work to be done, is gratefully acknowledged. 
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Einige neue Verfahrensprinzipien mit Wirbelschichten* 


E. Wicke 
Institut fir Physikalische Chemie der Universitat Hamburg 


( Received 30 October 1956) 


Abstract—-Quantities of heat locally supplied to or produced in fluid beds are very rapidly 
dissipated. Likewise, gases introduced at a given spot or formed inside the fluid bed by a chemical 
reaction are thoroughly mixed with the fluidized material. Hence, it is an obvious step to utilize 
these properties of the fluid bed for stopping rapid chemical reactions and to examine whether 
unstable intermediate products can thus be isolated. This may be achieved by two methods : 
(1) The reaction starts outside the fluid bed, which acts as quenching medium ; 

(2) The reaction proceeds inside the fluid bed with direct participation of the fluidized material. 


These two principles are explained by means of three practical examples. 


Résumé 
4 l'intérieur du lit sont trés rapidement dissipées. De méme, des gaz introduits en un point donne 
ou formés 4 lintérieur du lit fluidisé par réaction chimique sont intimement mélés au matériau 
fluidisé. Par suite, il est logique d’utiliser ces propriétés du lit pour arréter les réactions chimiques 


Les quantités de chaleur fournies en des points isolés d'un lit fuidisé ou dégagées 


rapides et pour tenter d'isoler s'il y a lieu les produits intermédiaires instables. 


L’auteur propose deux méthodes : 


(1°) 
un inhibiteur. 


Amorces de réaction extérieures au lit fluidisé, lorsque celui-ci se comporte comme 


2°) Développement de processus réactionnels a l'intérieur du lit fluidisé avec participation 


directe du matériau fluidisé. 


Ces deux principes sont illustrés par 3 exemples pratiques. 


1. EINLEITUNG 
Den Begriff der * Wirbelschicht " gibt es seit dem 
Jahre 1921, als Frirz Winker erstmalig die 
“ kochende * Bewegung einer glihenden Fein- 
kohlefillung beobachtete, die von Luft aufwarts 
durchstrémt wurde. mit 
diesem neuartigen hydrodynamischen Zustand 
legte er in den Patenten zur Vergasung von 
feinkérniger Kohle oder Koks nieder, die zu der 
bekannten Entwicklung des Winkler-Generators 
fihrten, aber auch schon die Anwendung auf 
Réstprozesse u. dergl. erwihnten. Bereits 1926 
wurde auch fiir katalytische Crackverfahren dieses 
Prinzip der friheren I. G. Farben-Industrie A. G. 
patenticrt. Der eigentliche grosse Ausbau des 
Wirbelschicht-Crackens begann jedoch erst etwa 
12 Jahre spater in den Vereinigten Staaten, als 
dort die Wirbelschicht unter der Bezeichnung 


Seine Erfahrungen 


Kntschei- 
dend war die kurz darauf gemachte Feststellung, 
dass feinkérnige Stoffe im Wirbelzustand ahnlich 
wie Flissigkeiten durch Rohrleitungen geférdert 
werden und mit Hilfe Standrohren auch 
Druckdifferenzen gehalten werden kénnen [1]. Aus 
dieser Feststellung, die auf den ersten Blick 
vielleicht unbedeutend erscheinen mag, 
entwickelten sich in USA die grossen fluidized 
Crackverfahren mit Kreislauf des Katalysators 
zwischen Reaktor und Regenerator, wie sie unter 
den modernen Einrichtungen der Erddlraflinerien 
auch in Europa errichtet wurden. Inzwischen ist 
“ Fluidization "’— die Wirbelschichttechnik 
mitt Feststoff-Kreislauf oder mit einmaligem 
Feststoff-Durchsatz — auch fiir andere chemische 
Prozesse Bedeutung geworden ; ihre 
Entwicklung ist noch keineswegs abgeschlossen. 


“ fluidized bed" neu entdeckt wurde. 


von 


etwas 


die 


yon 
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Neben der hohen hydrodynamischen Beweglich- 
keit—dem “ Fliessvermégen "’—des Wirbelzu- 
standes liegen seine besonderen Vorteile in der 
grossen Grenzilache zwischen dem feinkérnigen 
Wirbelgut und der Gasstrémung sowie in der 
hohen Intensitat der Warmedurchmischung 
innerhalb der Wirbelschicht. Diese starke 
Warmedurchmischung beruht ebenso wie der 
hohe Warmetibergang an Heiz- und Kuhlelemen 
ten auf dem Warmetransport durch die wirbelnden 
Partikeln, der die Warmekonvektion der Gas- 
stroOmung um Grdéssenordnungen tbersteigt und 
fir sehr schnellen Ausgleich értlicher Tempera- 
turunterschiede sorgt. Auch die Einmischung 
eines etwa seitlich oder von oben eingefihrten 
Gasstrahls in das Wirbelgut erfolgt verhaltnis- 
miassig schnell, doch fehlt hier ein analoger 
Transporteffekt der wirbelnden Partikeln, sodass 
diese stoffliche Durchmischung langsamer vor 
sich geht als z.B. bei turbulenter Strémung im 
leeren Rohr. 

Immerhin legen diese Eigenschaften der Gas- 
wirbelschichten nahe zu versuchen, ob man mit 
ihnen schnell verlaufende mehrstufige chemische 
Reaktionen abbrechen und dabei instabile Zwi- 
schenprodukte erfassen kann. Von drei Versuchen 
dieser Art, die sich durch das Arbeitsprinzip und 
die apparative Anordnung voneinander unter- 
scheiden, soll im folgenden die Rede sein : vom 
Wirbelschicht-Tauchbrenner, von der Wirbel- 
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Ass. 1. Tauchbrenner in der Wirbelschicht, Prinzip der 
Anordnung. 








schicht-Funkenentladung und von der heiss-kalten 
Wirbelschicht. 


Ass. 2. Gesamte Versuchseinrichtung mit Wirbelschicht- 


Tauchbrenner. 


WIRKBELSCHICHT-TAUCHBRENNER 

(GEMEINS. MIT F, Frrvine) [2] 
Der Wirbelschicht-Tauchbrenner arbeitet nach 
dem Prinzip, die schnelle chemische Reaktion als 
Flamme vorgemischter Gase ausserhalb der 
Wirbelschicht anlaufen zu lassen und die heissen 
Flammengase durch die Wirbelschicht schnell 
abzukiihlen bzw. abzuschrecken. Als Modell- 
reaktion wurde die unvollstandige Verbrennung von 
Methan mit Sauerstoff gewahit, die in dem bekann- 
ten Verfahren von Sacusse und BartTHotomé 
durch Abschreckung der Flammengase mit einge- 
spritztem Wasser das instabile Zwischenprodukt 
Acetylen liefert. Als Wirbelgut bewahrte sich 
Siliciumearbid von 50 bis 70 « Korngrésse, als 
Wirbelgas wurde Stickstoff verwendet. Abb. 1 
zeigt schematisch die Anordnung des in die 
Wirbelschicht eingehangten Tauchbrenners (s), 
der hier als Kelchbrenner gezeichnet ist, mit der 
angedeuteten Flammenfront (i). Die Temperatur 
der Wirbelschicht wurde in der Regel durch eine 
Kiihlschlange (k) niedrig gehalten; sie musste 
aber natirlich genigend oberhalb des Taupunktes 
der Rauchgase liegen, um Kondensationen zu 
vermeiden. 

Die gesamte apparative Anordnung ist in 
Abb. 2 dargestellt. Die Wirbelschicht (a) von 
8cm Ruheschichthéhe befand sich in einem 
Rotosiltiegel (6) von 12 cm Innendurchmesser. Als 
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Anstrémboden war eine keramische Filterplatte 
(d) eingekittet. Die Kihlschlange zur Tempera- 
turhaltung der Wirbelschicht auf ca. 100°C ist 
hier nicht eingezeichnet. Fir eine Reihe von 
Versuchen mit héheren Abschrecktemperaturen 
konnte der Wirbelschichtbereich von aussen 
beheizt (¢) und auch der von unten zustrémende 
Stickstoff vorgeheizt werden. Ein Thermoelement 
(t) diente zur Messung der mittleren Wirbel- 
schichttemperatur. Als Tauchbrenner (s) ist hier 
der Typ eines Bunsenbrenners mit zylindrischer 
Misch- und Brennkammer eingezeichnet. Die 
Siliciumearbid-Schicht wurde mit 350 bis 400 em*® 
NTP Stickstoff je Sekunde homogen aufgewirbelt, 
das entspricht einer Lineargeschwindigkeit von 
8,5 bis 4cm,sec unter Normalbedingungen im 
leeren Rohr. 1,6 bis 1/8 dieser Menge wurde an 
Brenngasgemisch dem Tauchbrenner zugefihrt. 
Der Brennerrand tauchte 3 bis 5cm unter die 
Wirbelschichtoberflache. Dies zur 
Abschreckung der Flammengase, brachte aber 


genugte 


noch keine vollig gleichmassige Verteilung der 


Konzentrationen tiber den Schichtquerschnitt. 
Daher wurden in Abzugskanal (g) zwei 
Turbulenzsiebe (m) eingebaut, die sicherstellten, 
dass durch den Saugstutzen (n) zuverlissige 
Analysenproben in die Glaskolben (0) und (p) 
gezogen werden konnten. 


den 





Ass. 3. Kelchbrenner aus Quarzglas (Mage in mm). 


Von verschiedenen untersuchten Brennerarten 
bewahrten sich am besten die Typen des Kelch- 
brenners und des Bunsenbrenners. Abb. 3 zeigt 
einen ganz aus Quarzglas gefertigten Kelchbrenner. 
Die getrennt zugefiihrten Gase Methan und 
Sauerstoff mischten sich in der Kugel (a) und der 
turbulent durchstrémten Kapillare (6). Der 
Flammenraum war als Diffusor ausgebildet ; ein 
eingeschmolzener Quartzstab (d) diente zur Stabi- 
lisierung des Flammensitzes. Die Brenngase 
durchstrémten Flammenzone in 0,02 bis 
0,03 sec; unmittelbar unterhalb der Flammen- 
front setzte die Kihlwirkung der Wirbelschicht 
mit Temperaturgradienten von etwa 400 grd ‘cm 
ein. Nach einer Abkihlzeit von etwa 0,05 sec, 
d.h. ungefahr 2mm vor der Brennermiindung, 
hatten die Rauchgase bereits die mittlere Wirbel- 
schichttemperatur angenommen. Bei einer Reihe 
von Versuchen wurde die Brennermiindung mit 
einem Nickeldrahtnetz (e) (Maschenweite 0,8 mm) 
abgedeckt, um die Flamme vor  stossenden 
Partikelwolken bei heftiger Bewegung der Wirbel- 
Dadurch wurde die Abkihl- 
zone aus der Brennkammer heraus in die Brenner- 
mindung verlegt, doch blieb dies auf 
Abgaszusammensetzung ohne Einfluss. 


die 


schicht zu schiitzen. 


die 


30 


26) 


22 











Ass, 4. Mit dem Kelchbrenner nach Abb. 3 erhaltene 
Acetylenausbeute,n abhangig vom Methangehalt im Brenn- 
gas-Gemisch (O.-CH,-Gemisch). 


Die Acetylenausbeuten dieses Kelchbrenners bei 
einer Wirbelschichttemperatur von ca. 100°C sind 
in Abb. 4 dargestellt, als Funktion des Methange- 
haltes im Brenngasgemisch. Ausbeute bedeutet 
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hier Anteil des Kohlenstoffes, der aus dem um- 
gesetzten Methan in Acetylen iibergegang:n ist. 
Die eingezeichneten Unsicherheitsspannen geben 
die Fehlerstreubreite der Acetylenanalysen an. 
Die Punktgruppe 1 wurde mit frei im Rohr 
haingendem Kelchbrenner, 10cm oberhalb der 
Wirbelschicht, erhalten. Auch dabei treten 
schon erhebliche Acetylenausbeuten auf; durch 
Eintauchen des Brenners in die Wirbelschicht 
werden sie etwa verdoppelt (Punktgruppe 2). 
Mit steigendem Methangehalt im  Brenngas 
steigen die Acetylenausbeuten etwa_ linear 
an, wie man es zu erwarten hat, solange die 
Flammentemperaturen noch nicht unter den 
ginstigen Bereich von etwa 1300°C absinken. 
Ausser dem Acetylen enthalten die Abgase neben 
Stickstoff und Wasserdampf erhebliche Mengen 
Wasserstoff und Kohlenoxyd sowie geringe Anteile 
an Kohlendioxyd und nicht umgesetztem Sauer- 
stoff und Methan. Die héchsten Ausbeuten von 
26 bis 30°, liegen in dem Bereich der Werte, die 
fir den Acetylenbrenner von SacusseE und 
BaRTHOLOME angegeben werden. 














Ass. 5. Bunsenbrenner aus Stahlrohr mit Quarzglasglocke 
(Mae in mm). 


Uber einen Methangehalt von 60%, hinaus 
konnte in dem Kelchbrenner keine stabile Flamme 
aufrechterhalten werden. Dies gelang jedoch mit 
dem in Abb. 5 skizzierten Bunsenbrenner. In 


einem Stahlrohrstutzen (a) wurden die Brenngase 
gemischt und konnten von aussen (h) vorgeheizt 
werden. Den Abschluss zur Brennkammer bildete 
eine Quarzwolleschicht (d) mit einem Drahtnetz (5), 
auf dem die Flamme aufsass. Die Brennkammer 
selbst bestand aus einer tibergeschobenen Quarz- 
glasglocke, deren Miindung wieder ein Nickel- 
drahtnetz (e) abdeckte. Die Flammentemperatur 
von etwa 1300°C reichte bis dicht an dieses Netz 
heran, darauf folgte ein sehr steiler Temperatur- 
abfall von etwa 1800 grd/cm. _ Bereits 4mm 
jenseits der Brennermiindung hatten die Rauch- 
gase nach nur etwa 0,01 sec Abkihizeit die 
mittlere Wirbelschichttemperatur angenommen. 
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Ass. 6. Mit dem Bunsenbrenner nach Abb. 5 erhaltene 
Acetylenausbeuten, abhingigvom Methangehalt im Brenn- 
gas-Gemisch. 


Die mit diesem Brenner erzielten Acetylenaus- 
beuten bei 100°C Wirbelschichttemperatur zeigt 
Abb. 6. Die Werte des Kelchbrenners wurden 
nicht ganz erreicht, obwohl hier der Methangehalt 
im Brenngas bis auf 62,5°%, gesteigert werden 
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Ass. 7. Acetylenausbeuten des Bunsenbrenners bei héheren 
Wirbelschichttemperaturen. 
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Bei Vorheizung auf 650°C brannten 
» Methan noch mit 
stabiler Flamme, doch ergaben auch diese keine 


konnte. 


27 0 


sogar Gemische mit bis zu 67 


héheren Acetylenausbeuten. 

Mit dem Bunsenbrenner wurde 
erhéhter Wirbelschichttemperatur 
Abb. 7 zeigt die Ergebnisse bei Abschrecktem- 
peraturen von 400, 600 und 800°C. Die Acetylen- 
ausbeuten hier erstaunlich 
hoch ; erst bei 800°C fielen sie erheblich ab. Das 
bei 400°C und insbesondere bei 600°C auftretende 
Maximum ist offenbar darauf zuriickzufiihren, dass 
Flammen- 

Bereich 
absinken. Wirbel- 
schichttemperatur wurde keine Andeutung eines 


auch bei 


gearbeitet. 


blieben auch noch 


bei héheren Methangehalten die 


temperaturen unter den optimalen 


In den Versuchen mit 100 
solchen Maximums gefunden; durch giinstige 
Flammenstabilisierung _ sollte 
Acetylenausbeute dort eventuell noch steigern 


man daher die 


kénnen. 
Die Untersuchungen ergaben somit, dass Wir- 
durchaus sind, schnelle 


belschichten geeignet 


Reaktionen abzubrechen und instabile Zwischen- 


produkte zu fassen. Fir die zu diesem Nachweis 


benutzte Modellreaktion, die Acetylenbildung aus 
Methan durch Teilverbrennung mit Sauerstoff, 
ist dieses Verfahren technisch allerdings kaum 
von Bedeutung. Dazu Anteile 
Acetylens im Abgas zu gering (0,3 bis 2,2 Vol.-°,) 
und die Durchsatze im Brennerquerschnitt zu 


sind die des 


klein (einige Prozente der Durchsatze im Accty- 
lenbrenner nach SacussE und Bartuo_ome£). Man 
kiénnte zwar beide Nachteile dadurch vermindern, 
dass man unter erhéhtem Gasdruck arbeitet und 
anstatt Stickstoff als Wirbelgas Wasserdampf 
benutzt. Dieser kénnte bei der Temperatur- 
haltung der Wirbelschicht selbst erzeugt und aus 
dem Abgas unter Warmeriickgewinnung wieder 
herauskondensiert werden. Doch sollte ja kein 
neues Acetylenverfahren entwickelt sondern viel- 
mehr gezeigt werden, dass bei Abschreck vorgangen 
Wirbelschichten mit gutem Erfolg und in breiten 
Temperaturbereichen Verwendung finden kénnen. 


WIRBELSCHICHT-FUNKENENTLADUNG 
(GEMEINSAM MIT 
H. Scumipt vu. M. Buscumanyn) 


8. Die 


Eine zweite Arbeitsweise besteht darin, die zu 


unterbrechende Reaktion innerhalb der Wirbel- 
schicht ablaufen zu lassen, wobei dann auch das 
Wirbelgut selbst an der Umsetzung beteiligt sein 
kann. In diesem Falle sind Temperaturgradienten 
in der Wirbelschicht erforderlich, die bekanntlich 
nur schwer aufrecht erhalten werden kénnen. Bei 
elektrisch leitendem Wirbelgut sind jedoch durch 
Funkenentladung 6rtlich zeitlich 
Temperaturgradienten leicht zu 
erzeugen. Die eine Elektrode kann z.B. als 
Drahtnetz oder Drahtspirale auf den Anstrém- 
boden gelegt, die andere als Spitzen- oder Bi- 
schelektrode in Rohrmitte in geeignetem Abstand 
oberhalb des Anstrémbodens angebracht werden. 
Befindet sie sich etwa 1 cm oberhalb der Schicht- 
oberflache, so geschieht bei ruhender Schicht auch 


und schnell 


wechselnde 


nach Anlegen einer Wechselspannung von einigen 
KV zunichst noch nichts. Erst wenn die Partikeln 
nach Aufwirbelung und Schichtexpansion in die 
Nahe der Spitzenelektrode aufsteigen, setzt ein 
Funkenregen ein, zahlreiche Teilentla- 
dungen, Korn zu Korn die 
Wirbelschicht durchsetzen. Besonders feinteilig 
und gleichmassig aufgespalten ist dieser Funken- 
regen bei feinkérniger, graphitischer Kohle bzw. 
bei Kokspartikeln als Wirbelgut. Bei aufgewir- 
belten Metallpulvern teilt sich die Entladung 
weniger stark auf; die einzelnen Funken werden 


indem 


von springend, 


kraftiger und harter und durchsetzen nur einen 
Teil des Wirbelschichtraumes in der Umgebung 
des nachsten Abstandes zwischen den Elektroden. 
In jedem Fall glihen jedoch die von den Funken 
auf ihrem Zickzackweg getroffenen Partikeln 
momentan hell auf. In dieser Weise entstehen 
die erwahnten steilen Temperaturgradienten in 
den Ortlich und zeitlich schnell 
Einzelentladungen, insbesondere an den Brenn- 
flecken auf den Partikeloberflachen. Dort 
kénnen chemische Umsetzungen  thermisch 
geziindet werden, die anschliessend durch Abkih- 
lung sehr schnell wieder abbrechen. 

Mit feinkérniger Aktivkohle als Wirbelgut 
wurden einige Vorversuche zur Umsetzung mit 
verschiedenen Wirbelgasen in der Funkenent- 
ladung durchgefiihrt. Die Gase strémten mit 
Zimmertemperatur an, die mittlere Temperatur 
der Wirbelschicht lag infolge der Warmeentwick- 
lung durch den Stromfluss von 50 bis 100 mA 


wechselnden 
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zwischen 100 und 200°C. Mit Wasserstoff als 
Wirbelgas wurde Acetylen gebildet, mit Stickstoff 
Dicyan, mit Ammoniak Cyanwaserstoff, Kohlen- 
dioxyd setzte sich zu Kohlenoxyd um. Die 
Ausbeuten dieser endothermen  Reaktionen, 
gerechnet als Verhaltnis der chemisch gebundenen 
Energie zur elektrisch aufgewandten, waren 
natirlich gering; im besten Falle, bei der 
Reduktion des Kohlendioxyds, kamen sie an 10°%, 
heran. 

Gréssere stoffliche Umsatze liessen sich erwar- 
ten, wenn die endothermen Prozesse — ahnlich 
wie bei der Teilverbrennung des Methans — mit 
exothermen gekoppelt sind. Als Modellreaktion 
fiir einen solchen Fall wurde die Vergasung von 
Kohle mit Wasserdampf und Sauerstoff gewahlt 
[3]. Mit steigendem Sauerstoffgehalt im Wirbelgas 
sollte man hier den Bereich von der elektro- 
thermischen Arbeitsweise — mit reinem Wasser- 
dampf sowie geringem Sauerstoffzusatz — bis zum 
autothermischen Verhalten iiberstreichen kénnen, 
bei dem die chemische Warmeerzeugung durch 
die exotherme Verbrennung tberwiegt und 
sich die Temperatur der Wirbelschicht daher 
selbstandig einstellt. 





Abgas 
(He CO,CO,) 


NM 
Thermo- 
elernente 
monte! 
Vor - 


heizung 
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Ans. 8. Wirbelschichtreaktor aus Quarzglas mit Hoch- 
spannungs-Funkenentladung. 


Die fiir diese Untersuchungen benutzte Appa- 
ratur ist in Abb. 8 skizziert. Das Wirbelschicht- 
rohr aus Quarzglas von 40mm _ Durchmesser 
besass eine eingeschmolzene Quarzfilterfritte als 
Anstrémboden ; die beschriebene Elektrodenan- 
ordnung wurde mit 2 bis 8 kV Wechselspannung 
betrieben. Als Wirbelgut wurde durchweg 
Aktivkohle Supersorbon von 0.1 bis 0.4mm 
Korngrésse verwendet und zu 1 bis 8 em Schicht- 
héhe eingefillt. Mit Wasserdampf sowie Wasser- 
dampf-Sauerstoff-Gemischen wurden die Schichten 
aufgewirbelt und auf 3 bis 7 cm expandiert. Es 
wurde mit Elektrodenabstanden zwischen 1,5 und 
5cm, mit Entladungsstrémen zwischen 50 und 
150 mA gearbeitet. Um stérende Kondensationen 
zu vermeiden, wurden die Wirbelgase auf 250 bis 
350°C vorgeheizt. Zur Warmeisolierung war das 
Reaktionsrohr lediglich von einem etwas weiteren, 
der Lange nach in zwei Halften geteilten Glas- 
zylinder umgeben, sodass man die Vorginge in 
der Wirbelschicht gut beobachten konnte. Ther- 
moelemente und Thermometer dienten zur Tem- 
peraturkontrolle. Nach Kondensation des 
umgesetzten Wasserdampfes wurde das Abgas, 
das im wesentlichen aus Wasserstoff, Kohlenoxyd 
und Kohlendioxyd bestand, analysiert bzw. in 
einer kleinen Flamme laufend verbrannt. 

Massgebend fiir das Verhalten der Wirbel- 
schicht und die Gasproduktion war erwartungs- 
gemass der Sauerstoffgehalt im Frischgas. Mit 
wachsendem Sauerstoffgehalt stieg die mittlere 
Temperatur der Wirbelschicht stetig an _ bis 
schliesslich zur sich selbst unterhaltenden Ver- 
brennung der Kohlefiillung. Eine Ubersicht iiber 
den ganzen Bereich der verschiedenen Zustande 
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Ass. 9. Temperaturverlauf in Aktivkohle-Wirbelschichten. 
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gaben Messungen der Temperaturverteilung lings 
der Schichthéhe. Sie wurden mit einem ver- 
schiebbaren Thermoelement durchgefiihrt, das 
in einem an der Spitze besonders diinnwandigen 
Quarzréhrehen steckte. Die Ergebnisse sind in 
Abb. 9 zusammengestellt. Bei Aufwirbelung mit 
reinem Wasserdampf ohne angelegte Hoch- 
spannung (Kurve I) lag die Temperatur an der 
Oberseite des Anstrémbodens bereits etwa 50° 
unter derjenigen des Frischgases, blieb langs der 
Wirbelschicht konstant und begann erst im 
Oberflachenbereich der Schicht abzusinken. Bei 
Zufuhr elektrischer Energie durch Funkenent- 
ladung (Kurve II) lag die Temperatur der Schicht 
tiber der des Frischgases, die Temperaturver- 
teilung blieb jedoch gleichférmig. Wurden nun 
steigende Anteile des Wasserdampfes durch 
Sauerstoff ersetzt, so verschob sich das Tempera- 
zunichst noch gleicher Gestalt 
Von etwa 


turprofil bei 


weiter zu héheren Temperaturen hin. 
15 bis 20°, Sauerstoff ab bildete sich jedoch ein 
Temperaturgefalle von unten nach oben aus, das 
bei 25°, schon recht ausgeprigt war (Kurve III). 
Gleichzeitig begann die Wirbelschicht, sich durch 


die chemische Warmeerzeugung _ selbstandig 
glihend zu erhalten und von der elektrischen 
Energiezufuhr unabhangig zu werden. Oberhalb 
30°,, Sauerstoff war dieser autothermische Zu- 
stand erreicht, mit einer mehr oder weniger aus- 
gepragten Orydationszone am Fusse der Schicht. 
Das héchstliegende Temperaturprofil in Abb. 9, 
Kurve IV, stellt dieses Verhalten etwas iber- 
spitzt dar. Wegen des schnellen Abbrandes bei 
50°, Sauerstoff musste hier eine grébere Kohle- 
kérnung bis 2 mm verwendet werden, bei der die 
grésseren Partikeln sich in der Nahe des Anstrém- 
bodens anreicherten; daher dort die stark 
ausgepragte Temperaturspitze der Oxydations- 
zone. Hier war somit der Zustand eines thermischen 
Wassergas-Generators erreicht, allerdings eines 
schlechten Generators, da wegen des kleinen 
Schichtquerschnittes und der mangelnden War- 
meisolierung hohe Sauerstoffgehalte im Frischgas 
erforderlich waren und die Temperatur unmittel- 
bar oberhalb der Oxydationszone steil absank, so 
dass sich keine wirksame Reduktionszone aus- 
bilden konnte. 

Fir den gesamten Bereich von reinem Wasser- 


dampf bis zu 60°, Sauerstoffgehalt wurden die 
Abgase nach Menge und Zusammensetzung be- 
stimmt. Die Gasproduktion, welche ohne Sauer- 
stoff nur einige em*/sec, d.h. wenige Prozente des 
durchgesetzten Wasserdampfs, betrug, stieg mit 
wachsnedem Sauerstoffzusatz starker als linear 
bis auf etwa 35 cm®/sec an. Einen Uberblick tiber 
die Zusammensetzung der (trockenen) Abgase 
gibt Abb. 10. Mit reinem Wasserdampf entstanden 
bei mittleren Temperaturen der Wirbelschicht um 
500°C etwa 52°, H,, 36% CO, 9%, CO,, 2% CH, 
sowie 1%, O, und Spuren von C,H,. Ausserdem 
bildeten sich héhermolekulare Verbindungen, die 
sich im Kondenswasser listen, insbesondere 
Ameisensaure. Der Methananteil konnte durch 
Verwendung sehr feinteiliger Kohle (mittlerer 
Korndurchmesser 60 x) und Zusatz von Nickel- 
oxyd oder Raney-Nickel als Katalysatoren zur 
Wirbelschicht bis auf 5°, erhéht werden. Die 
gesamte Energieausbeute dieser elektrothermi- 
schen Wasserdampfvergasung bei niedrigen Wir- 
belschichttemperaturen betrug 5 bis 10°. 


. 
onan OQ, 
a 375 5 
% O 
H,O—-O,-Vergasung einer Aktivkohle-Wirbel- 
schicht. 


Ans. 10. 


Bei Sauerstoffzusatz zum Frischgas stieg die 
gesamte Gasproduktion durch verstarkte Bildung 
von Kohlenoxyd und Kohlendioxyd erheblich 
an. Priméar erzeugt der Sauerstoff durch Angriff 
an den heissen Brennflecken der Funkenentladung 
Kohlenoxyd; dieses verbrennt jedoch, be- 
schleunigt durch die katalytische Wirkung der 
anwesenden Wasserstoffverbindungen, zu einem 
erheblichen Teil zu Kohlendioxyd nach; daher 
die steile Zunahme dieser Komponente in Abb. 10 
bei steigendem Sauerstoffzusatz. Die Wasser- 
dampfzersetzung wird durch die lokale Warmeer- 
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zeugung beim Angriff des Sauerstoffs auf die 
Kohle verstarkt ; doch verbrennt ein erheblicher 
Teil des Wasserstoffs wieder zu Wasserdampf. 
Diese Nachverbrennung des Wasserstoffs, zusam- 
men mit dem starken Anstieg der gesamten 
Abgasmenge, verursacht das steile Absinken des 
Wasserstoffgehaltes (wie auch des Methange- 
haltes) im Abgas. Mit weiter steigender Sauer- 
stoffzufuhr und auf 600 bis 800°C wachsenden 
Wirbelschichttemperaturen verbrennen immer 
gréssere Anteile des Kohlenoxyds zu Kohlen- 
dioxyd; ausserdem greift hier der Sauerstoff, 
unabhangig von der Funkenentladung, die ge- 
samte Oberflache der Kohlepartikeln an und 
erzeugt dabei im Primarprozess bereits erhebliche 
Mengen CO, neben CO. Daher fallt der Kohlen- 
oxydanteil im Abgas gegeniiber dem Kohlen- 
dioxyd immer starker ab. Erst bei verhaltnis- 
miassig hohen Saurestoffzusitzen und _selbst- 
glihender Wirbelschicht deutet der Wiederan- 
stieg der Wasserstoffausbeute einen Ubergang zu 
andersartigem Verhalten an. Dies hangt jedoch, 
wie schon an den Temperaturkurven in Abb. 9 
erlautert, mit dem vollstandigen Verbrauch des 
Sauerstoffs in einer verhaltnismassig schmalen 
Oxydationszone und der Ausbildung einer darii- 
berliegenden Reduktionszone, also den tblichen 
Verhaltnissen rein thermischer Gaserzeugung, 
zusammen, 

Die aus diesen Ergebnissen zu ziehenden 
Schliisse kénnen kurz folgendermassen formuliert 
werden : Die Abschreckwirkung einer Kohlewir- 
belschicht reicht aus, um Verbindungen, die sich 
in der beschriebenen Funkenentladung mit 
geeigneten Wirbelgasen bilden, wie z.B. Methan, 
Ameisensaure, Acetylen, Cyanverbindungen usw, 
zum Teil zu niedrigen Temperaturen herunter zu 
retten, bevor sie thermisch wieder zerfallen. Sie 
reicht jedoch nicht aus, um bei Sauerstoff-haltigem 
Wirbelgas die Nachverbrennung von Wasserstoff 
und von Kohlenoxyd in wasserstoffhaltiger 
Atmosphare in starkerem Ausmass zu verhindern. 


4. Hetss-KALTE WIRBELSCHICHT 
(GEMEINSAM MIT M. BuscHMANN) 


Parallel mit den Versuchen an der Funken- 


durchsetzten Wirbelschicht wurde geprift, ob es 
nicht doch médglich sei, in einer Wirbelschicht 


auch ein stationdres Temperaturgefalle zu erzeugen 
und aufrechtzuerhalten. Hierzu wurde das 
Wirbelgas vor Eintritt in die Schicht stark vorge- 
heizt - auf 800 bis 1100°C —und in den Ober- 
flachenbereich der Wirbelschicht eine Kihl- 
schlange eingehingt. Abb. 11 zeigt die apparative 
Anordnung. Das Wirbelschichtgefiass mit Frisch- 
gaszufuhr war wieder aus Quarzglas hergestellt. 











Ass. 11. Heiss-kalte Wirbelschicht. 


In den Frischgasstutzen wurde ein Pythagorasstab 
mit aufgewickelten Heizdrahtwendeln eingebaut. 
Oberhalb dieses Heizelementes erweiterte sich das 
Rohr diffusorartig zum Wirbelschichtraum. Am 
Grunde dieser Erweiterung sass ein Quarzwolle- 
pfropfen als Anstrémboden, 5 bis 8 cm oberhalb 
tauchte die Kupferkihlschlange in den Ober- 
flachenbereich der brodelnden Wirbelschicht ein. 
Zu-und Ablauf des Kihlmittels waren durch eine 
Schliffkappe am Kopf des Apparates gefiihrt ; 
ebenso ein Thermoelement, das langs der Rohr- 
achse verschoben werden konnte, um die Tem- 
peraturverteilung tiber die Schichthéhe zu messen. 
Solche Temperaturprofile zeigt Abb. 12. Sie 
wurden erhalten mit feinkérnigem Calciumnitrid 
(Ca,N,) als Wirbelgut und Stickstoff als Wirbel- 
gas, der auf 800 bzw. 1000°C vorgeheizt worden 
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2 3 a 5 6 

cm uber Anstromboden 
Ass. 12. Temperaturverlauf in heiss-kalter Wirbelschicht. 
war. Quarzsand als Wirbelgut zeigte ahnliche 
Ergebnisse. Bei geringerer Gasstroémung und 


ruhender Schicht ergaben die Messungen im 


Zustand 
Kihlschlange die Vorheiztemperatur des Gases, 


stationaren bis einige cm vor der 
dann folgte ein sehr steiler Temperaturabfall. Mit 
Aufwirbelung Schicht 

Temperaturgefalle 


Quarzwollepfropfens 


beginnender der brach 


dieses __ steile zusammen ; 
unmittelbar 


wurden im Wirbelzustand Temperaturen gemes- 


oberhalb des 
sen, die bereits erheblich unter der Gasvorheizung 
lagen. Trotzdem verblieb in der Wirbelschicht, 
wie Abb. 12 zeigt, noch ein Temperaturabfall von 
Hundert Graden. Wir haben 
Zustand daher eine “ heiss-kalte ’’ Wirbelschicht 


mehreren diesen 


genannt. Sie erscheint geeignet fur Umsetzungen, 


in denen das Wirbelgut in verschiedenen Tem- 
peraturbereichen fiir eine Reaktion genutzt und 
wieder regeneriert wird und eventuell gleichzeitig 
gasférmige Zwischenprodukte im Strom der 
Wirbelgase abgeschreckt werden sollen. 

Als Modellreaktion fir Fall 
wurde die Umsetzung von Calciumnitrid mit 
Wasserstoff gewahit, in der ein Teil des Nitrids 


in Caletumhydrid umgewandelt wird : 


Ca,N, + 8H, 7 8CaH, + N, 


einen solchen 


Bei etwas héheren Temperaturen lasst sich das 
Hydrid durch Stickstoff wieder in das Nitrid 
zurickverwandeln. Bei der Hydrierungsreaktion 
tritt intermediar Ammoniak auf, d.h. sie verlauft 
zum Teil tiber die beiden Stufen : 


Ca,N, + 6H, > 8CaH, + 2NH, 
+ 8CaH, + N, + 3H, 


Diese Verhaltnisse wurden bereits im Jahre 1905 
von Haber aufgeklairt [4], bei seinen ersten 
Versuchen zur Ammoniakgewinnung aus den 
Elementen. Beim Uberleiten von Wasserstoff iiber 
Calciumnitrid in einem Porzellanschiffchen stellte 
600°C Ammoniak im 
Gasstrom fest, doch waren die Prozentsaitze so 


er oberhalb abziehenden 
gering —entsprechend der ungiinstigen Lage des 
Ammoniakgleichgewichtes bei diesen Tempera- 
turen 

Es wurde nun frisch hergestelltes Calciumnitrid 
in der beschriebenen heiss-kalten Wirbelschicht 


dass er diesen Weg nicht weiter verfolgte. 


mit einem Stickstoff-Wasserstoff-Gemisch aufge- 
wirbelt, das auf etwa 1000°C vorgeheizt war, und 
das Abgas nach Staubabscheidung auf Ammoniak 
Die 


Ammoniakgehalte lagen zwischen 5 und 10°,, mit 


analysiert. Abb. 13 zeigt die Ergebnisse. 


Gasdurchsatz 
Damit war das Ziel dieser Modell- 


steigendem erwartungsgemiass 


abnehmend. 


"le NH; m Abgas 





WO 200 300 400 
enr¥sec H>-No- Gemisch 
Ans. 13. Ammoniakausbeuten in einer heiss-kalten Wirbel- 
schicht mit Calciumnitrid als Wirbelgut. 


versuche erreicht: mit der Arbeitsweise der 
heiss-kalten Wirbelschicht lassen sich Zwischen- 
produkte, die bei héheren Temperaturen durch 
Reaktion mit dem Wirbelgut intermediar auf- 
treten — wie im vorliegenden Falle das Ammoniak 

durch schnelle Abkihlung zu einem Teil vor 
dem thermischen Zerfall bewahren und auf nor- 
male Temperaturen herunterretten. 

Auch diese Art der Ammoniakdarstellung geht 
jedoch Uber den Charakter einer Modellreaktion 
nicht hinausund hat insbesondere keine Bedeutung 
fiir eine technische Ammoniakgewinnung. Abge- 
sehen von dem Warmeaufwand fiir die Gasvor- 
heizung halt die Wirksamkeit des Calciumnitrids 
nicht lange genug an. Die wechselseitigen Umwand- 


lungen in Hybrid und Nitrid fihren zusammen mit 
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der heftigen Wirbelbewegung der Partikeln zu 
starkem Abrieb; es bildet sich viel Flugstaub, 
der ausgetragen wird. Ausserdem nimmt das 
Ausmass der Nitrid-Hydrid-Umwandlung allmah- 
lich ab. Wie Abb. 14 zeigt, tritt bei der erstma- 
ligen Hydrierung einer frisch hergestellten 
Calciumnitrid-Probe oberhalb 300°C ein kraftiger 
Wasserstoffverbrauch ein. Es ist hier der Gas- 
druck in einer Kreislaufapparatur aufgetragen, in 








i i i s 
300 400 500 600 70 
°C 
Apps. 14. Wasserstoff-Umsatz mit Calciumnitrid. 
Kurve I: erste H,-Behandlung, 


Kurve II: Druckabfall in Hg-Atmosphire nach mehrma- 
liger H,- und N,-Behanfdlung. 


der das Calciumnitrid unter Wasserstoff langsam 


auf die langs der Abszisse angegebenen Tempera- 


turen aufgeheizt wurde. Nach Riickbildung des 
Nitrids durch eine analoge Stickstoffbehandlung 
wurde bei der nachfolgenden Hydrierung weniger 
Wasserstoff umgesetzt, und nach mehrmaligem 
Wechsel zwischen Hydrierung und Nitrierung sank 
der Wasserstoffumsatz auf den durch Kurve II in 
Abb 14 dargestellten geringen Wert. Das Calcium- 
nitrid hat daher fir die Ammoniakgewinnung 
keine geniigende Dauerwirkung. 

Die voranstehend beschriebenen Untersuchun- 
gen haben somit zu keinem direkt anwendbaren 
technischen Verfahren gefiihrt. Sie haben jedoch 
gezeigt, dass Reaktionsgemische mit Wirbel- 
schichten nach verschiedenen Methoden wirksam 
abgeschreckt und dabei instabile Zwischenpro- 
dukte gefasst werden kénnen. 

Die Arbeiten wurden apparativ und finanziell 
unterstutzt durch die Deutsche Forschungsgemein- 
schaft und den Fonds der Chemischen Industrie, 
denen ich an dieser Stelle meinen besonderen 
Dank zum Ausdruck bringen méchte. 
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Abstract—-Approximate laminar boundary layer solutions for mass transfer across the plane 
interface between two co-current parallel fluid streams are derived in the form 
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where p and q are not constants but are functions of U,/U, and pgpg/p)). 


It is suggested that the boundary layer theory 


provides a more realistic physical picture than 


either the Lewis and Warrman stagnant film theory or the Hicsre Penetration Theory. 


Résumé 


L’auteur présente des équations approchées basées sur la théorie des couches limites 


a& écoulement laminaire exprimant I'échange massique a travers l'interface plane de deux co- 


courants fluides, séus la forme 
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Cette théorie donne une image physique plus fidéle que les théories de Lewis et Warrman 
du film stagnant ou théorie de la pénétration de Hicsre. 


1.1. INTRODUCTION 


In 1912 Lanomurr [1] suggested that the resis- 
ance to heat transfer between a solid body and 
neighbouring air could be taken to be that of a 
thin layer of stagnant air in contact with the hot 
body, so that the heat lost was conducted through 
the thin stagnant layer of air. This model for 
convective heat transfer was taken up by Lewis 
and Wurman [2], [3], [4], [5] and applied to 
convective mass transfer. 


1.2. Lewis anpD WuitMan THEORY 


In considering the heat transfer between a drop 
of water and a hot gas, WuitmMan wrote [3], “ It 
is well known that the layer of gas against any 
solid or liquid surface is almost at rest compared 


to the mass of gas outside, and it is this so-called 
film through which heat must pass by the slow 
process of diffusion ... Similarly, a liquid film 
exists on the surface of the water-drop, and the 
combination of gas and liquid films in series gives 
a two-film resistance ... The rate of absorption 
of matter from gas to liquid . . . is a process 
paralleling that of heat transfer.” (This analysis 
into two resistances in series could easily be 
modified to incorporate an interfacial resistance 
if this should prove necessary.) The concept of 
the stagnant film, while valuable in its day, has 
now been replaced in heat transfer studies by the 
concept of the boundary layer. It is the object 
of this paper to show how the boundary layer 
concept can be applied to mass transfer between 
fluid phases. 
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1.3. Hicapie PENETRATION THEORY 


The model of unsteady state diffusion into a 
stagnant fluid was proposed by Hicsie [6] for 
the prediction of mass-transfer coefficients for 
falling liquid films. The theory predicts that at 
time t, the mass transfer coefficient k is given by 
the following equation 


(2) " 


Also, if two fluids are in contact and a solute 
diffusing from one to the other — then 


“ne (2): ail (2) (2) 


nt nt 


It follows from the fact that k,/k, is independent 
of time that the interfacial concentrations (a,), 
and (ay), are also independent of time and 
therefore constant throughout the period of 
diffusion. 


1.4. DANCKWERTS’ MODIFICATION OF 
THE PENETRATION THEORY 


The Hicsre Penetration Theory was revived in 
1951 by Danckwerts [7] but in a modified form. 
In a liquid maintained in turbulent motion by 
stirring at a steady rate, DANCKWERTS supposed 
that the motion of the liquid would continually 
replace with fresh surface those parts of the 
surface which have been exposed for a finite 
length of time. Then, if the mean rate of pro- 
duction of fresh surface (per unit area of surface) 
is constant and equal to f, and the chance of an 
element of surface being replaced within a given 
time is assumed to be independent of its age, and 
if the rate of absorption into fresh surface is given 
by the Hicsre relationship, then 


k= (Df) (3) 


The fact that k is proportional to the square root 
of the diffusion coefficient is not significant since 
it is in essence assumed in the derivation. Further- 
more, the theory is valueless if no method for 
predicting f can be found. However f could not 
be predicted without a detailed knowledge of the 
fluid dynamics of the system, and if a detailed 
knowledge of the fluid dynamics of the system 
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were available then the Danckwenrrts’ theory 
would not be needed. 


1.5. Gas-Fitm—BounDARY LAYER 
THEORY 


Due to the similarity of heat and mass transfer 
[8] it is always possible to write down a mass- 
transfer relationship if the corresponding heat 
transfer relation is available and if the mass 
transfer is small, i.e. at low concentrations. The 
heat transfer between fluid and a flat plate has 
been studied both theoretically [9], [10] and 
experimentally [11], [12], the agreement between 
theory and experiment being satisfactory. 

For laminar flow in the boundary layer 


Nu = 0-331 Pr’ Re} (4) 


In a corresponding mass transfer problem, say 
the evaporation of mercury from an amalgamated 
surface into an air stream, 


+ /Uor\4 
kr _ 9.381 (4 (=) (5) 
D oD) \u 


The experimental results of Maxweit and 
Srorrow [13] for evaporation of mercury confirm 
the above equation, which is a convenient approxi- 
mation for mass transfer when a gas is flowing 
over a liquid. It can only be an approximation 
since the liquid will be set in motion under the 
tangential stress of the gas. 

If a liquid is flowing over another then the 
shear stresses set up by their relative motion 
will determine a velocity distribution which will 
in its turn determine the rate of mass transfer. It 
is this prob!em which is to be discussed in what 


follows. 


2.1. Bounpary Layer THEORY 


When a fluid flows past a flat plate placed parallel 
to the direction of flow of the fluid, the influence 
of the viscosity of the fluid is confined to a thin 
layer called the boundary layer (see Fig. 1). In 
this thin layer the velocity of the fluid increases 
from zero at the plate to the full velocity of the 
fluid. Pranvti [14] therefore suggested that the 
field of flow could be divided for the purposes of 
mathematical analysis, into two regions: the 
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thin boundary layer near the wall in which friction 
must be taken into account, and the region out- 
side the boundary layer where the forces due to 
friction are small and may be neglected, and 
where therefore the classical theory of ideal fluids 
offers a very good approximation. 

With this simplification Pranpt.L was able to 
replace the Navrer-Sroxes differential equations 
of fluid motion — which are difficult to solve — by 
the simpler equation (two-dimensional case) 




















Fic. 1. Boundary layer on a flat plate and velocity distri- 
bution curves in the fluid ahead of and at a point on the 
plate. 


combined with the equation of continuity 


yv 
v 
= 0 


~ (7) 


AF 


For steady flow, equation (6) simplifies to 


(8) 


Furthermore, in the case of a fluid flowing past 
a flat plate the pressure gradient is zero and 
hence equation [8] may be further simplified to 
du du yu 


20 OP em (9) 


“w— 
d y dy 


2.2. Sim1LarRiIty oF VELocITY PROFILES 


An important simplification can be introduced 
into the equation of motion if a “ similarity ” 
solution can be found. “ Similar” solutions are 
those in which two velocity profiles at different 
co-ordinates differ only by a scale factor in u and 
y. The velocity is made dimensionless by making 
the stream velocity U the scale factor, and 
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distance normal to the plate is made dimension- 
less by making the boundary layer thickness the 
scale factor. The partial differential equation (9) 
can then be reduced to an ordinary differential 
equation. For the case of two liquids in co-current 
flow, KeuLeGcan [15] and Lock [16] have shown 
the existence of a similarity solution. A corollary 
of the similarity of the velocity profiles is that 
the velocity at the interface is constant. 


2.8. Von KAaARMAN MOMENTUM 


EQUATION 

It is convenient to regard the boundary layer as 
having a definite thickness 5 so that at the plate 
where y= 0, u=0 and at the edge of the 
boundary layer where y = 5, u = U, the velocity 
of the fluid stream. In actual fact the velocity u 
will approach U asymptotically but the error 
introduced is small [17]. 

If the equation of motion (9) is integrated with 
respect to y between the limits 0 and 4, 


dy + (wi), — | 


0 0 


AY 


du 
u 


From the equation of continuity (7) 


which is equivalent to 


. 8 

JF? > 

on dy — - 

2 y v2 [ud 
0 0 
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“dy 


ou 
0 
dys e 
Hence 
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d , du 
[ (U —u)udy v (*“) 
or dY/o 


0 
and if 


5 
: | o(U —u)udy (10) 
Mr 


. * 


This is the von KaRMAN momentum equation 
(18). 

Approximate solutions for the flow of a fluid 
over a flat plate can be obtained by inserting 
into [14] various approximate expressions for the 
velocity distribution. 


8.1. Co-cURRENT PARALLEL FLUID 
StreaMS—Lock’s APPROXIMATE 
SOLUTION 


Two parallel co-current streams, the upper of 
uniform velocity U, and the lower of uniform 























Fic. 2. Hydrodynamic boundary layers at the interface of 
co-current streams. 


velocity U, are brought together at O (see Fig. 2). 
It is assumed that the interface remains plane. 
The a-axis is the interface and O is the origin. 
The y-axis .is taken normal to the interface. Let 
the velocity of either phase parallel to the x-axis 
be u; ma, and py the viscosities of the upper and 
lower phases, p, and p, their densities. Due to the 
shear stresses arising from the difference between 
U, and U,, two laminar boundary layers are 
formed wherein the velocity changes from the 
interfacial value uy to the stream values U, and 
U,. Let 8, and 4, be the thicknesses of the 
respective layers. The shear stress at the inter- 
face is ty and »,, vg the kinematic viscosities of 
the two phases. 

The considerations which led to the derivation 
of equation (10) for the flow of a fluid over a flat 
plate apply equally in the present case only the 
boundary conditions are different. 

Thus, for the upper and lower streams 


ujudy (11) 


| (U,—u)udy (12) 


. 


0 


= Up 
po du! o 


It then remains to choose suitable velocity 
distributions to insert in equations (11) and (12). 
KEULEGAN [15] has utilized quartic and sextic 
polynomials for the special case of the lower fluid 
at rest (U, = 0), Lock [16] has utilized a quartic 
polynomial for the general case (U, +0) and 
Porrer [19] has investigated the use of a sextic 
polynomial for the general case. Here the quartic 
polynomial will be used. Thus for the upper and 


0 
)+(5) Jos 


8.2. BouNDARY CONDITIONS 


lower streams respectively 


where r = 


Some of the boundary conditions can now be 
written down 
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r (15) 


— 8,, u= U,, 





dy 


These boundary conditions merely assert that 
the velocity must equal u, at the interface and 
U, or U, at the edge of the respective boundary 
layer, and that the velocity profile must pass over 
to the stream velocity without discontinuity at 
the “ edge ” of the respective boundary layers. 

The shear stress should be continuous at the 
interface, so that 


(>) - = 
m dy/ 0+ “e = ry 


If the conditions at the interface are inserted 
into the equation of motion (9) it reduces to 


Be (fe 
dy? 0+ ° dy? 7) 
du 


= 0 and v 
z 


(16) 


(17) 


since 0 at y= 0. It should be 


noted that if mass transfer is occurring, the latter 
statement v = 0 at y = 0 will only be approxi- 
mately true, but the more true as the solutions 
are more dilute and the rate of mass transfer 
correspondingly less. 

These boundary conditions are met by the 
assumed velocity distributions (18) and (14). 


SoL_uTion In Terms or ASSUMED 
VeLocity PROFILES 

If u from equations (18) and (14) is substituted in 

the momentum equations (11) and (12) and 

integration performed with respect to y 


du - d +2 
n(x), male aa 


(3). 


2 [UA —8,MA—r\742-+ 1157) /680) (19) 


r)(74+115r) /630] (18) 


= U, (1 — 7) 2/8, (20) 
O+ 


(“). =U,(A—r)2(—&) — (21) 
therefore on integrating equations (18) and (19) 
with respect to z, and rearranging, 


5° = (4v, 2/U,){630/(74 4+ 115r)} (22) 


8? = (4v,2/U,) {680/(74A + 115r)} (28) 


Hence 
744 115r vg 


h 24 
74A + 115r yy (24) 


52 
52 = 


From equation (16), utilising (20) and (21) 
by/d + A He/ Hy 


Ps 


25) 
5,/3, T Pe/ Fi ( 


Equation (24) becomes 


(8)/()- 
w/e + A 


(26) 


(2) / (2) + 115) + 189, "98 


Equation (26) can be solved by successive approxi- 
mation for any A and jy py/,p,. Then 


du 
To = fy | ), 


== p, U, (1 — r) 2/8, 
whence 


lm? (Tene) (—*) " 


$30 (27) 


To = 
p,U;* i 


ma 
= ¢€ j 
“1 


where c is the shear stress coefficient defined by 
equations (27) and (28). 

Note that when the interface is at rest, i.e. 
r = 0, then c = 0.343. This compares well with 
the exact solution for flow over a flat plate where 
¢ = 0-881 [20]. 


"1 


(28) 
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In Fig. 83 values of c for A= 0 and A= 0-1 
are plotted against p.p./p,p, and in Fig. 4 the 
velocity profiles for A = 10 at various values of 
HgP2/!)p, are presented in terms of the dimension- 
less parameters u/U, and », where 


. = (U,/2v,2)*y, 2 = (U,/2vgr)ty 
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Fic. 3. Shear stress coefficients vs. zgp9/~,p, for U,/U, = 
01 and U, = 0. 


The values shown in Figs. 3 and 4 were actually 
computed from the sextic polynomial solutions. 


4.1 Two-RESISTANCE THEORY APPLIED 
To MOMENTUM TRANSFER 


For mass-transfer between two fluid phases 
Mass 
Anca xTime ~*1l41—(t)o} = ka { (4a), — 4a} 
= K,{A, —(Aj).} (29) 


where k, and k, are the gas and liquid film coeffi- 
cients and K, is an overall coefficient. 


If 
(A,), —_ m Ag (30) 


(81) 
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Fic. 4. Dimensionless velocity profiles for U,/U, = 10, 
and values of pgpo/u\p; from 0-01 to 100-0. 
Symbol MePe/#P1- 

A-001; B-010; C-03816; D-—1-00; 
E-3-16; F- 10-00; G - 100-0. 


Similar relations to these can be written down 
for momentum transfer between two fluid phases 
in relative motion. 

Momentum 
~ Area x Time 
= kg*(pytop2U2) = 
where k,* and k,* are the respective momentum 


transfer film coefficients and K,* is an overall 
momentum transfer coefficient. It will be observed 


1* (p,U, —pyto) 
K,* (p,U, — p, Ug) 


"7 p,U4p,U, 
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Also, 
1 
K,* 
where 
m* = p,/ Py 


Substituting for 7+, from equation (31) 
74 + ~) (Ss) 


yee)" 


*. jr al . 
Ke = 0, (SF) (3 
l-—r 
“Pe Fe Bi Beek. 
be =U, | 


It will be noted that while a particular 
solution of the flow problem depends on A and 
HePs/H,P;, Values of k,* must additionally incor- 
porate the term p,/p, which cancels out however 
when m* /k,* is evaluated. 

From Locx’s solution for air flowing over still 
water at 10°C, taking p,/p, = 0-001247 


k,*/kg* = 17-29 ; k,*/K,* = 1-0215 


As m* is very small, this case in momentum 
transfer is the analogue of the case of a very 
soluble gas in gas-absorption and it can be seen 
that the gas-film momentum transfer coefficient 
equals the overall momentum transfer coefficient 
with an error of two per cent. 


BounDARY Layer THEORY 
FoR Mass-TRANSFER 


5.1. 


The differential equation for molecular diffusion 
in a fluid in two-dimensional motion under steady 
conditions and with insignificant density changes, 
is 


(82) 


where a is concentration and D the diffusion 
coefficient. 


5.2. SIMILARITY OF CONCENTRATION 
PROFILES 
Just as before (Section 2.2) an important simplifi- 
cation can be introduced if the concentration 


profiles are similar at different co-ordinates 2. 
This will be so if the concentration can be expressed 
in such a way as to follow the same boundary 
conditions as the velocity ratio u/U in the viscous 
boundary layer problem, and this can be done if 
a/A be taken to represent the concentration. 
Comparison of equations (9) and (82) establishes 
that the velocity and concentration distributions 
will be identical provided that » = D, i.e. pro- 
vided that the Scumipt number equals unity. It 
follows that if the velocity distributions are 
similar at different co-ordinates z then the con- 
centration distributions will be similar. As a 
corollary the concentrations at the interface will 
be independent of the co-ordinate 2. 


5.8. Co-CurRRENT FLuID STREAMS, 
AuTHOR’S APPROXIMATE SOLUTION 


Suppose now that two streams are flowing co- 
currently as in Fig. 2 and that a solute is being 
transferred from one phase to the other by 
diffusion. It is assumed that the solutions are 
dilute and identical in density and viscosity with 
the pure liquids. If the solutions are dilute then 
the mass-transfer will be small. Modifications 
must be made in the theory if the concentrations 
at the interface and in the stream are widely 
different. This point will be discussed later. 

On each side of the interface there will now be 
two boundary layers, the viscous boundary layer 
of thickness 5, (or — 5) and the diffusion boundary 
layer of thickness A, (or — A,). The thickness of 
the diffusion boundary layer is the distance 
measured normally from the interface in which 
the concentration changes from the interfacial 
value to the stream value. 

If equation (82) is integrated with respect to y 
from y = 0 + to y = h where h is a length which 
is greater than the thickness of either the viscous 
or diffusion boundary layers 


da 
1) ty — — D (=) 
=) 7 . dy 0+ 


since a vanishes for y = h. The subscript 1 has 


y 
been introduced to refer to conditions in the 
upper fluid. 


(38) 
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a 


da, d 

—! dy = — (a, — A,)d 

at hand Jue i) dy 

0 0 

where A, is the concentration outside the 


boundary layer) 


~ Ef lo 


da, d 
— dy = | v — (a, — A,)d 
v dy y J oy‘ 1 1) dy 
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- [v (a, — 4} J (a, — Aj) oy 


Also v vanishes for y= 0 and a, = A, for 
y = h so that the terms with square brackets is 
equal to zero. 

Hence, substituting in equation (33) and making 
use of the equation of continuity (7) 


A 


d as da, 
se | ma eo (>).. 
0 


The limit h can now be replaced by A, since 
outside A, the term (a, — A,), disappears. 


4 


d 7” 1 
2 u(A, — a,) dy = D, () (34) 


Similarly, for the lower stream 
dag 
— ay) dy = Dy =) (35) 
oY Jo- 


Equations (34) and (85) correspond to the 
momentum equations (11) and (12). 


BouNDARY CONDITIONS FOR 
Mass TRANSFER 


5.4. 


The concentrations must equal the interfacial 
values (a), and (dg)y at the interface, and the 
stream values A, and A, at the “edge” of the 
concentration boundary layer, and pass over 


without discontinuity from their values inside the 
boundary layer to the stream concentration. Thus 
y=0-+, a, = (d)o 
y=, a =A), 
g=uo—, 


— A 


ay = (4)o ‘ 
Ody _g, BH _ 
dy dy? 
It will be noted that it has been taken that the 
concentrations remain constant along the length 
of the interface. This point has been discussed 
previously. 

An additional boundary condition is that the 
rates of mass transfer immediately on each side 
of the interface should be equal, i.e. 


» (F),.- (%), 
dy /o+ oy 


since (a,)) is 


»= a, = Ag, 


(37) 


At the constant 


da, 


interface, 
: = 0 at y = 0 and alsov = 0 aty = 0. It has 
OF 

already been pointed out that this will be true 
only if the rate of mass transfer is low, i.e. with 
dilute solutions. It follows that 


(=),. ain (53) 


5.5. CONCENTRATION PROFILES 


Quartic polynomials which meet the boundary 
conditions are 


upper stream : 


lower stream : 


“o-o(-) a 
where i, = (@,)9/Ao; ty = (Gq)o/A 


SOLUTION IN TERMS OF 
ASSUMED PROFILES 


5.6. 


On substituting the velocity distributions from 
equations (13) and (14) and the concentration 
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distributions from equations (38) and (39) in 
equations (84) and (35) respectively, and inte- 
grating with respect to y, assuming A, < 5, and 
A, < 4,, i.e. Scumipr numbers greater than unity 
as will be the case with liquids. 

Ai 


. 


1, = | (A, —a,) udy = A, U, (1 — i) 8, » 


0 
3 j2 2 8 44 L s\| 
pra tan lis” 3+ sap %} | 


1, =| (A,—a,) udy = A,U, (1—i,) (— 44) » 
0 


3 4 1 s| 
_ + xo 41 
140? 3%} (41) 
where o, = A,/3,; o2 = — A, 5. 
For the upper stream it follows that 


(I,)- b, (“*) ~ D, 2A, (1—i,)/, 


oy o-+ 


) 
a 
If equation (22) is differentiated with respect 

to z then 

_ 2x 680», 
‘dz 744+ 115r U, 

This is substituted in equation (42) which can 

then be integrated since o, will be independent 


630 o,* 
74+ 115r 


{3 2 
lio "+" r) |—o,-— 


Similarly for the lower stream 


630 o, 
74A4 115r 
{3 2 3 


 -4-(l— Batic a tee 
110° T°") 15 2 ian? 


=—-=N,e' (44) 
Ys 


From these equations it is possible to calculate 
o, and o, and hence A, and Ay. 
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5.7. Mass TRANSFER COEFFICIENTS 


Local values of the mass transfer coefficient can 
then be calculated. Thus, for the upper stream 


° A — ) 
k, A\(l—%) = — D, (+)... at ae ae 


k, = 2D,/o,8, 


where from equation (22) 


i 
— (as ay) (2) 


2 (* 7 aks 9 (2#)' 


= “ (45) 


Similarly, for the lower stream 


kyr 1 (= + a | (Y=) 


. 630 


= 46 
D . oy ( ) 


“2 

For the case of mass-transfer from a flat plate 
to a parallel stream, u = 0 and r = 0. In this 
case equation (43) reduces to 

2 3 1 - 
(5a —_ ao + w= Nz; (47) 

A solution for (47) which is accurate to within 

5 per cent is [21] 


o, = Ng,’ 


Equation (45) then reduces to equation (5) 
with the difference that the constant is 0-848 and 
not the exact value of 0-331. 

If r = 1 then the upper liquid is moving at 
uniform velocity. On substituting r = 1 into 
equation (43), there remains 

o00 o_ c—- Ns,’ 
74+ 115 10 


5 = Ns.* 


If this value of o, is substituted in (45) and r 
put equal to 1 


= "a (i) (TE) 
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Since the time of contact t = U,/a 
| 4 
ky = 0-548 (7) 


This result compares well with the exact result 
for stagnant diffusion, equation (1) 


; j 
k, = (2) + 0-564 (7) 


mt 


(It might at this point be suggested that the 
approximate method given here would be a 
useful method for introducing unsteady state 
heat conduction or diffusion to students deficient 
in mathematics.) 


5.8. VARIATION IN EXPONENT OF THE 
Scumipt NUMBER 


It has been shown that for flow where the 
interfacial velocity and hence r is zero the mass 
transfer coefficient is proportional to the Scumipt 
number raised to the power of one-third, and 
also that for stagnant diffusion where r = 1, the 
mass transfer coefficient is proportional to the 
ScumMipt number raised to the power of one-half. 
The question arises as to what the exponent 
should be in those cases where r takes a value, 
between 0 and 1, or greater than 1. 

If suitable values for r are inserted in equation 
(48) then solutions can be found for the relations 
between o, and Sc,. In liquid systems the 
Scumipt number will be large and the thickness 
of the concentration boundary a small fraction 
of the viscous boundary layer. It is then possible 
to neglect the last two terms on the left hand side 
of equation (43). 

To determine the exponent on the Scumipt 
number, a solution has been assumed of the form 


oO; = Se,~* 


and a solution found for a range of values of r 
between 0 and 1, assuming Sc, = 100. These 
solutions will not of course be exact but only 
accurate to about 5 per cent over the range of 
Scumipt numbers. The asymptote when r> 1 
can be found by assigning a large value to r, 
say 100,000. Then equation (43) becomes 
approximately 


630 o? 8 
——— ——— | —- X 100,000 — 100,000 x 
115 = 100,000 (10 

4 o7} — o7} + 
15 140 
This simplifies to 
189 
_ o; 
115 


630 2|2 38 3 1 4 of 
o — CO, - o, + go = NN 49 
15 | Fe ‘140 * * 180 ‘| sn (49) 


If Sc, = 1, then it is necessary that o, = 1. 
This can be seen to be true for equation (43). 
If the solution to equation (49) is taken as 


= \ n 
= Se, 
and a value of n sought for Sc, = 100, then 
n 0-55. 
In Fig. 5 values of the exponent n are plotted 
against the relative interfacial velocity r. 
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Fic. 5. Exponent on Schmidt number vs. relative inter- 
facial velocity, r. 





5.9. LIMITATION OF SOLUTIONS 
to Se> 1 


The solutions obtained are limited to the case 
when the diffusion boundary layer is not greater 
in thickness than the viscous boundary layer. 
This will be the case with liquids. When the 
diffusion boundary layer is thicker than the viscous 
boundary layer the algebra is more complicated 
because two different expresionss have to be used 
for the velocity. Squrre [21] states that the 
solution for heat transfer between a fluid and a 
flat plate is unaltered provided that the PranpDTL 
number is greater than 0-5, so it would therefore 
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be correct for the present solutions to say that 
they are accurate to Scumipt numbers greater 
than 0-5. 


LIMITATION TO Low RaTEs oF 


Mass TRANSFER 


6.1. 


When a fluid flows past a flat plate it is obvious 
that the transverse velocity v must be zero at the 
wall if no mass transfer is occurring. If however 
there is mass-transfer then the intrusion of this 
mass into the flowing fluid alters the velocity pro- 
file from what it would have been without mass 
transfer. If the mass transferred is small, then 
the disturbance is small and can be ignored as 
in the present theory. If the mass transfer rates 
are large then a correction must be applied. The 
appearance or disappearance of mass at the 
wall makes the flow problem analogous to that of 
“ blowing ” or “ sucking In 
aerodynamics one of the methods [14] proposed 
for controlling the flow of gas over the foil is to 
suck the gas through the (porous) foil at a con- 
trolled rate. In this case then the transverse 
velocity v is not zero at the wall, ru» +0. The 
and 


> 


* boundary layers. 


hydrodynamics solution of ScHLicuTinc 
BussMANN [22] were utilized by Mickey et al. 
[23] to provide exact solutions for mass transfer 
between a fluid and a flat plate. Eckert and 
LieBLEIN [24] have solved the problem for the 
specific case of water vaporisation. SPALDING [25] 
has utilized the momentum integral approximate 
method for the general case. The solutions so 
obtained show that the rate of mass transfer is 
an important parameter. 


APPLICABILITY OF PENETRATION 
THEORY TO Co-cURRENT FLOW 


6.2. 


If two liquids are flowing co-currently at different 
velocities, the interface will move at a constant 
velocity. In the case of liquids the Scumipt 
numbers are high and so the thickness of con- 
centration boundary layers will only be a fraction 
of the thickness of the viscous boundary layer 
and the variation in velocity over the concentra- 
tion boundary layers small. This leads to the 
conclusion that unless the viscosities, densities 
or velocities of the two phases are widely different 


PoTTEerR 


it will be sufficiently accurate to calculate the 


mass transfer coefficients on the basis of the 
Penetration theory, if the time of contact is 
obtained from the interfacial velocity and the 
length of the interface. This conclusion would be 
invalid for countercurrent flow but no solution 


has yet been obtained for this condition. 


7.1. Discussion 


Interfacial instability 


The question of the stability of the interface has 
not been raised in the present study but the 
solutions given are only applicable as long as the 
interface is in fact stable. Interface (or free 
surface) stability has been investigated mathe- 
matically by Lessen [26], Yin [27] and Cima- 
RULLI [28]. It has been shown that instability 
sets in at low Reynolds number. Yin, for a case 
of vertical flow of a fluid with a free surface, 
demonstrated mathematically that instability sets 
in at Reynolds number of 6, a theoretical value 
which compares well with experimental values of 
the order 10 to 25. Naturally at such low Rey- 
nolds numbers instability means only the onset 
of waves at the free surface and not the onset of 
turbulence. It seems clear however that attention 
must turn to the fundamental study of unstable 
interfaces including the effect of waves on the 
rate and mechanism of mass transfer. In this 
connection the recent study by Jackson [29] 
should be noted. He measured the thickness of 
falling films and concluded that liquids with 
viscosities greater than that of water give values 
of the film thickness which are less than for true 
viscous flow and that the wave motion occurs 
when the Froude number exceeds unity. A 
surprising result was obtained by Srirpa and 
Hurt [30] who measured the rates of solution in 
falling liquid films of slightly soluble solids coated 
on the wall of a vertical tube. They found the 
rate of solution in apparently laminar films to be 
up to two and a half times the theoretical pre- 
diction. On the other hand, Davipson and 
CULLEN [81] have clearly demonstrated that 
absorption of a pure gas into a ripple-free liquid 
falling over a spherical surface can be predicted 
from the laws of flow of a fluid in laminar flow. 
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7.2. TurRBULENT Bounpary LAYERS 


The present theory cannot be extended easily to 
deal with the problem of turbulent boundary 
layers. A more extensive knowledge is first 
required of the decay of turbulence near liquid/- 
liquid and gas/liquid interfaces. At one limit, 
say a gas flowing over a liquid, the situation 
approximates to a fluid flowing past a solid 
surface ; at an intermediate stage, say a liquid 
flowing over another liquid nearly identical in 
density and viscosity, the situation is very 
complex indeed ; finally, at the other limit is a 
liquid with a free surface. 


8.1 CONCLUSION 


In this paper the analogy, first developed by 
Lewis and Wuaitman, between heat and mass 
transfer for a system of two fluid phases is carried 
a stage further; the simple stagnant film theory 
is replaced by a boundary layer theory which 
takes into due account the velocity distribution 
arising from the relative motion of the two 
phases. The applicability of the new theory is 
likely to be reduced considerably by the tendency 
to instability of a fluid/fluid interface, but the 
same is true of the stagnant film and Penetration 
theories. 
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Temperature stability of fixed-bed catalytic converters 


H. E. Hoerscuer 
Department of Chemical Engineering, The Johns Hopkins University, Baltimore 18, Maryland 


Abstract—The temperature stability of a fixed bed catalytic reactor is examined for several 
different kinds of operating conditions. The appropriate differential equation containing the 
generation function does not have a simple analytical solution for any but the most trivial 
cases. A generalized method of attack is shown permitting an analysis of the problem for any 
form of the generation function. The method is illustrated in detail for one form of a wall boundary 
condition. Solutions for another case are indicated. 


Résumé—L’auteur exmine, pour différents modes opérationnels la stabilité de la température 
dans un réacteur catalytique a lit fixe. L’équation différentielle adéquate contenant la fonction 
de génération na pas de solution analytique sauf pour certains cas tout 4 fait simples. L’auteur 
préconise une méthode générale d’attaque, permettant l'analyse du probléme pour une certaine 


forme de la fonction de génération. 


La méthode est illustrée en détail pour un cas limité aux parois. 


des solutions pour d'autres cas. 


TEMPERATURE stability in a fixed bed catalytic 
reactor is a problem of obvious importance in the 
design of such pieces of equipment. The tempera- 
ture distribution in a fixed bed catalytic reactor 
has been studied analytically by AmuNDson [1]. 
These papers, although recent, are already classic 
in that they present a complete analytical des- 
cription of the temperature distribution between 
fluid phase and solid surface in a packed bed 
catalytic reactor for the most general case. These 
solutions are necessarily complex and even 
machine computation is difficult. The following 
paper is an attempt to set up a generalized 
technique for estimating the axial temperature- 
length relationship in a fixed bed catalytic reactor 
for any form of the heat generation function. The 
axial temperature is always expected to be the 
extreme (either highest or lowest) value and, as 
such, constitutes an important design criterion. 
The method to be used, while general for any 
form of a generation function (to be defined later), 
is limited to the prediction of the mazimum or 
minimum temperature at any axial position. It 
does not give a detailed temperature-position 
function throughout the bed. As such the resulting 
equations can only be considered to be design 
expedients. They are analytically simple for 


I] indique d’autre part 


certain cases and capable of rapid and simple 
machine solution for most others. 

The problem, then, is to determine that com- 
bination of variables (velocity, bed constants, 
and generation rate) such that a given reactor, 
with a known wall boundary condition, will 
operate at or less than a prescribed fixed maximum 
temperature throughout its length. A similar 
problem would be to determine at what reactor 
length the extreme (maximum or minimum) 
temperature — assumed to be at the axis of the 
bed - reaches a prescribed design value. 

The applicable differential equation for this 
problem is as follows : 

oT oT . 
ag tH SVT + AP) (1) 

This equation makes use of a generalized bed 
constant, K, which is the effective thermal 
diffusivity for the bed. The assumption is implicit 
herein that the gas phase temperature and the 
solid surface temperature are identical. If it is 
desirable to fasten one’s attention on the fluid 
phase only, the constant K should be replaced 
by the eddy diffusivity. Similarly, if the solid 
surface temperature is most important, the 
thermal diffusivity for the solid must be used. 
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In the above equation it has been assumed that 
the radial velocity components are negligible. The 
function, 8(7), is the generation function and is 
a function of temperature. It is desirable to make 
two further assumptions as follows : 


1. Longitudinal conduction is assumed negli- 
gible compared to lateral, or radial, con- 
duction. 


The length to diameter ratio is very much 
greater than one, and hence, the tempera- 
ture profiles may be assumed similar (or 
affine) along the bed length. 


The second statement implies that the tem- 
perature profile at any length may be superim- 
posed upon the temperature profile at any other 
length by changing only a scale factor. It is not 
possible to validate this assumption from existing 
experimental data. However, the assumption 
seems wholly reasonable, in particular in the 
vicinity of the tube axis where (it will be shown) 
the method to be developed is most nearly correct. 

It is convenient to introduce a dimensionless 
temperature and a dimensionless radial position. 
The dimensionless temperature, 7, will be defined 
most conveniently after consideration of the 
boundary conditions of interest to the problem. 
The dimensionless radial position, p, will be 
defined as follows : 

r 


in 


and is zero along the axis of the cylindrical bed 
and unity at the wall. 

Equation (1) then becomes : 
8 tte he 1d 
x CR? \dp* sp Dp 
In eq. (3) the axial conduction term is assumed 
small — as previously stated —- and dropped. The 
generation function §,(r) is related to A(T) 
through a constant term only, depending specific- 
ally only on the form chosen for r. Note, finally, 
that each term in eq. (3) has the dimensions of 
reciprocal time. 

The complete solution of eq. (3), for £,(+) other 
than constant, is difficult and calculation for 
specific problems is frequently equally difficult 


) + B,(r) (3) 


even by machine methods. For the present 


problem, set 
r= 7,, (@, 0) -f(p) (4) 


where +,, is the maximum (axial) + and is a 
function only of z, and @, and f(p) is a suitable 
Note that 


function of p, 


| {(p) dA 
4 
A 


(5) 


where y is a constant. 
Then, introduce eq. 
through by dA, divide through by A, 
integrate term by term. 
Proceeding, eq. (3) becomes : 


(4) into (3), multiply 
and 


[ s) dA 


[ro dA 


Ly 4 





, 


- A) dr 
coefficient of “ and —_™ 


Tm 
are 
0 ax 


Note that the 


r 


. 7 K 
both constant. The coefficient of a is another 


constant, A, which, like y, depends on the form 
of the function f(p). Thus : 


[ B,(r) dA 
—_ kK ToyA } a 
30 y R* A 
The integral term on the right hand side must 
be evaluated separately. It will be a new function 


h(r,,), and will not be a function of p; thus : 


re K Atm ny ) 8 
ala  » R? Si (8) 


dT 


Eq. (8) defines a +,,— @—a surface in space. 
At any point on this surface, the normal will have 


184 





Temperature stability of fixed-bed catalytic converters 


dr. 


direction numbers —", - 1, and, from an 


inspection of the differential equation (8), it 
appears that at this point there will be a tangent 
vector having components 


1, u, (“er +h (r~)| (9) 


Further, all solution surfaces have a normal which 
is perpendicular to this vector. This normal will 
have direction numbers as just stated. Therefore, 
all solution surfaces are tangent to this vector. 
Thus to map the characteristics of the solution 
surfaces, the following system of total differential 
equations must be solved. 


dé dr = 
Se ‘u . 


dr,, 
A ay 
yk 


(10) 


+ h (4%) 


Using eq. (10), it is possible to map out the 
linear elements of which a particular solution 
surface (defined by specific values of the con- 
stants) is made. The +7, — @ characteristics 
defined by eq. (10), when projected on the z — @ 
surface, are straight lines with slope = u. 

The following special cases will be considered 
as examples. They are intended only to illustrate 
the manner in which the resulting equation — eq. 
(10) — may be used. The first case to be examined 
will be for a particular wall boundary condition, 
namely the condition under which no heat is lost 
through the wall to the reactor. The boundary 
conditions are specified and lead to a particular 
form of the temperature profile. This case is 
probably not particularly realistic for a well- 
constructed homogeneous reactor but does serve 
to illustrate the application of the method to a 
practical problem. Case I-B provides the same 
boundary conditions but that the 
generation function is more realistically defined 
in terms of an exponential in the reciprocal of 
temperature. Finally, the industrially realistic 
problem of a constant heat flux at the wall is set 
up and a solution indicated. 


assumes 


Case I — The “adiabatic” reactor 
The dimensionless temperature, 7, will be defined 


as 


T 


* = 7, (x, 0) f (p) (11) 


TV. 

There are four boundary conditions which may 

be written for the temperature — radial position 

(or +—p) relationship, these being the value of 
dT 

1 and the value of — at 


dp 
These four permit a third 


r at p= 0 and p 


p=0 and p= 1. 
degree polynomial to be used to describe the 
temperature profile, viz.: 


S (p) 


The boundary conditions are : 


A 4 Bo Cp? + D,* 


0 (centre), 
sT 


or 


(1) @ p 
p = 0, 


p = 1 (wall), 


Using these boundary conditions one obtains : 


tT = tT. (1 — Bp? + 2p°) (14) 


a 


The form of the profile is shown in Fig. 1. The 
form of the equation for the temperature profiles 
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satisfies all boundary conditions. At this point — 
for the first time -—it is necessary to decide on 
the cross section shape. For a circular cross 
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section, the constants y and A are found as 


follows : 
iy 


1 
7¥= | 2p f(p) dp = 0-8 
0 


1 
2 
1 [ (21+ 22) dp =o 
; dp pop ; 


Thus, the @ — +, relationship defined by eq. (10) 
becomes : 





(16) 


‘ 


where A(r,,), to be designated as the 
integral,” is given by : 


* instability 


B,(r) dA 


h(r,.) = 6 


17 
7 (17) 


A definite form for the generation function is now 
necessary. 


Case 1-A: 
Assume 
A(T) = aT 
Then 
B,(r) = aT, (7 + 1) 


and 





Br) = a(r + 1) 


The instability integral will then be 


| Patra i 
h( 7.) = ar oe = * Jo + 1)dp (19) 
Y Y 


0 
Using eqs. (14) and (15) it may be shown that 


Wry) = 97 + = (20) 


= 06 
Eq. (16) then becomes : 


(22) 


Fig. 2 shows the +, — @ relationship for the 
three possible types of thermodynamic situations. 
Furthermore, from eq. (10) it is also possible to 
compute the axial distance downstream at which 
any prescribed temperature is reached. 


Case I1-B: 


Assume the generation function, (7'), to be more 

realistically defined as follows : 
AT) = ae S/*T (23) 

Then 

a E 1 

Pele) = 5 =P - RT, 7 +1 


The instability integral becomes 


1 
h(r,,) = pi | vex - 
- 0 


E 1 
> a -| dp (24 
RT, |” ti 


Let «= E/R'T, and introduce the 


equation for +: 


profile 


Br p* 
TmP (1 
1+T» 





3 


2 » 
- . 
me (1 — e) will always 
+ To 3 : 
In particular, for the region 


’ 3 
Note that the group 


be less than unity. 
near the tube axis, this group will be very small. 
As a first approximation, it will be assumed 
negligible in order to obtain some information 
concerning the identity of the instability integral. 
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a=0 
Heat of reaction=O 





a<O 
Endothermic™~ reactions 











6 


Thus : 
h(t,) = Fr 
; wy (26) 
3 T,m€ 2 2 7,, € 


[ree |- arene + app] & 
0 


The positive exponential term may be expanded 
in a series, dropping terms in p® and higher with 
the following results : 


1 
_—— me —~ B8Tme 
x -" 1+Tm | fe etre” 


h(r,,.) = Ty 





i, 


\ (28) 





“ 


The last is the first correction term in the indicated 
series integration of eq. (26). For values of « of 
the order of 10, which will be shown to be wholly 
reasonable, the correction term contributes less 
than 20 per cent to the value of h(r7,,) and the 
importance of the correction term diminishes with 
increasing «. Thus : 


a R’ i=tm (1 + 7,,)? 
6Ey To 


h(r,,.) = 


(29) 
and the @ — 7,, relationship is given as before : 


mo e/l+Tm d- 
m 


(30) 


The result of this integration is quite complex. 
The following three equations are quite close 
approximations. The equation for « = 10 was 
checked by graphical techniques and found to 
represent the true @ — +, relationship within 
the permissible limits of the numerical integration. 

7) 2-62 . 


ne of Bl 

7.58 (<2) 
aR 

fore =5 and e. >t 





«= 157, >0-1) 
6 


Fig. 3 indicates the form of the 7, — 6 Ey 
= 5,€= l0 ande = 15. \ aR’ 


Thus, for a value of 


relationship for « 


0 
Note that ———- = 
6 Ey 
—F) 


\a (< eI -70 and e« 


Thus, if « is large, the reactor is likely to be 
unstable, that is, the time required to reach an 
unstable operation will be small. Conversely, if 
a is small, the reactor is likely to be stable for 
large time ranges. If « is negative, it is evident 
from an inspection of equations (31) that 7,, 
will continually decrease with increasing times. 
Finally, if « is zero, 7,, will be constant. These 
possibilities correspond to the thermodynamic 
cases of a large exothermic heat of reaction, the 
unstable case only for long times, the endo- 
thermic reaction, and finally, the case of zero heat 
of reaction. 
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Case II - The constant heat flua case 


Assume that the rate of heat removal (or 
addition) through the walls of the reactor is 
constant at g BTU /hr ft® of surface area. The 
boundary conditions determining the equation 
for the temperature profile parallel those shown 


as eq. (16), viz.: 


(a) @p=0 


Evaluating the constants in eq. (15) one obtains 


io tema} ae] om 


Let 
qR 


TV. 


~ 
*- 


Then 
T = 7, (1 — 8p* + 2p") + 2(p? — p*) 


For this case, 


Sr +2 


Y= 


available 


The @— +7, 
through eq. (10), viz.: 


relationship is again 


= 


This case may now be studied for various forms 
of the generation function as before. 


SUMMARY AND CONCLUSIONS 


A generalized approach to the problems of 
temperature stability of a fixed bed catalytic 
converter has been presented. The method permits 
examination of the conditions under which a 
converter of given length operating on an exo- 
thermic or endothermic reaction of known thermo- 
dynamics and kinetics will maintain an axial 
temperature within some arbitrary limits of a 
design optimum. The method appears to be 
applicable to a wide variety of problems in this 
area and to permit relatively easy exact solution 
for any given problem by analytical methods or 
by numerical or machine computations. The 
most difficult problem presented is the evaluation 
of a definite integral. The solution can often 
be approximated to an acceptable degree of 
accuracy by classical techniques. In addition, 
such problems are readily handled by modern 
computational methods. The method presented 
herein has been used on an industrial problem 
and found to yield meaningful information. 
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NOTATION 


A = cross sectional area 
E = energy of activation 
h (t,) = the instability integral defined by eqs. (7) 
and (8). 
r = radial distance 
R’ = gas constant 
R = tube radius 
T = temperature 
T,, = axial temperature, a function of z only. 
T~ = wall temperature 
wu = linear velocity 





Temperature stability of fixed-bed catalytic converters 


v = a constant defined by eq. (28) = thermal diffusivity 
= a constant 
Greek symbols : = time 

x = a constant, defined by eq. (18) = dimensionless radial position 

B(T) = the generation function defined by eq. (1) = dimensionless temperature 
B,(r), Bolt) = modified forms of the generation function = dimensionless maximum or axial tempera- 
defined by eqs. (3) and (18) ture 
y = a constant, defined by eq. (5) = dimensionless maximum or axial tempera- 


« = defined as E/RT,,. ture atz = 0 


x axial distance 
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Oan \§. ; , 

«| in which the resistance R of 
OaN —g 

the filter plus cake and the resistance Ry of the filter have been based on a simplified Fannine 
equation, has been derived for blocking filtrations of the type in which the filter is blocked 
strongly by the particles of the filtered suspension. Such a filtration will occur in many technical 
filtrations, preceding an actual cake filtration and also in the filtration of very dilute suspensions. 
Experimental proof of the validity of the equation is given for different and differing systems. 


Abstract—A general equation of the form R R,( 


It has again been confirmed that for a certain system the type of filtration is governed by 
the concentration of the suspension. Starting from very dilute suspensions an increase in con- 
centration shows a gradual change from a strong blocking filtration to a cake filtration. An 
indication has been obtained that a third type of filtration, the cake blocking filtration may 
exist. 


Résumé—L’auteur donne a partir de |équation simplifiée de FANNING une équation de la forme 


0 


N 
R= Ry 7) s dans laquelle 2 est la résistance globale filtre + gateau et Ry la résistance du filtre 
a 


N—-g 

seul. Cette équation a été calculée pour des filtrations a pré-filtre (blocking filtration) appartenant 
au type de celles o& le filtre est bloqué fortement par les particules de la suspension filtrée. Une 
telle filtration est fréquente en pratique, et précéde la filtration de suspensions trés diluées. 
La vérification expérimentale de cette équation est donnée pour des systémes variés et suscep- 
tibles de modifications dans le temps. 

Il est & nouveau démontré que pour certains systémes le mode de fitration est imposé par 
la concentration de la suspension. En partant de suspensions trés diluées et eu faisant croitre la 
concentration on peut observer une évolution progressive de la filtration a pré-filtre fortement 
bloqué jusqu’é la filtration avec gateau. On met aussi en évidence un troisiéme type de 
filtration: la fitration a gateau et pré-filtre (cake blocking filtration). 


INTRODUCTION AND THEORY ance ” of the “ cake” (so-called for convenience 
sake) increased with the amount of substance on 
the filter, the more so the nearer an absolute 
blocking was reached. For a cake filtration the 
specific resistance of the cake is independent of 
the amount filtered. A survey of the results 
obtained has been presented graphically in Fig. 12 
of the paper cited. 


IN previous work [1] the phenomena occurring 
during the filtration of very dilute suspensions of 
“ puritine "* over nylon cloth have been des- 
cribed. It appeared that when suspensions 
containing about 1 mg of solid substance per litre 
were filtered a blocking of the filtercloth took 
place. The blocking decreased with an increase 
of the concentration, ranging from an almost An equation for an absolute blocking filtration 
absolute blocking to no blocking (cake filtration). has also been given. Expressed in the total 
For the blocking filtrations the “ specific resist- resistance Rt or in the differential specific 


* An active carbon coal. 
+ For the significance of the symbols see the Nomenclature. 
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resistance r, of the “ cake,” this equation could 
be given as : 

a (a ): 

aN \aN — g. 

For the other types of blocking filtrations found, 
no mathematical expressions have been given. 
Partly because the physical phenomena were not 
clear to us, which was to some extent due to the 
fact that the amount of experiments on which 
the observations had been made were not sufficient 
and only based on one system. Therefore the 
work on the blocking filtration has been expanded, 
in particular to more defined systems. Some of 
the results obtained on blocking filtration will be 
given in this paper. 

Before entering into the experimental details 
some considerations will be given, the choice and 
the treatment of which are determined by the 
systems investigated and by the present opinion 
obtained therewith. 

Let us consider first an ideal system consisting 
of a horizontal filter with a large number N of 
circular pores per unit of surface by means of 
which a suspension of concentration c is filtered, 
which contains spherical particles of weight a and 
of a diameter larger than the diameter of the pores 
and smaller than the distance between thé 
centres of two pores. 

After a certain time @ a volume V has been 
filtered off. This amount V will be considered 
to have flown through in S steps, each step filtering 
an equal amount V’. Therefore : 


V'S=V (1) 


If the filter has a surface A, after the first step 
the number of pores completely blocked will be 
i and therefore the blocked surface will be : 
a 


,? 


én’ if it is assumed that each particle reaching 


an open part of the filter blocks one pore com- 
pletely. 

The surface of the filter A’ still open after 
filtering off V’ therefore will be : 


yee (2) 


V'e 
aN 


If the calculated amount of substance deposited 
per unit of filtering surface after filtering off a 
volume V is g: 


(3) 
Combination of the equations (2) and (8) gives : 


A= A(il— 5.) (4) 
\ SaN} 

In the second step again a volume V’ will be 
filtered. Because part of the pores are blocked, 
the phenomenon that the solid particles will not 
always follow exactly the streamlines of the 
fluid must be taken into account. Let us assume 
as a limiting case that the path of the particles 
over the whole length of the container over the 
filter does deviate less from the vertical than the 
distance between two pores, whereas the liquid 
can move freely. In that case the concentration 
c’ of the second portion to be filtered will be less 
than ¢ and is equal to : 


(5) 


This phenomenon will be called the dilution effect.* 
The still open surface A”’ after the second portion 
has been filtered will be : 


A"=A' _ Ve = 4’ V c A’ -a(1 4 (6) 


aN aNA  SaN 


Combination of (4) and (6) yields : 


” &g 
A" = Aj{l — —— 
( aan) 
Using the same reasoning in the following steps, 
after S steps the free surface A* can be expressed 


by: 


From the filtration data, assuming laminar flow, 
the resistance of the filter R, or of the filter with 


* The phenomenon will also be encountered in the first 
step. Its magnitude will depend on the number of pores of 
the filter and the occurring deviation of the particles from 
the vertical. By the conditions chosen it has for con- 
venience sake been eliminated in the first step. 
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cake R are defined by the special form of the 
FANNING equation : 
1dV 1 AP 


Ad@- » Rior R,) (8) 


For turbulent flow the form : 


1dV 1 (AP)°*s 
Ad@ 7°" Rior R,) 
has been used (see later). 

(In both these equations R and R, are inde- 
pendent of the “filtration pressure’’ AP and the 
viscosity ». 

Considering that the open filter by definition 
has a resistance R, : 

4 4a 
RR, 


and therefore : 


(9) 


(10) 


S. N 
R ;( =} : 
SaN — g, 
The differential specific resistance rz per gram of 
solid substance defined by : 


(11) 


Ve 
in which R, = R — R,, is equal to : 
dk, dk 


dg ~ dg 


, R SaN S+1 
— Hs —<) 

. R 
For all values of S: (rg)p59 = - f. 


The quantity of substance which actually 
blocks the pores g, is smaller than the calculated 
quantity g. As: 


(18) 
it follows that : 


2o = aN (2 = (1 —! (14) 


A second case will now be considered in which 
more than one particle will finally block one 


pore. 


Let us assume that all the pores participate at 
the same time to the same extent in the blocking 
operation. If laminar flow is present in the pores 
and if reduction of the surface of each pore is such 
that the diameter (or the equivalent) decreases 
from D, to D after a volume V has been filtered 
off, then : 

21281 an 
7 D,*N 
128L 


; (16) 
n D*N 


Now if it is assumed that each particle covers a 
; 1 
surface equivalent to 3 7D regardless of the 


actual diameter of the pore, so that : 


from which : 


bp — D? 


combination of (15), (16) and (17) and introducing 


, (7) 


gives : 
naN ) (18) 


R= R, ( 

4\naN —g 

If the flow in the pores is turbulent the same 
equation is valid. In_ that however 


1-38 
D,; 


case 


The partially blocked filter with uniformly 
blocked pores can be visualized to consist of two 
parts, a completely blocked part and an open 
part of resistance R, and of surface A). 

From the considerations given, it is apparent 
that A, can be expressed mathematically as : 


(19) 


The amount of substance filtered has covered a 


Vel 
rf: Sus D2. 
pore surface o — 5 om 


real blocked surface of 


This is equivalent to a 
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l “—— 
3 7D, N 


The real open surface of flow A, therefore is 
equal to : 
'. 2 
404-25 
aN D,? 


Its resistance is : 


n( 


If again the dilution principle is applied to the 
case under investigation and using the con- 
siderations given above we find that after the 
first volume V’ has been filtered off the corres- 
ponding resistance R’ is defined by : 


N 3 
R’ (, me r:| (21a) 
— = _- on 
R, naN — — 


A second portion V’ will give a filter with resist- 
ance R” for which : 


R = (21b) 


; , aw : , 
In this equation n’ = (5) in which D’ is the 
s 
diameter of the pores after filtering V’. Moreover : 


Because : 


J V'e 
(D,)? — (D’?? = D? — 


* AaN’ 
we 2 
“he: (x sts) 


Combination of the equations (21) gives : 


R" SnaN \* 
Sa (sav) (33) 


also : 


Repeating this procedure S times yields for 


R= R': 
SnaN \*% 
= § 1.00 23 
san (<2. -_ -) (28) 


The same equation holds for turbulent flow, in 
which again 


The cases described will not often be encoun- 
tered. The form of the pores and of the particles 
will in general prevent complete blocking. In the 
next case it will be assumed that one pore will 
only be blocked partially by one particle and that 
the cake formed on the partially blocked pores 
will have negligible resistance compared with 
that of the blocked filter. The total volume of V 
will be filtered off in h portions of V’. After the 
first portion there will be an open filter of surface 
A’ (resistance R,) and a partially blocked filter 
of surface 

Mi nin sot oo Mls 
aN haN 
The resistance of this part R,’ can be found by 
applying the equations (15) and (16). One 
obtains that : 
’ D,\* 
R,' = R, (=?) 


and because 


2 
, n 


? n—1 


As far as flow is concerned the part of the filter 
of surface A, and resistance R, is equivalent 
with a filter of surface A,', and resistance R, in 
which 


, 1 ; 
A, —_-— As, 


Because : 
- ; 
i A + A, 


A' +A, 


it follows that : 


Upon filtering a second portion V’, through the 
open filter (with surface A’) will flow an amount of: 


”. By this procedure a second portion 
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A’ Ve. 
of the filter of surface A,’ = - will 


A’ +A'’ aN 
be blocked. Again this is equivalent to a filter 


with resistance R, of surface A,” = A,’ J 

Therefore after the second portion has been 
filtered off, the surface of the free filter is A”, 
the blocked filter has a surface A,’ + A,'’, which 
is equivalent with a filter of resistance Ry, of 
surface: A,’ + A,”. Therefore : 


” 


Aj" +... Af). (24) 


A solution of equation (24) can be approached 
in the following manner. We have seen that 


R, t — 
- haN  phaN 


substituting g/haN =m. Also: 


m 
1—m-+ — 
Pp 
Continuing, it can be derived that : 

m 
= ls _— . >a. . = 


1—m-4 


_ mip 
m 
1—m+— 
Pp 


m 
1—m+ ——m + 
P 


Therefore after the filtration of h portions : 


m 


th=ty—m+Pandt,= 1. (25) 
Ry 


thy 
A direct solution of this series has not yet been 
found by the present author. 
An approximate solution takes the form : 


ty 
As: 
i—¢ 


and because : 


_ m), etc. 


oth 


it follows that : 
4 =1- a —p) {(1 - 


m\* 


P 
After substituting : 


j= a and m ifs 


therefore 


R g \* 
a ' ee 
ee ee (2 stn) ' 


If p is large compared with unity (for absolute 
blocking p = ©) as in the cases investigated, the 
number A decreases to unity. For large values of 


(27) 


p an approximative solution of (27) is: 


q 4 
ee : . 
R aN ~ paN 


This equation can also be derived directly from 
(25). For large p : 


th = ty 


R aN 


> ae 2S § =o ates 
paN raN 
Along the same lines as given before it can be 
derived that if the dilution effect is applied to 


this case, that : 
= ( a= sin) 
SraN 


The equations (11), (18), (23), (28) and (29) are 
valid for filtrations where the filter is blocked 
totally or partly but to a great extent. Moreover 
it has been assumed that if a cake is formed over 
the blocked pores its resistance compared with 
the resistance of the filter is negligible. For every 


(29) 
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cake filtration it is possible and in most technical 
cases very likely that such a type of blocking 
filtration preceeds the cake filtration. 

Taking into account that even in the case of 
complete blocking of the pores, the dilution effect 
needs not always to take place, due to the fact 
that the particles will always follow the stream- 
lines of flow to a certain extent ; that unavoidable 
vibrations of the apparatus influence the move- 
ment of the particles and that the slurry will not 
always be homogeneous, the five cases mentioned 
can be represented by one general equation of 
the form : 


and the equation for rg : 
7 1 
nH 8. (Ny 
aN O \OaN —g 

In these equations O and S can range from unity 

to any value greater than 1. In general, even 

with one particle blocking one pore O + S, the 

limiting case being that O = S = 1, absolute 

blocking of a pore by one particle without 

dilution effect. Another mathematical limiting 

ease will be when S and therefore also O will 
reach infinity. 

For S = O = o the equations (30) and (31) 
transform to : 


It will appear later that before very large values 
of S and O are reached cake filtration occurs. 
For one system, increasing S and O assuming 


the ratio > to be constant, will result in smaller 


values for R and rz at the same values of g, 
showing therefore an apparent decrease in block- 
ing, although in effect the blocking of the pores is 
the same. 

In order to check the two equations (30) and 
(31) it is necessary to know the values of OaN 


or O if aN is known and S. These values can be 
found from the R-V, R-g or r,-g, data in princi- 
pally the same manner. From the R-g (r,-g) 
values, e.g. curves of log R versus log g are con- 
structed. On this curve three points (1, 2, 3) are 
chosen such that log R, — log R, = log R, — 
log Rs. 
For these points : 


log (OaN — g,) — log (OaN — g,) = 
log (OaN — g,) — log (OaN — gs) 
Therefore 
OaN = bs — 818s 
2g2 — (g: + Bs) 





From the values g,, g, and g, therefore OaN can 
be calculated. Via this a plot of log R (or log rq) 
versus log (OaN — g) can be constructed, which 
ought to give a straight line with slope 
— S(or —(S + 1)). 

For rather large values of S and O the log r,- 
logg line can appear to be almost straight, 
especially towards the higher values of logg, 
having nevertheless a slope + zero, therefore not 
representing a cake filtration. In such a case it is 
impossible by lack of accuracy to determine the 
value of OaN. It is therefore impossible to decide 
whether these filtrations belong to the filter- 
blocking type, or that they represent a third 
type, that of a cake-blocking type governed by a 
law of the form : rz = const. (g)’. 

It seems worth while to try to indicate the 
factors governing the magnitude of the para- 
meters S and O. Excluding the vibrations of the 
apparatus for which it is nigh to impossible to 
predict anything in this respect, there must 
amongst other things be an influence of the 
concentration of the slurry, of density differences 
between particle and liquid, of the size-spectrum 
of the particles and of the surface tension between 
particle and fluid. For the parameter O there is 
in general also an influence of the form of the 
pores and of the particles. As far as the con- 
centration is concerned it may be realized that 
each particle will be surrounded by a volume of 


liquid of“. Therefore the distance between the 
c 


L»\ 1/8 
particles will be directly proportional to (?) . 
c 
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The slurry can be thought to be built up of layers 


_» A ; , a\'/8 
containing ( ) particles on a distance ot (¢) . 
a c 


The path each particle has to travel from the layer 


; a\?/8 
nearest to the filter is : - The larger the 
e 


concentration, the smaller will be the length of 
this path. The smaller this length will be the 
less chance a particle will have to come into one 
of the streamlines the 
pores and the bigger will be S and O. In connection 


directed towards open 
with the foregoing the ratio between the amount 
of particles in one layer and the number of pores 
in the filter N will also be important. In the 
region smaller than unity, the larger this ratio the 
greater will be the chance that a particle will 
enter a pore, therefore the smaller will be S and 
O. An 


shorter path, the particle has to travel therefore 


increase in concentration results in a 
in a larger S and O. 

It however also signifies a bigger number of 
particles per layer per pore which will tend to 
Which of the two factors will 
prevail cannot be predicted, although because the 


decrease S and O. 


number of particles per pore is very small in 
blocking filtration, the path influence seems to be 
the largest. 

Increase of the density difference between pore 
and particle will tend to increase S and O. 

A small number of pores will give rather curved 
streamlines near the pores which causes through 
the inertia that the path of the particles will show 
a greater tendency to deviate from the stream- 
lines, increasing S and O. This will also depend on 
the size of the particles, on the density difference 
and on the filtration velocity. The whole situation 
seems rather complex. 


EXPERIMENTAL 
(With the collaboration of: J. P. Branrs (Br.), H. S. 
Joncerrer (J.), and H. J. M. van Sogst (v.S.). 


Besides the experiments already published* several other 
types of filtration apparatus (with horizontal and vertical 
filters) and other systems have been investigated, to check 
the above considerations. Not every type of apparatus 
used will be described here. Only one set-up, as used by 
H. J. M. vAN Soest (Fig. 1) will be discussed. 


* Experiments by H. v. p. Haas (v.d.H.). 


) Heating 
lelement 

















Detail of C 


Fic. 1. Sketch of the filtration apparatus. 


The horizontal filter (C) with the aid of two parted 
metal flanges (see detail) was mounted between two glass 
containers A and B provided with two ground glass 
flanges. Tubes Lem above and under the filter were 
connected with either a vertical or an inclining manometer. 
Two overflows of which one was adjustable enabled to 
establish the pressure difference over the filter. The 
resistance of the whole set-up could be changed by means 
of two stopcocks. The suspension to be filtered was pre- 
pared and stored in a container of 301. (D) made of 
stainless steel 316. One of the difficult problems was to 
prepare a homogeneous suspension. Therefore a smaller 
container of glass (E, 101.) was installed which could be 
stirred better with our equipment and which was filled 
by syphonic action from the large vessel. Overtlowing 
suspension was assembled in the vessels (F) and (G) and 
returned periodically to the container (D) by means of 
compressed air. The filtrate was assembled in the reser- 
voir (H) provided with a calibrated burette, which enabled 
the volume to be read. In some cases a rotameter (/) was 
used. By means of a small displacement pump the filtrate 
was pumped back via a filter (K) to the supply reservoir. 
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The filter proper was placed on a foam rubber cushion 
and the connections with the other parts of the apparatus 
were made of flexible tubing to reduce the influence of 
vibrations as far as possible. All the measurements were 
carried out at a constant temperature, maintained in the 
normal way in the supply-vessel. 

Prefiltered and partly degassed water was always used 
in order to prevent a filter effect caused by impurities 
and (or) air bubbles. 

Special precautions had to be taken with the horizontal 
filters to eliminate the air trapped under the filter. In 
most cases with the filtration of suspensions in water a 
small amount of wetting agent was added to promote the 
formation of a homogeneous suspension by breaking up 
particle-clusters. 

Two liquids have been used: water (p = 0-998, 
1 0-915 cp at 25°C), and a mixture of benzene and 
carbontetrachloride (p = 1:19 and » = 0-803 cp at 25°C). 

Three different sorts of solid substances have been used : 

(a) Two pearl polymers : methylmethacrylate (p = 1-19 
at 25°C); polystyrene (p = 1-045 at 25°C), 

These polymers were graded rather closely by careful 
sieving in the presence of a moist strip of filter paper 
followed by rolling the particles over an inclined glass- 
plate to obtain particles as spherical as possible. The 
mean diameter of the particles was determined by weighing 
about 3,000 particles and by calculating the diameter 
from the mean weight of one particle, assuming that it 
was a sphere. 

(b) “ Puritine ” (as used before) as such or fractions 
obtained by sieving. The fractions are indicated by the 
measures of the two corresponding sieves. 

Five different filters have been used. Two of these were 
nickel filters, made by electrolitical deposition, with a 
very regular pattern of circular pores. The three others 
were nylon filters with more or less rectangular pores. 
Some of the characteristics are given below. 





| Pore 
Indicated | diameter | Thickness 
in mv | in “ | 
= 
| appr. 100 
110 


< 65 - 


Nickelfilter I 
Nickelfilter II 
Nylonfilter I : | 80 
Nylonfilter II | 50 x 40 - 
Nylonfilter III i 15 x 25 - 


260 
280 





The pores of the nickel filter were not quite cylindrical, 
the opening at one side (the prefilt-side) being rounded, 
the other opening being rather sharp. 

The surface of the filters used ranged from 0-8 — 15 em*, 
‘ , “ AV dV r 
lhe magnitudes measured were : } o( vy te rT ) aa 
(AP),, from which R, Ry and g could be calculated, 
because A, c and » were known. 


RESULTS AND DISCUSSION 


For the nylon filters with their rather small pores it was 
not difficult to measure in the laminar region, without 
impairing the accuracy of AP too much. For the nickel 
filters this was not possible and the measurements had to 
be carried out in the turbulent region. To eliminate 
uncertainty for the very thin filters used, the transition 
region between laminar and turbulent flow has been 
determined. At the same time the powers of AP and 7 as 
used already in equation (9) for turbulent flow have been 
established. This has been carried out for the two nickel 
filters (1, 2) and one nylon filter (3) with both the liquids. 
The resistance of the filter Ry calculated from equation (8) 
and the Reynolds number Re calculated from : 
(Re) = weD, 

;, ' 4 dV 
in which : ne, 

AN nD,* d@ 
have been determined at different AP and plotted in 
Fig. 2 as log Ry and log Re versus log AP. From the 
FANNING equation AP = 4f }pw* L/D) it follows that 
for circular pores of diameter D, : 


For the laminar region : 


, j 16 
For the turbulent region (assuming f= (x a)) : 


1 

2-B 1-B 

32,))-8 Lb \— 2-B 
36) AP 


4 
aN (p= 7 28 


1 
2-B 


»D,8 
Re pe. 
B2L 9? 


For rectangular openings analogous expressions have been 


(35) 


used. 

From Fig. 2, using the slopes of the lines as given there, 
based on equations (32), (33), (34) and (35), it can be 
calculated either based on the log Re lines or on the log Ry 
that B = 0-20 (as is often found) and that the 
transition mentioned for the filters used 
between Re = 3 to 7. 

All the measurements with the nylon filters have there- 


lines, 


region lies 


fore been carried out in the region for which the Reynolds 
number was < 3, the resistances have been calculated 
with aid of the equation (8). 

For the nickelfilters the measurements were carried out 
in the region for which the Reynolds number was > 7. 
The resistances have been calculated with the aid of 
equation (9). 
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Fic. 2. Log Ry and log Re as a function of the filtration pressure AP. 


The actual experiments can be compiled most con- 
veniently in two series. The first series will contain the 
experiments carried out with the best defined systems, 
viz. the filtration of pearl polymers over nickel filters. 
The second series will contain the filtration of puritine 
over nylon filters. 

In both series it has been found as before that an increase 
of the concentration caused the transition of a blocking 
filtration to a cake filtration. In the first series however 
this transition was far less reproducible than in the 
second series and more abrupt. In general the repro- 
ducibility in the first series, with the much larger particles 
was rather poor. It often happened that using one system 
with one concentration a blocking filtration was found in 
one experiment a pure cake filtration in another experi- 
ment. Two reasons mentioned before which are difficult 
to eliminate stand out: vibrations and lack of homo- 
geneity. It must be kept in mind that the realization of a 
blocking filtration with rather large particles is a very 
delicate operation necessitating the use of a very dilute 
suspension. The lowest value used in this respect was a 
concentration of 1-8 particles per litre. 

In the series to be presented here the cake-filtrations 
have not been incorporated. Moreover the few experiments 


in series two, where the logr, — logg lines appeared to 
be straight (see before) have been excluded. The experi- 
ments are presented in table form and characterized by 
the parameters S and O or OaN, which includes the first 
and almost most important conclusion that the initial 
stages of all the filtrations as carried out, giving a blocking 
filtration, can be represented in this way and therefore 
can be correlated by means of equation (30). 

Deviations have been found but these are not incor- 
porated. In most cases these deviations could be attributed 
vw some disturbing influence, such as the filter effect, 
actual changes in concentration in the supply vessel, 
sudden observable extra vibrations in the building, forma- 
tion of broken particles by the propellor or by the stirrer, 
small air bubbles in the water and leakage of the filter 
flanges. 

Deviations also occurred near the moment the whole 
filter was blocked with particles. 

In order to give an impression about the determinations 
a few worked out examples, taken from the two series will 
be given. 

In Fig. 3 log rg — logg curves for puritine-nylon are 
presented and the corresponding log rz — log (OaN-g) 
lines in Fig. 4. 
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Fic. 3. The log rg-log g curves for puritine-water filtered over nylon at different concentrations. 


l:c=1x10% 2:c=2x10% 8:c¢=8 x 10° 


In the Figs. 5 and 6 the log R — logg lines and their 
corresponding log R — log (OaN — g) lines for polystyrene 
in water are given. 

The data of the first series are compiled in Table 1, 
those of the second series in Table 2. 

The dilution effect as introduced is clearly marked. In 
some cases it has been possible to give an additional 
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Fic. 4. The log "a log (OaN-g) curves for the system 


puritine-water derived from the lines in Fig. 3. 


l:ec=m=1x10% 2:¢=2x 10° 8:c=8 x 10°% 


4:c=6x10% 5:c=10x10°% 6:c=10 x 10°% 


confirmation in the following manner. After the filtration 
the container above the filter was drained quickly and 
the filter taken out carefully. 

The number of pores blocked by one particle each was 
counted. For experiment v.S.7 (Table 1) this appeared 
to be 4,400. From the volume filtered it was calculated 
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Fic. 5. Two log R-log g curves for polystyrene-water over 
a nickel filter. 
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Table 1. Filtration of spherical particles through a filter with circular pores. 


Systen: : 


Methylmethacrylate—benzene and carbontetrachloride (Ap = 0). 


Filter No. 1, N = 400, D, = 260 », placed 


vertical. 





Br 


Experiment 6 





cin 107 g/cem® 
ain log 

D particles in » 
Particles /litre 








System : Polystyrene-water. Filter No 


400, Dy 260 ww, placed vertical. 





Experiment 


cin 107 g/em® 
ain log 

D particles in » 
Particles /litre 


Br 





System : » water 


Filter No, 2, N 


467, D, 289 «u. placed horizontal. 





Experiment 


cin lug cn 
ain 105g 

D particles in » 
Particles /litre 


v.S. v.S, 


13 





that for this case “s 


was 4,680. If the mean real con- 


centration as far as blocking is concerned is c, for this 
experiment © = 0-94. 


From equation (14) (using the different parameters O 
and S) it can be seen that : 


caN 


“ 


From this 0-93, which is in good agreement. 


By 


counting it was found that from the 280 particles retained 
in the apparatus 100 were already deposited on the 
filter. A photegraph (Fig. 7) of this filter giving the 
underside and the top shows the phenomenon. 
Inspection of the results given, especially of the para- 
meters S and O (or OaN) shows that the conclusions to 
be drawn must be of a quantitative nature. In general it 
can be said that S and O increase with the concentration. 
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Table 2, Filtration of puritine in water over nylon-cloth with horizontal filter 





v. di. v. S 


Experiment v.d.H v.d.il 
= ~ — “ . - 2 1] 
cin 10% g/em® | 1 2 3 ‘ 5 
Fraction in total | total total | 6 j 60-90 60-90 
N 1,000 | 4, | 2,250 2,250 
OuN in Wr g as (| | 80 
s 18 6-0 
OaN/S in 10° ¢ +4 13 





Experiment 


cin 10°% g/em? 





Fraction in 4 . 90-150 150-220 
we 2,2: 2,2. 2.250 2.250 
OaN in 10°% g 2: : 100 
Ss : 3-7 
OaN/S in 105 g 27 

VOL, 

6 
L9 56/ 57 Experiment 


2 35 
4(60- 90) 4(60-90) 4(60--90) 4(60-90) 4(60-90) 
4(90- 150) \(90-150) \(90-150) 4(150-220) \(150-220) 
N 22,50 2,250 2,250 2,250 2,250 2,250 
OaN in 10°5 g 45 50 50 30 0) 60 
s 41 +1 2-2 1-0 1 1-3 
OaN/Sinio*g | 11 12 23 30 


cin 10°8 g/em3 | 
Fraction in yu 





§(60-90) 
§(150--220) 





Experiment » 8 v.S | v.S 
: ; | 
cin 10°% g/em3 1 10 &.f.2 : 10 ; 1 5 


Fraction in » | 050 | 0-50 0-50 0-50 f c 50 | Of 0-50 0-60 0-60 0-60 


N 1750 |1,750 1,750 1,750 ,7! .75 1, 1,750 4.000 = 4,000 | 4,000 
OaN in 10° g 160 /1,110 | 100 1,100 | 8 12 | 20 
s 15) 6-0 S. ‘ | 14 2-0) 2-3) 1-4 
OaN/Sinio*g) 107 | 185 50 5 | 78 4 52} 14 





For a polydispers system with irregular formed particles for which O = S = 1, can be realised up to 90 per cent 


this effect was found to be more pronounced than for a of the total blocking of the pores (Exp. v.S. 29). 


monodispers system. Apparently in the last case by-effects 
have a great influence on the results. This also holds for 
the transiticn from blocking filtration to cake filtration. 
This transition, taking place upon increasing the con- 
centration, in the systems investigated, always occurred 
before large values of S and O had been reached. Appar- 
ently at a certain concentration the distance of the particles 
is such that bridge formation over the pores overlaps the 
dilution phenomenon. Under optimum and carefully 
chosen conditions and with some luck with regard to the 
non-controllable influences an absolute blocking filtration 


There is an indication that for monodispers systems this 
optimum lies at a certain concentration, depending on N 
and particle size and Ap, below which S decreases, above 
which S increases with increase of the concentration. The 
possibility of such a behaviour has been indicated in the 
discussion and can be found by plotting S and versus 


1 \ c 
—~ 10) . 
a8 a 
In some cases with the puritine-water system for low 


concentrations S has been found to be <1. This is an 
impossibility. The reason must be found in some 
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Fic. 6. The log R-log (OaN-g) lines derived from Fig. 5. 


uncertainty with regard to the zero point of the volume of 
the filtrate. A change in the zero point for V will change 
the S found. Although in this respect a correction is 
possible by displacement of the R-V' lines such that 
Ryo = Ry, for low concentrations with a relatively 
low value of 7 at low values of }’ some uncertainty will 


remain. With the fractions containing small particles (the 
original puritine consisted for 85% of the 0-50 » fraction) 


OaN , = 
~ has been found to be constant. This signifies that 


R 
In these cases Tdy_9 = =: 


concentration. For the 0-50, puritine fraction with 
filter No. 3 (rather large openings) it could be observed 
that some substance passed the filter. This effect has the 
same influence as the dilution effect and will increase the 
magnitude of S and O. 

A marked influence of the position of the filter, of N 
and of Ap has not been found. Filtrations under still more 
idealized conditions, e.g. in a building without vibrations, 
must be carried out to find these influences. It is question- 
able whether such experiments are of importance for 
technical filtrations. 


O-=S. independent of the 


CONCLUSION 
A blocking filtration in which the filter is blocked 
by the particles and in which a strong blocking 
occurs can be correlated by a general equation 
r \s 

of the form R = R, (can .) both for laminar 
JaN —g 

and turbulent flow. 
The parameters S and O in general increase 
with an increase of the concentration of the 
prefilt. After a certain concentration has been 


reached the filter-blocking filtration changes to a 
cake filtration. This change is gradual. There is 
an indication that a third type of filtration may 
exist, viz. the blocking of the cake. If this should 
be reality, this filtration stands in between the 
filter blocking filtration and the cake filtration. 

A blocking filtration is defined as a filtration 


in which = HS) whilst in a cake 


filtration rz + f(g). The occurring of a filter 
blocking filtration and a cake filtration coincides 
with the deposition of the cake in and over the pores 
of the filter. If the pores in a blocking filtration 
are not blocked strongly an analogous equation 
between R and R, of a somewhat different type 
can be derived. This equation has not yet been 
checked experimentally nor the phenomenon been 
realized. It seems very likely that a moderate 
blocking will practically never occur. 

Either the particles enter the pores of a filter 
and than almost always will give rise to a great 
increase of the resistance (a decrease in diameter 
of 4 already signifies a 16 times larger resistance) 
or they form a cake over the pore. This seems 
the more true because particles, very small with 
respect to the diameter of the pores will pass the 
filter. 


NOMENCLATURE 


= weight of a particle Mg 
= surface of the filter em? 
= power of Reynolds number in friction factor 
= concentration of prefilt gem 
= diameter of pore em and « 
= equivalent diameter of surface covered by 
one particle 
friction factor 
amount of solid substance per unit of filter 
area 
parameter 
= length of a pore 
eS 
~ haN 
= factor 
= number of pores per unit of filter area 
= parameter 
= factor 
= pressure difference over the filter 
gem 


em and yu 


gem 


em and yu 


' sec”? (dyne em™) 
= factor 
= differential specific resistance of one g of 


solid substance em g™ or (em g™!)? 
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Fic. 7. The nickel filter with polystyrene particles, seen 
from the underside (left) and from above (right). 
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total resistance of filter and solid substance = term of series 
em” or g° em-? ’ = volume of filtrate 
= resistance of cake em”! or g°° em-?** velocity of fluid in a pore 
= resistance of filter em™! or g?™ em"! density 
= Reynolds number viscosity 
= parameter = time of filtration 


If no indication is given the magnitudes are dimensionless. 
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(Received 5 February 1956) 


Abstract—Preliminary tests have demonstrated the valuab!e versatility of mercury evaporation 
as a technique for the study of mass transfer systems. The use of a sensitive ultra-violet absorp- 
tiometer for measuring the concentration of mercury vapour in gas streams enables studies 
to be made of systems with hydrodynamic control, the results agreeing with published correlations. 
The sensitivity of the detector aJlows measurements of transfer from very small surfaces, leading 
to surveys of positional effects. The choice of copper alloys as backgrounds for amalgamation 
gives great versatility in the study of transfer from any selected shape and also provides the 
possibility of introducing selected interfacial reaction kinetics into the process. 


Résumé—Des essais préliminaires ont montré la souplesse qu’offre évaporation du mercure 
comme technique pour l'étude des phénoménes de transfert de masse. La mesure de la con- 
centration du mercure dans la phase gazeuse par absorption dans l'ultra-violet permet l'étude 
des phénoménes contrélés par l'écoulement gazeux : les résultats concordent avec les corrélations 
publiées. La sensibilité du détecteur permet des mesures de transfert a partir de surfaces extré- 
mement petites et donne une idée sur les effets de position. Le choix d’amalgames a base d’alliages 
de cuivre offre plus de possibilités dans l'étude des transferts a partir d'une forme géométrique 
choisie et per.net ainsi d’introduire la cinétique des réactions interfaciales dans le phénoméne 


d'échange. 


Tue study of evaporation systems for the correla- 
tion of mass transfer coeflicients is restricted when 
using liquids such as water by the difficulties of 
supplying liquid evenly to any specific shape of 
surface and supplying energy to maintain a 
system controlled by the resistances to the 
evaporation mechanism, and not by the energy 
transfer system from the fluid above or the 
material beneath the liquid. These restrictions 
can be avoided by the use of a liquid of low 
vapour pressure giving low evaporation rates 
requiring low energy transfer rates. Mercury 
surfaces appear to fulfil this requirement and to 
have the great advantage that the geometry of a 
desired system can be obtained readily in copper 
or brass and transformed to a mercury system 
by amalgamation. Once the clean mercury 
surface is formed it can be used for a short time 
without any significant change in the surface as 
evaporation proceeds. The work reported here 
consists of the initial tests to check the possible 


use of mercury transfer system and to assess their 
limitations. In using other low rate mass transfer 
systems, e.g. 2-naphthol into water, napthalene 
into air, the available shapes are limited and 
relatively fragile, but the main problem becomes 
the technique for measuring accurately very low 
concentrations of the transfer component in the 
moving fluid phase. With mercury in minute 
amounts in a gas stream the ultra-violet absorp- 
tiometer provides a powerful analytical tool. 
Instruments designed specially for the detection 
of mercury vapour in gas have hitherto been rare 
but the recent availability of the Hanovia 
apparatus has provided a ready and most sensitive 
method of analysis, which now makes it possible 
to use the mercury system to study many mass 
transfer problems and to exploit the ease of 
making an evaporative surface of any desired 
form. As the instrument is uncommon a brief 
comment on its operation and sensitivity has 
been included as an appendix. 


*Present address: Levington Research Station, Levington, Nr. Ipwich, Suffolk. 
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Apparatus for mercury evaporation transfer from shaped surfaces. 
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Fic, 2. Transfer Tubes and shapes. (a) Rods. (b) Spheres. 














205 





R. W. 


APPARATUS 

To enable tests to be made on transfer from rods and 
spheres a glass apparatus was constructed with the flow 
system shown in Fig. 1. Air or nitrogen was fed through 
a blower F, a Variac-controlled electric heater H, and a 
suitable rotameter R, or R, to one of the glass pipes of 
internal diameter A, lin.; B, 2in.; C, 4im. Mercury 
vapour evaporated from a surface within A, B, or C 
left the pipe with the gas stream. A continuous sample 
from this stream was passed through a rotameter Ry 
and if necessary blended with fresh air from rotameter 
R, before flowing through the absorptiometer M for 
analysis. All exhaust gas containing mercury vapour 
from vents V was discharged outside the laboratory 
building. The major experimental adjustments consisted 
of controlling the input gas stream to pipe A, B, or C 
for both rate and temperature, and controlling the sample 
dilution ration with R, and R, when necessary to maintain 
the vapour concentration in the gas stream to the analyser 
within the desired range for accurate measurement. The 
sample rate through the analyser was always above the 
limit for avoiding oxidation of the mercury vapour in 
the analyser (see appendix). 

Gas entering pipe A, B, or C passed through three 


Maxwe.tt and J. ANDERSON SToRROW 


screens of 30-mesh copper gauze to remove major swirls 
in the flow approacning the evaporative surface. With 
limited space available for the apparatus it was not 
possible to provide adequate calming sections. The major 
dimensions in Fig. 2 are given for the experiments 
discussed, as the disturbance in the fluid meeting the 
evaporative surfaces may have been significant. Brass 
rods were carried axially in the pipes by screwed extensions 
located by the perspex carrier baffles b, and bg, the baffles 
being “ sliding fits’ in the appropriate glass pipe. The 
baffles were all perforated, b, cut away in quadrants to 
give the maximum opening and bg with holes near the axis 
to provide mixing of the outlet gas stream before it reached 
the sample withdrawal tube. By “ painting" with 
mercury, any selected part of the rods or spheres could 
be tested for evaporation rate, allowing length effects 
and circumferential position effects to be studied, so long 
as the evaporative sections were kept a few diameters 
away from baffle b,. Nonevaporative sections of the shapes 
were lacquered to prevent unwanted deposition of mercury 
during the preparation of the desired shape of surface 
and also during experiments from the gas stream on 
to copper surfaces downstream from the evaporative 
surface. 
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Detector scale reading 
iN) 
2) 





Run Material Temp*C Re 
1 copper 18-19 2260 
2 brass 18-19 2260 
3 copper 17-5 2780 

brass 17°5 2780 
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Fic. 8. Effect of time on mercury evaporation rate from copper and brass rods of 0-5 in. diameter in 1 in. diameter tube. 
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EXPERIMENTAL RESULTS 
The results were expressed in terms of mass 
transfer coefficients k based on the evaporative 
area and the logarithmic mean of the partial 
pressure difference (p* — p) 
length of the evaporative section. 
are given in terms. of Sherwood numbers or 


along the axial 
Correlations 


Colburn j factors in conjunction with Reynolds 
and Schmidt numbers based on the physical 
properties of the carrier ges es the smal] mercury 
content of the strecm had a negligible effect 
on these properties. The diffusion coefficient for 
mercury vapour into the gas stream was assumed 
to be the same <s that for mercury into nitrogen 


[5]. The data were corrected for experimental 
conditions by the form 


D= D(z.) 5 


Detector scale reading 


oa 
° 


The mercury surfaces were prepared on copper 
or brass backgrounds by rubbing mercury on to 
surfaces which had been 
allowed to stand in mercuric salt solution. 
mercury 
series of tests. 


and 
The 


short 


cleaned in acid 


surface wes renewed for each 
The effect of time on the transfer 
rate is demonstrated by the variation of the 
The data for 
the copper base show a small fall in transfer 
rate with increasing time, the data being plotted 


outlet gas concentration in Fig. 3. 


as scale readings from the absorptiometer on the 
lower sensitivity range shown in the appendix. 
Only the copper shapes show this low rate of 
decrease over a period of hours. The data for 
brass bases of the same dimensions in Fig. 3 
show a characteristic deviation from the ageing 
effect on copper. Brass bases show a slow initial 
fall, and agreement with the copper-base data, 
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Fic. 4. Initial ageing lines for brass rods } in. diameter in a tube 1 in. diameter with air between 14 and 17°C. 
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followed by an abrupt change to the exponentially 
decreasing form which appears to approach an 
asymptotic value after long evaporation periods. 

From these initial studies, the evaporation 
process was clearly sensitive to flow, showed 
reproducibility between shapes of the same 
dimensions and also differences in ageing 
functions depending on the background metal. 
The last suggested that the mercury system could 
well provide a tool of wider scope than antici- 
pated, with probably reaction kinetics systems 
available at the bounding surface. The sensitivity 
of the measuring apparatus and the consistency 
of measurement of smal] mercury concentrations 
also showed that the transfer system could be 
used on small transfer areas providing an unusual 
versatility for purposes of survey over surfaces. 
The tests now reported were undertaken primarily 
to check that the experimental apparatus could 
produce results comparable with accepted correla- 


tions from other systems and retain sensitivity 
to all major variables. Data from copper and brass 
surfaces were extrapolated to zero time from 
ageing plots, as in Fig. 4, comprising only the 
initial falling rate periods. They show below 
that as a good approximation the gas concentra- 
tion at zero time may be taken as that corres- 
ponding to a nominal mercury surface, with the 
slow ageing effect to be considered as a minor 
modifying function. 


Evaporation from spheres 

With this assumption, that the zero time data 
were for a pure mercury surface, the mass transfer 
from spheres is exemplified by the data in Fig. 5 
for the complete surface of a lin. diameter 
brass sphere in the pipe C. A_ satisfactory 
correlation is obtained using the physical charac- 
teristics of air, with reproducibility for two 
series of tests on separately prepared surfaces, 
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Fic. 5. Transfer from 1 in. diameter sphere in 4 in. diameter tube showing agreement with accepted correlations. 
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and agreeing with Frorss.ino’s equation for the 
evaporation of water, aniline, nitrobenzene, and 
naphthalene from spheres [6] : 


k.d, i , V, d, p . La ’ 
p-2 |r roar (=) (S)) 
This equation is confirmed by the data of other 
authors over the experimental range in the present 
work. 

These data were for a spherical surface complete 
except for the relatively small disturbance caused 
by the holding rod behind the sphere. The 
influence of position on the sphere may be 
assessed by testing a band around the sphere 
at various positions along the tube axis. The data 
on Fig. 5 for a band from @=0 to @ = 42° 
on the nose of the sphere (Fig. 2b) show that the 
apparatus is sensitive to the position effect 
despite the smallness of the transfer surface. 
The reproducibility with the small area is worse 
than with the complete sphere as must be 
expected, but the full sphere data compare very 
favourably with the bands of transfer data 
normally accepted as justifying correlations, 
such as the Froessling equation. The value of 
the technique is also shown in the large range of 


the abscissa, limited at the high rates by the 
blower used. The extreme sensitivity is exem- 
plified by the point at the lowest abscissa included 
as a single test at the very low flow rate of 
0-5 litre /min. 


Evaporation from rods 


The apparatus in Fig. 2a provided the j,, 
data in Fig. 6 showing agreement for different 
diameter ratios d,/d, with 7 = f/2 derived from 
the appropriate momentum transfer equation 
for laminar motion in an annulus [4] 


f= 16 (7-5). 4 (/ay, (2) 


(d/? vai d,*) 
log, d,/d, 


— (43 + a], 


#d,/d,) = (a, — a)? / 
and for turbulent flow in the annulus {1] 
—0-20 0-10 
f = 0-055 (“%-2) (“>“) - (8) 
B d, 


The equation (8) for 7 was found to require a 
modification in the diameter ratio parameter to 
correlate heat transfer data [8] 
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Fic. 6. Transfer to air flowing in annuli between rods and 
tubes of various diameters, showing effect of diameter ratio d,/d,. 
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(4) 


The j lines for turbulent flow on Fig. 6 were 
based on equation (4). The difference between 
j values from equations (3) and (4) is significant 
for the diameter ratios in the present work. 
The data confirm the validity of equation (4) 
for relatively large d,/d, ratios. Fitting smooth- 
ing cones, of length equal to eight diameters 
of the rod, to the upstream face of the long 
cylinder (Fig. 2a) did not appear to affect the 
the 
expected 


transfer coefficients significantly, though 
flat end of the cylinders would be 
to show increased transfer by disturbing boundary 


layer formation, particularly in the laminar flow 


region. This effect is being studied further in 
tests on the distribution of local transfer coefti- 
cients along the length of the cylinders. 

The j,, data agree very well with the narrow 
band of j for a wide range of d,/d, in laminar 
motion and also show the expected differences 
in the turbulent flow region, agreeing well with 
the appropriate equation (4). 


Ageing effects on mercury surfaces 

In Fig. 3 the outlet gas concentration at fixed 
flow rate show two periods of falling transfer 
rate. The first falling rate period appears with 
both copper and brass base surfaces and may 
the 
reduction in 


reflect an oxidation reaction on 


giving 


mercury 


surface progressive the 
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evaporative surface. With a brass foundation 
this period is followed by a more rapid exponential 
function with an apparent asymptotic value for 
transfer rate after long periods which is only a 
small fraction of the initial rate. This period may 
be concerned with the diffusion of zinc from the 
brass base to the mercury-gas interface, there 
affecting the oxidation and transfer processes. 
It was observed that over the periods shown 
for the tests the mercury surfaces remained 
“mirror clean”’’ and the apparatus is capable 
of sensing changes which are invisible. These 


processes are being studied in more detail but a 
preliminary comment is considered valuable 
in that the transfer system has obvious posibilities 
for studying overall transfer mechanisms involving 
moving phase dynamics and surface reactions, 
with modified surface reaction forms available 
by the choice of the amalgamation background. 
Should further studies of mercury condensation 
on metal collecting surfaces also prove successful, 
the mercury transfer system may provide an 
unusual versatility of great value. 

The presence of oxidation reactions affecting 
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Fic. 8. Influence of air flow rate on ageing functions with annulus of d, = } in., d, = 1 in., at 18°C with brass surface. 
Each curve is labelled with its constant flow rate in litres/min. the full line represents similar data at 36°C. 
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the ageing rate is demonstrated by the data on 
Fig. 7. For all gas mixtures at constant flow 
rate and temperature, the first falling rate 
period varies little though showing an ageing 
with oxygen which does not appear with cylinder 
nitrogen. The second period shows a much 
faster ageing which increases with increasing 
oxygen concentration. In these tests the enhanced 
oxygen concentrations were obtained from 
oxygen injection to the blower through a rota- 
meter, the content being controlled by rates. 
For the pure oxygen and pure nitrogen data 
cylinder gases were fed directly to the heater H. 
The variation of apparent break-point with 
oxygen concentration may be merely fortuitous 
as the data have been pinned to an arbitrary 
zero of time for each test. The mercury film 
thickness may be a major variable controlling 
the incidence of the second falling rate period. 

The dependence of the ageing functions on 
flow rate is shown in Fig. 8 which includes data 
for air flowing through an annulus at rates from 
17 to 145 litre/min corresponding to Re, from 
760 to 6,450. The variable equivalent to transfer 
rate is given by the product of outlet concen- 
tration and flow rate, and further analyses of the 
effect are required, but the apparatus discri- 
minates successfully. The curve for data at a 
higher temperature suggests the ability to provide 
analyses of the effect of this variable to assess 
its influence on the hydrodynamic and reaction 
mechanisms. The curves on Fig. 8 also show the 
need for a mechanism model which gives an 
inversion of curvature rather than an abrupt 
break-point separating the two falling rate 
periods. Further work on these problems is 
now proceeding. 


CONCLUSIONS 

The preliminary tests on the transfer of mercury 
vapour into gas streams have demonstrated that 
the system provides a valuable technique. By 
choosing the background metal for amalgamation 
and the gas for the fluid phase, it may provide 
a choice of control by hydrodynamic, surface 
reaction, or diffusion mechanisms. The hydro- 
dynamic mechanism has been shown to agree 
closely with established correlations and the 


system can be used with ease to study the 
effect on local transfer coefficients of position on 
the selected surface shapes. The sensitivity 
of the detector apparatus is such that the survey 
of mass transfer on surfaces by amalgamation of 
a small test area appears more versatile than 
previous techniques and may provide a powerful 
tool for fundamental studies. The selective 
introduction of diffusion and surface reaction 
kinetics may extend the usefulness of the 


technique to reproduce a variety of processes 
occurring on any geometric form of surface, thus 
covering a very large number of problenws. 


APPENDIX 
Calibration of mercury vapour detector 

The study of the proposed transfer systems rests 
entirely on the use of an extremely sensitive detector for 
mercury vapour in a gas stream. The availability of a 
mercury vapour lamp as source of radiation of limited 
spectrum for the absorptiometer is a fortunate chance, 
and although the instrument can be used for measuring 
organic vapours in gas the high sensitivity to mercury 
is of major value. For stability the mercury vapour lamp 
must be supplied from a controlled stabilised voltage 
source and a balancing bridge circuit is required for 
measuring the output from the receiving photocell. 
The Hanovia instrument now incorporates both these 
features and although developed primarily for limit 
indication in toxicity surveys, its reproducibility, 
sensitivity and stability were exceptional for this type of 
absorptiometer. The instrument used in this work had 
two ranges of sensitivity. The calibration on the lower 
sensitivity range is shown in Fig. 9, the data showing the 
reproducibility after considerable use. Although the 
experiments discussed were concerned with evaporation 
from areas down to one square inch the lower sensitivity 
range was used throughout. The higher sensitivity would 
have given full scale deflection for a concentration of 
100 microgram mercury per cubic metre of air, which 
could lead to more detailed surveys with either smaller 
areas or higher velocities than discussed above. 

The use of such low vapour concentrations entails 
difficulty in calibration. The analysis of gas samples by 
chemical absorption methods was considered unsatis- 
factory and the calibration was based on the known 
function of vapour pressure of mercury with temperature 
[2]. It may be noted that the present calibration fell 
between two curves based on absorption of mercury in 
dithizone reagent which differed by a maximum of 30%. 
Air taken from outside the laboratory building passed 
through a saturating tube and was mixed with fresh air 
through calibrated rotameters before flowing through 
the absorptiometer. The calibration of the detector 
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Fic. 9. Calibration curves for mercury vapour detector. 


thus rested on the measured flow rates and presumed 
saturation of one stream. The saturation concentration 
was provided by using a glass column packed with rolled 
fine-mesh brass gauze with an amalgamated surface, an 
extra column being added occasionally to check that 
doubling or quadrupling the surface did not effect the 
concentration of the outlet gas within the sensitivity of 
the detector. A successful saturator for gas rates up to 
2 litre/min comprised a column of 2 ft length and 3 in. 
diameter packed tightly with 30-mesh gauze, equivalent 
to about four times the minimum saturation column. 
The only limit on calibration and experimental flow 
rates through the absorptiometer is a minimum below 
which the mercury vapour is noticeably oxidized while with- 
in the ultra-violet absorption tube, the meter reading then 
depending on flow rate. In the present apparatus the 
limiting rate was 5 litre/min but an acceptable limit 
should be determined for any other detector as it clearly 
depends on residence time and hence on the geometry 
of the absorption tube. 
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NOMENCLATURE 

d = Diameter: sphere d,; tube, internal d,; rod, 
external d,; annulus, effective d,=d,—d, 
(cm). 

D = Diffusion coefficient: Dg at standard pressure 
Py and temperature Ty (cm?/sec). 

k = Evaporation transfer coefficient. [g/see cm? 
(g/c.c.)]. 
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= Partial pressure of mercury vapour in gas; Dimensionless groups 
p* in equilibrium with liquid mercury (atm.). = Friction factor. 


= Total pressure (atm.). Momentum transfer factor ; j = {/2. 


° a . k f 
= Absolute temperature (°K). = Mass transfer factor ; j,, = p° (Sc)*/8, 


ial Vad, p 


" = Mean velocity in annulus (cm/sec). a ne ghee, Me 
: ) : : " 
m 


= Mean approach velocity for sphere (cm/sec). Vd 
,; annulus, Re, = ——*?. 
= Density of gas (g/c.c.). . ’ i 


= Viscosity of gas (g/cm sec). y Schmidt number: Sc = yp, D. 
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Diffusion in a ternary gas system with application to gas separation 


J. J. Keyes, Jr.* and R. L. Picrorp 


University of Delaware, Newark, Delaware, U.S.A. 


Abstract— Measurements are reported of the separation effected by diffusion of a gaseous mixture 
of hydrogen and nitrogen through two organic vapours for the cases of (a) stagnant vapour and 
(b) countercurrent flow of the vapour relative to the hydrogen and nitrogen. Comparison of the 
data with classical diffusion theory is facilitated by the use of simple geometry in the diffusion 
zone, and by careful control of the experimental boundary conditions. 


In general, the data agree well with theory based on a one-dimensional diffusion model, for 
both stagnant and flowing vapour. The role of countercurrent vapour flow in effecting large 
enrichments is particularly emphasized. For example, the measured separation factor for a 
hydrogen-nitrogen mixture was increased about thirty-fold when the mass flow of vapour was 
raised from zero to about three times the rate of feed flow. 


VOL. The theoretical analysis is extended to determination of performance characteristics for a 
6 practical diffusion stage separating isotopes of low natural abundance of the light component. 
1956/57 


Résumé— I] s‘agit d'un compte-rendu de mesures de séparation effectuées par diffusion d'un 
mélange gazeux de Hg, et Ny a travers deux vapeurs organiques dans deux cas : 

(a) vapeur au repos. 

(b) écoulement a contre-courant de la vapeur avec le mélange No — Ho. 

La comparaison des données avec la théorie classique de la diffusion et facilitée par lusage de 
considérations géométriques simples dans la zone de diffusion et par un controle soigneux des 
conditions expérimentales aux limites. 

En général, les données concordent bien avec la théorie basée sur le modéle unidimensionnel, a 
la fois pour la vapeur au repos et la vapeur en mouvement. Les auteurs ont particuliérement 
insisté sur le réle de écoulement a contre-courant pour effectuer d'excellents enrichissements. 
Par exemple, le facteur de séparation mesuré pour un mélange Hy—Ng est accru environ 30 fois 
pour une élévation du flux massique de la vapeur de 0 4 3 fois le taux d’alimentation. 


L’analyse théorique est étendue a la détermination des caractéristiques de fonctionnement 
pour un étage pratique de diffusion comme de la séparation d’isotopes rares de composés légers. 


INTRODUCTION various practical stage configurations, in which 
the overall separation depended upon both the 
diffusive enrichment and concentration gradient 
effects normal to the direction of diffusion. This 
complicated comparison of the data with any 
theoretical treatment of multicomponent diffusion. 

Benepicr and Boas [1] have published an 
excellent analysis of the mass diffusion process 
for two special cases : 


Ir nas been recognized [1, 2, 5, 8, 9] that partial 
separation of the components of a binary gas 
mixture can be effected by diffusion at constant 
total pressure through a third gas or vapour, 
provided there exists a difference in diffusivities 
of the components with respect to the vapour. 
The third component is commonly called the 
separating agent. 

Past experimentation with this process, which (1) Isotopic mixtures for which the diffusivities 


has been referred to in the literature as ‘“* mass with respect to separating agent of the com- 


diffusion,” “ atmolysis,” and “ sweep diffusion,” ponents to be separated are nearly identical. 


emphasized the performance characteristics of (2) Non-isotopic mixtures for which the mutual 


* Present address : Union Carbide Nuclear Company, Oak Ridge, Tennessee, U.S.A. 
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diffusivitity of the components to be separated 
equals that of one of the components with 


respect to the separating agent. 


Other methods for handling multicomponent 
diffusion calculations, such as those of CoLBURN 
and Drew [3] and WILKE [12] are not sufficiently 
general to cover the important case of counter- 
current flow of separating agent. Since the 
treatment of Benepict and Boas is also somewhat 
limited in application, it was felt that, in addition 
to experimental measurements of the diffusive 
enrichment effect, a more general analysis of 
ternary diffusion would be useful. Comparison 
of experiment with theory for well defined diffuser 
geometry should, it was felt, provide a firmer basis 
for further application of mass diffusion. 

Accordingly, experimental techniques and con- 
ditions were designed to reproduce as nearly as 
possible one-dimensional diffusion which occurs, 
for example, within a single “ pore’ of a mass 
diffusion partition, where the primary enrichment 
effect occurs. Two practical cases are treated : 

(i) Diffusion of the binary mixture into a stag- 
nant separating agent. 

(ii) Diffusion of the mixture against the direction 
of flow of the separating agent. This situa- 
tion, which might be termed “ countercurrent 
gas extraction,” results in greatly enhanced 


and R. L. Picrorp 


separation at the expense of throughput, a 
consequence of the effect of the counter-flow 
in selectively retarding diffusion of the heavy 
component, permitting the light one to 
progress against the flow. In previous work, 
this particular role of separating agent flow 
has not been adequately emphasized or 
evaluated. 

The theoretical analysis with which the data 
are compared is derived from the classical theory 
of MaxweELt and Sreran [10], as employed also 
by Benepict and Boas. As will be evident, this 
simple approach is indeed adequate to describe 
one-dimensional diffusional separation. 

The results will also have application to such 
problems as diffusion controlled catalytic reactions 
wherein one of the products must diffuse away 
from a surface in opposition to the inward 
diffusion of the reacting gases. 


APPARATUS 
Several requirements to be fulfilled by the 
apparatus were recognized : 
(1) Steady flow of diffusing components had to 
be provided. 
Convective effects had to be minimized. 
Boundary conditions assumed in the mathe- 
matical derivations had to be accurately 
reproduced in the experimental procedure. 


(2) 
(3) 
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Fig. 1 is the flow diagram. The diffuser con- 
sisted of a four-foot section of 1 in. copper tubing 
loosely packed with glass wool (97°, free volume) 
this design having been found to minimize con- 
vection. The tube was jacketed for steam 
heating, and maintained at 100°C. Feed consisting 
of a binary gas mixture was admitted near one 
end of the diffuser at 100°C, 1 atm., where it was 
divided into two streams, the product diffusing 
toward the separating agent inlet, and the 
raflinate discharging close to the feed point. The 
separating agent, in this work a condensable 
vapour, was circulated at the product or discharge 
end of the diffuser by means of a reboiler-con- 
denser combination indicated in Fig. 1. Because 
the rate of vapour circulation greatly exceeded 
the diffusion rate of the product stream, the 
concentration of the binary mixture at the point 
of product discharge was reduced essentially to 
zero. As will be seen, the assumption of negligible 
concentration of the diffusing mixture at the 
product end of the diffuser effects considerable 
simplification in the mathematical analysis, thus 
facilitating comparison with the data. 

Product and raffinate streams were quantita- 
tively freed of vapour by means of the condensers 
and cold traps indicated, and the residual binary 
mixtures were analysed by measuring their 
thermal conductivity relative to that of the feed 
mixture. Since the thermal conductivity cells 
gave vapour-free analyses, a method for deter- 
mination of vapour concentration was developed. 
This involved sampling directly from the diffusion 
tube, condensing at constant volume, and applying 
the perfect gas law to compute the concentration 
of the third component from the changes in 
pressure and temperature. 

Flows of all streams were controlled by means 
of micrometric needle valves operating as critical 
orifices. The extremely small gas rates involved 
in this work (about 100 std. cm® per hour) were 
measured by determining frequencies of bubble 
formation in a non-volatile oil. 

All runs were made at essentially atmospheric 
pressure. Oil-pumped tank nitrogen and pure tank 
hydrogen gas were used as feed materials. The 
separating agents used were analytical reagent- 
grade organic vapours. 


EXPERIMENTAL PROCEDURE 
AND RESULTS 


Diffusion through stagnant separating agent (Table 1, 
Case 1) 


No raflinate stream was removed: all of the feed gas 
diffused through the tube. This is equivalent to operation 
at total reflux, since no net material leaves the system. In 
addition, there is no net flow of separating agent. Of 
course continuous production of enriched mixtures is not 
possible, but there is a variation in composition of the 
feed mixture along the diffuser in order to maintain the 
proper diffusion rates of the two components. Fig. 2 
illustrates this case. It is desired to compare measured 
and calculated compositions, utilizing the theory outlined 
in the Appendix and summarized in Table 1. 
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Fic. 2. Illustration of operation with diffusion of gases A 


and B through stagnant vapour C. 
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Fic. 3. Illustration of operation with countercurrent 
diffusion of gases A and B through vapour C. 
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where yp is the vapour-free mole fraction of light com- 
ponent in the product stream (this is also equal to the feed 
composition, z, as is apparent from material balance con- 
siderations), and y, is its mole fraction at the inlet, just 
inside the feed point, on a vapour-free basis. Fig. 4 com- 
pares a few of the measured and calculated separation 
factors. The feed mixture was 45% hydrogen, 55% 
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Fic. 4. Comparison of observed and calculated separation 

factors for diffusion of hydrogen-nitrogen mixture through 

stagnant vapours. ¢g = 0 (stagnant separation agent) ; 

z= 45% He; —-—-— calculated for cyclohexane ; 
- — calculated for methanol. 


nitrogen (z = 0-45). Measured and calculated vapour 
compositions are compared in Fig. 5. 
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Fic. 5. Comparison of observed and calculated vapour 

compositions for diffusion of hydrogen-nitrogen mixture 

through stagnant vapours. q = 0 (stagnant separation 
agent); z = 45% Hg. 


2. Countercurrent flow of separating agent (Table 1, Case 
Ia, II and IV) 


If flow of separating agent is maintained opposite to 
the direction of diffusion of the feed mixture, the heavier 
component is selectively retarded, and the separation 


factor is increased above the value for stagnant separating 
agent (Case I). For this case, a fraction of the feed is 
discharged as raffinate stream, the remainder diffusing 
through the tube ; see Fig. 3. The separation factor is : 


_ vp (l — wn) 
y, (1 — yp) 
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The data for 45% hydrogen, 55% nitrogen as feed mixture 
are summarized in Figs. 6 to 8. Theoretical points obtained 
from the parent relationship given in Table 1 are also 
plotted in Fig. 7. Runs were continued for a time (about 
16 hours) sufficient to establish composition equilibrium 
in product and raffinate streams. Material balances closed 
t. within 5-10 %. 
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SUMMARY OF THEORY 

Solution of the Maxwe.t-Sreran [10] differential 
equations to obtain the steady-state composition 
variation along the diffuser is outlined in Appendix 
1. Simplification of the general solution, equations 
(d), is effected by assuming negligible concentra- 
tion of the components to be separated at the 
product end of the diffuser, a condition which was 
fulfilled also in the experimental work. Table 1 is 
a summary of solutions for the separation factor 
a including some detailed special cases of interest. 
The influence of the important parameter gq, 
defined as the ratio of mass velocities of separating 
agent and of light component, N./N,, will be 
emphasized in the discussion to follow. 

Case I is applicable to the first phase of the 
experimental program, for which g = 0. 
seen to depend upon the binary feed composition, 
z, the three binary diffusivities, D,, etc., the 
mass velocity of light component, N,, the 
diffusion path length, L, total pressure p, and 
temperature, 7. For the special case J (a), where 
components A and B approach infinite dilution 
in separating agent, C, the expression for « 
simplifies to a ratio of diffusivities, D4¢/Dg¢. 
Fig. 9 illustrates the relationship between inlet 
and discharge concentration for the system 
hydrogen, nitrogen, carbontetrachloride. 

The influence of separating agent flow, g + 0 
is given by Case IJ when components A and B 


« is 
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Feed (or product) composition Z | 








oS 10 15-20 25 30 35 40 45 50 55 
Composition of light component % 
(vapour-free) just inside feed point 

Fic. 9. 
hydrogen-nitrogen 


Calculated product compositions for diffusion of 


mixture through stagnant carbon- 


tetrachloride vapour. 


approach infinite dilution in the separating agent. 
Specifically, as countercurrent flow of the latter 
becomes very great, g > 0, the separation factor 
is seen to increase exponentially with the product 
gN,L = NoL, when rpc > r4c- 

Cases III and I V are discussed under “Analytical 
Treatment for Special Cases,”’ p. 222 


Discussion OF EXPERIMENTAL 
AND COMPARISON WITH 
Phase A: ¢ =N, =0 
It is desired to compare the measured separation 
factors with values obtained from equation (2), 
Table 1. Observe that the tube length, L, and 
molar rate, N,, appear together as a product. 
Converting molar rate to volumetric rate : 


RESULTS 
THEORY 


N,L = 0-232 x 10°* y, FL 


where F is the vapour-free flow rate in em* hr at 
latm. and 20°C, y, is the vapour-free product 
composition. By retaining the flow-length pro- 
duct, extrapolation to other lengths and flow 
rates is facilitated. 

The data for cyclohexane and methanol as 
separating agents together with calculated curves 
equation (2), are plotted in Fig. 4. Note good 
agreement with theory over the range studied. 
Average deviation of the measured from the 
calculated separation factors is 3-2% for cyclo- 
hexane and 1-5°%, for methanol. 
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Table 1. Summary of theoretical relationships for the local separation factor, « 





Special cases of equations (d), ApreNpIx 1, assuming negligible concentration of the components to be separated at 
the product end of the diffuser. 


Note — A, B and C refer to light component, heavy component, and separating agent respectively. 
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Note that implicit solutions for N 4 and Nz» in terms of compositions for the case q = 0 have been obtained by E. R. 


GILLILAND, as reported by Suerwoop [10]. 





Case I (a). Ya, ¥B, > ° (Components A and B approach infinite dilution in separating agent). 





Case Ill. ¢g = No/Ng +; Ya, YB, > 0 (Components 


A and B approach infinite dilution in separating agent). 


- exp (— GrpcN4L) 


= - 4 
oe exp (— gr4c N4L) ‘*) 


Case II (a). x» (Countercurrent flow of separating 


agent). 


aia 


a = exp(—q(tac — tac) Na) (5) 


Case II (b). q¢ — «& (Co-current flow of separating agent). 


e = ] (6) 





Case 111. = Np/N4g <1 
Tac + 4p (Dilute mixture of component B in com- 
ponent A on a separating agent-free basis). 


q=No/Ng<-—1; 2” 


1 — exp (— m, N,L) 
c= —— 
1 — exp(— mg N,L) 





(7) 


™ =Tap + 4'BC 


Me = Tac (q + 1) 


Case IV. q = Nc/Np < l; a’ = N4/Np <i. 
(Dilute mixture of component A in component B on a 
separating agent-free basis). 

1 — exp (— mg Ng L) 





ee 
1 oo 
_ At om = 
¥+q4 E (myNgl)  exp(m, Ng i 


1 — exp (— mg NgL) 


’ 771 (10) 





¥B 





Case III (a). m, N4L and my N4L are large quantities : 


a = exp (Am N, L) (8) 


Am = m, — m,> 0 





Case IV (a). m, Ng L and mg NgL are large quantities : 


a = exp (Am Ng L) (11) 


Am =m, — m,> 0 
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The calculated values referred to are based 
upon three binary diffusion coefficients, D,,, 
Dac, and Dyce. From an average of values 
predicted by the semi-empirical method of 
GILLILAND [4], and that of Hirscure.per et al. 
[6], together with some experimental data, a 
value of D4, = 1-02 cm*/sec at 100°C, 1 atm. 
was chosen for the calculations : 

lap = a 8-0 x 10* (em) (sec)/g mole 

P Dap 
Diffusivities of hydrogen and of nitrogen through 
separating agent, D,- and Dg- respectively, are 
calculated from the GILLILAND equation. 


Calculated diffusivities and diffusional resistances 
of H, and N, in separating agents at 100°C, 1 atm. 
| 


| Dac = 0-392 Dac = 0-100 
cm? /sec em 

rac = 782 x 10* | rpc = 30-8 x 10% 
(cm/sec) 





Cyclohexane | 2 /sec 


“g mole 
0-681 Dac = 0-193 
em?; sec em?, sec 
= 449 x 10% | rp¢ = 15-9 x 10% 


Methanol 


| 
| 
| 





Extrapolating the curves of Fig. 4 to FL = 0, 
the limiting separation factors at infinite dilution 
are : 

Methanol: a = 3-45 
Cyclohexane : a = 3-62 


From equation (3), Table 1, the calculated values 
are : 


Methanol : = = 3-58 


Cyclohexane : a = Dac _ 3.99 
Dac 

A further test of the theory is a comparison of 
the observed vapour compositions at the inlet end 
of the diffuser with corresponding values pre- 
dicted by equation 2(a). Good agreement as 
shown in Fig. 5 is evidence of the effectiveness 
of the glass wool packing in eliminating mass 
convective mixing. 


Phase B : Countercurrent flow of separating 
agent: q< 0. 


For this case a fraction, h, of the feed gas is 
made to diffuse through the tube. In addition, 
the separating agent is made to flow counter to 
the direction of diffusion of the components of 
the feed mixture. It is the effect of this flow of 
separating agent that is to be emphasized in the 
following discussion. 

Let the volumetric flow rates (20°C and 1 atm.) 
of the product, feed, and raffinate streams be 
designated by appropriate subscripts: Fp, Fp, 
F,. For convenience, the flow-length product, 
F pL, is retained. It is desired to investigate 
the effects of FpL, of the separating agent/pro- 
duct ratio of g, and of the per cent recovery 
pp (= 100y, h)/z on a. 

The data for 45°% H,—55% N, feed mixture 
diffusing through methanol and cyclohexane as 
separating agents are plotted in Figs. 6 to 8. From 
Fig. 6, it is evident that large enrichments can be 
attained by increasing the ratio of separating 
agent to product flow, —g. This increase is 
achieved at the expense of reduced throughput, 
however. The values of — g that existed in the 
various experiments are listed beside the plotted 
points. For comparison, the lower, dashed, curves 
for g = 0 are replotted from Fig. 4. 

The effect of countercurrent carrier flow is also 
demonstrated in Fig. 7, in which « is plotted 
against — q directly, and the data compared with 
equation (1). The experimental points scatter, 
with maximum deviation of the measured from 
the calculated « of 25%. The averages are : 





Cyclohexane Methanol 


Average deviation of a: 10% 90% 


Variance of « : 16-0 % 150% 








In view of the uncertainties in measurement, 
these deviations do not seem excessive, and point 
to a general corroboration of the theory. 

At large values of — gq, « is seen to fall off 
slightly. This is a consequence of the increase 
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in flow of the diffusing components which resulted 
from the experimental method of maintaining 
flow. The calculated curve takes this increase in 
rate into account. 
Fig. 8 illustrates the rapid decrease in hydrogen 
recovery, 4, With increasing a, as is to be expected. 
General conclusions : 

(1) Large separation factors are associated with 
low product flow rate, small percentage 
recovery, and high separating agent-to- 
product flow ratio, although the correspond- 
ence is complicated by interdependence of 
the variables. 


Comparison of the data with theoretical 
calculations for countercurrent flow is within 
limits of experimental error. It is felt that 
this agreement, together with agreement for 
stagnant separating agent conditions con- 
stitute verification of the theoretical approach. 


TREATMENT FOR 
CASES 


ANALYTICAL 
SPECIAL 
Substantial simplification of the parent relation- 
ship is possible when one component of the feed 
mixture is present in small proportions. Case III, 
Table 1, is applicable when the heavy component, 
B, is dilute. To give some idea of the validity of 
this simplified relationship, anumberof calculations 
of « was made for a range of values of n, g, and N,, 
and the results compared with the accurate 
equation (1). It was concluded that a lower limit 
of — q of 2 should not be exceeded for the error 
to be less than 30%, with n as great as 0-1 (10% 
of the heavy component). 

Fig. 10, which is a graph of equation (7) 
illustrates the increase in « effected by (1) increas- 
ing the ratio, — g, of vapour to light component, 
at constant FL or (2) increasing FL at constant 
— 4. 

When the light component, A, is dilute, Case 
IV, Table 1, is applicable. Use is made of this 
relationship in Appendix 2 to derive an inte- 
grated expression for the performance of a large 
collection of parallel diffusion “ pores,”’ as would 
be employed, for example, in a mass diffusion 
stage of the type described by Marer [8] and 
Benepict [1). 




















Fie. 10. Calculated separation factors for diffusion against 

methanol or cyclohexane vapour-feed of a dilute mixture 

of nitrogen in hydrogen. —-— methanol ; — cyclo- 
hexane (n < 0-1). 


SUMMARY 


Experimental measurements of diffusion rates 
and composition gradients in a long tube diffuser 
lend support to the three-component diffusion 
equations obtained here from the Maxwe.t- 
STEFAN differential equations. Application to gas 
separation by extraction of the heavy com- 
ponent with counter flowing vapour emphasizes 
the role of the latter in effecting relatively large 
composition changes in a single pass. 


APPENDIX 1 
Theory of steady-state diffusion in a ternary gas mixture 
The generalized MAXWELL-STEFAN equations (10) are 
assumed to apply : 


_ % . 


de : (a)* 


Py 
D.. (4; u 


i) 


where p,; and u; are the partial pressure and diffusion 
velocity with respect to a space-fixed co-ordinate system 
of component i, 2 is the position coordinate in the direction 
of the u;, Dy; the ordinary diffusion coefficient of i and j, 
and p the total pressure. It is to be noted that viscous 
effects and thermal diffusion are neglected, and that the 

* It has been pointed out that equation (a) is a special 
case of a general relationship derived by Curtiss and 
HIRSCHFELDER [7] by assuming : 1 No thermal diffusion ; 
2. No pressure diffusion ; 3. No forced diffusion ; 4. Con- 
stant temperature and pressure. 
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total pressure is constant, Furthermore, Dy(- D;,) is 
assumed constant, independent of concentration, and is 
further assumed to be equal to the value for binary 
diffusion, as a first approximation. 


Following the nomenclature of Suxrwoonp [10], let 


(g — mole) Py RT : ; ' 
4 (sec) (cm?) = RT ity= . Dy’ a kind of diffusional 


resistance. Hence 


dpj y” , , 


Here the N; are all constant for the steady state. 

When equations (b) are solved simultaneously, the 
following pair of equations is obtained for the case of 
three components, A, B, and C : 


d? p 4/da* + K, dp4/dx + Kapa = Kay 
d® pp/da* + K, dpp/dx + Kgp4 = Kgp («) 


where the K’s are functions of the ry, and N;’s and are 
independent of z. 

Consider simultaneous diffusion of components A, B, 
and C in a tube of length L. the cross-section and mass 
fluid velocity being taken sufficiently small that com- 
position differences perpendicular to the direction of 
diffusion are negligible. The solution for this case will be 
referred to as the point or local diffusion equations. Sup- 
pose as boundary conditions that known compositions of 
A and B (y4,, ¥p,) be maintained at the discharge end of 
the diffuser.z = L. It is convenient to introduce the flow 
ratios : 

n= Np/Na, 4 = Na/Ne 


Solving equations (c) by the usual method for differential 
equations with constant coefficients, noting that y4 = y4,, 
Yp* YR, atx =Q@:; 

7 
(@ — mg — nd) 
va, (0 — ™)— 98,6 —-— 


wa, ' Am exp (m, Nig L) 





(0 — m, — nd) 
ya, (9 — ™) na ~ 2 —koor 


n+q-7+il1 





Am exp (me N4 L) 


@- 
(8 —m) |v ~™3)—¥n.d — 





¥B, = 4 Am exp (m, N4L) 

6- m,- ng 
n+q+1 
¢@ Am exp (mg NL) 


(9 —mg) Pac -™)—Y¥p,o — 





tater 





7 


Where : 
my, Me =(— K, +4/ K,? — 4Kq)/2N 4; Am 
K,/N4g =(n + I)rap t+ @ 4 


m) _ me 

lr4ac +(" +9) Tac 

Ky/N =(n +49 + V[tantac +" ap'ac +497 actac] 
@=(@+1)t4c + "tapi ¢=TaB— aC 


Equations (d) are the steady-state solution of MaxweE .u’s 
ternary equations corresponding to the boundary conditions 
given, and solved for the upstream compositions. 


Special Cases 
In order to render the solution more workable, y 4, and yg, 


will be set equal to zero, a case frequently encountered 
where two components are removed at a boundary by 
reaction, rapid dilution, etc. For this case, 


ne 1 0 - —=- nd 
YA, = n+q+1 | Amexp(m,N4L) 


0 — mg — nd 
™ Amexp(m, N4L) 


YR, _ n + q +] ¢ Am exp (mag Nu L) 


(9 — mg — nd) (8 — m,) 
- n 
exp (m, NL) + } 


1 | 1 eters 








Since, in most applications involving material transfer, 
it is desired to know the relative concentration change 
between two boundaries, the “ separation factor,” a, is 
introduced : 


¥4/¥B, 1p, 
= ~—_ . 


YA,/UB, nA, 
Hence, from equations (e) 


(@—m, —nO)(9 —mg) at (9 —mg—nd)O —m,) 





exp (m, NL) 


—_—— +ng¢Am 





- exp (m_N 4L) 
no @—m,—ndé O—m_—n¢d. 
exp(mgN4L) exp(m,N4L) ws 





Am 


This is the parent relationship for « as given in Table 1 
of the text, and from which the four detailed special cases 
are derived. 


APPENDIX 2 


Integration of the local diffusion equations for isotopic miz- 
tures 

For large scale separations, particularly separation of 
isotopes, economics dictate the use of multiple stages with 
large cross sectional area and short diffusion path. These 
requirements can be met by utilizing a porous partition to 
separate feed and separating agent streams, with diffusion 
occurring across this partition. 

If the diffusion partition is visualized as consisting of a 
very large number of parallel pores, each of the type con- 
sidered in this work, the local diffusion equations can be 


(f) 





J. J. Keyes, Jr. and R. L. Picrorp 


integrated over the whole area to obtain the overall stage 
enrichment and capacity factors. Benepict and Boas [1] 
have derived this integrated relationship for the case of 
isotopic mixtures where the diffusion coefficients D 4¢ and 
Dgc are nearly the same, with the added assumption of 
(a) balanced diffusion, wherein Ny = (b) 
¢=T4c —Tpc— 9. It is desired to outline a solution 
for isotopic mixtures which does not require either of 
assumptions (a) or (b), but which assumes negligible con- 
centration of feed components on the back or separating 
agent side of the diffusion partition. For further simpli- 
fication, let it be assumed also that the concentration of 
light component in the feed gas is small, as would be the 
case, for example, in the stripping section and adjacent 


Ny. oF 


to the feed point in a cascade used to separate isotopes 
having low natural abundance of the light component. 

It is desired to relate the overall stage enrichment factor, 
«, and the stage screen area, A, to the external flow con- 
ditions, namely the molal feed rate, F’, 
raffinate rate, R. 

The following assumptions will be made : 
(1) 


and the molal 


Mole fraction of light component A on a separating- 
agent-free basis is small in comparison with unity, 
say < 5%. 

Mole fractions of components A and B on the back 
side of the mass diffusion partition are negligible in 
comparison with unity. 

q = Nc/Np « 1, where Ne and Ny are the molal 
mass rates per unit area of separating agent and 
heavy component through the partition 

Dap * Dace Pw 

Feed gas to stage is free of separating agent. 
Countercurrent flow of feed and separating agent. 


Toto! oreo 





P Yep. Yap < * Seg?! 

product separating 
ogent 

Fidg IB, 


ne 


arn’ “er 
Raffinoate 





Fic. 11. Diffusion stage nomenclature. 


Fig. 11 illustrates the main stage variables and nomen- 
elature. The procedure for solution is as follows : 
(1) The integrated ternary diffusion equations for the 
single pore, Table 1, Case 11’, equations (9) and (10) 
are solved for the molal rates N 4 and N » as a function 
of diffusion coeflicients and upstream point composi- 
tions. 
Light 
differential material balances on 


eliminated from 
A and B, and the 
resulting equations integrated over the diffuser area 
to give relationships between flow rates and com- 
positions. 


component composition is 


(3) Results obtained in (1) and (2) are combined with the 
definition of an enrichment factor, and with an overall 
material balance on component B. 

Let the overall stage separation factor, based on product 


and feed compositions be defined as 


Yap / 24 


J 


YBp/ 1 — a 


For an even split of feed flow between product and 
raffinate streams, the usual cascade arrangement, it can 
be shown that In (2 fh €, 
consideration, « < 1. 


since, for the case under 


In Benepict’s notation, with y 


separability), the desired solution for the overall stage 
enrichment becomes : 
€ , (in 2R/ F’) (= F 
=in?2? — -— 
(in R/ F’) In 2 


The expression for the total stage area, A is : 


A (i R/F) 
F L TRC Pa In2 
In — 
In2Rk/F 
where L is the total diffusion path length, including screen 
and boundary layer thicknesses. 


For the diffuses 
boundary into the raflinate stream, there is obtained, 


separating agent which across the 


the dimensoinless stage parameters «/y, 





P-S/F 


Stage enrichment ratio, «/y 
Stage orea ratio 


Separating agent /feed ratio 





+ 


3456768 
Reject/feed ratio 








VO=NURVATS OFS = 


101) 1o1a 14) 
Ri 


Fig. 12. Calculated performance of an extractive diffusion 
stage for isotopes. 
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A/F’ Lrpc, and Py are expressed as functions of the N4, etc. = molar mass velocity of component, A, etc. 
raffinate /feed ratio, R/F’, alone. Fig. 12 illustrates the g mole 
relationship graphically. (sec) (em?) 
Pyin COTresponds to the amount of separating agent N,,/N 
j w , ; > H ai 7’ 4 B Ny 
which actually diffuses into the reject stream. The total ‘= N.IN 
= A RB 


carrier required to reduce the back pressure of A and B 
pressure, atmosphere 


to zero is, of course, infinite. However, it is reasonable = S/F; P’ =So/F 
to expect only a small change in « if the concentration of N IN ——— 
agent in the product, YCp is reduced to 0-8 or 0-9 mole ae hap A 

= Nc/Np 


fraction. Accordingly, lines are drawn representing these = RT/p Dap, ete. 

conditions and can be used as a first approximation to the (cm)(atm.) 
relative amounts, P of separating agent required. The gas constant = 82-] ————____ 
levelling off of «/y at large values of R/F’ contrasted ( masts) ('S) 
with the nearly linear rise of the area function would 
seem to indicate that the optimum R/F’ will not be very _ Dac (Pac — Dap) 

large, possibly 2 << R/F’ < 5. Dgc(Dac — Dap) 

: ) = stage separating agent flow rate, g mole/sec 
Acknowledgement—This research was conducted . te ao oe i agent wttn ditesion 
under an Atomic Energy Commission Pre-Doc- boundary, g mole /sec 
toral Fellowship grant. The authors wish to = absolute temperature, °K 
express appreciation also to the Union Carbide diffusion velocity, em /sec 


Nuclear Company for permission to publish * distance coordinate, cm 
mole-fraction of component A on a separa- 


Appendix II. ting agent-free basis 
- mole fraction of component A, etc. 


NOTATION = feed composition, mole fraction 
A= | stage area, cm? .= 
emer . on _ - ; = local separation factor = vi ( v1) 
D 4p, ete. = mutual diffusion coefficient for pair A-B, yp (1 — yp) 
; em? local separation factor at infinite dilution 
ete., — 
sec 


stage raffinate flow rate, g mole/sec 


: , separability = Dac -1 
volumetric flow rate (1 atm., 20°C) of diffu- 2 BC 


sing components, cm® /hr =T4ac — "Bc 

stage feed flow rate, g mole/sec m, — Mg 

fraction of feed gas diffusing stage enrichment factor 
(mn + 1)r4p+ @ + 1) t4c0+("h +9) TRO, =(¢+1)r4c+"r4ap 
(em) (sec) '4B —TAC 

'AB — "BC 


& 6& Da 


g mole 


(n+ q+ 1) [faptac + " Tac "BC Subscripts 


(em?) (sec)* 
(gz mole)? 
1 


light component 
heavy component 
= separating agent 
feed stream 
product stream 
raffinate stream 
inlet end of diffuser 
discharge end of diffuser 


TAC TBC): 


Kyp ‘AR Ny zcm 
K gnc [zr4c + (1 — 2) rpc) N4/s, em 
L = diffusion path length, em 


My, My = (— K, 4 VK; - tK») 2N Rs 


g mole 
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Abstract 
single parameter D, the * longitudinal dispersion coefficient.’ which is analogous to and has the 


The longitudinal mixing of fluids in flow can sometimes be characterized by a 


same units as the coefficient of molecular diffusion. The results of a study of this model show 
that a dimensionless parameter, the Peclet number can be used as the similarity criterion for 
longitudinal mixing. Also, it is pointed out that the obvious and direct method of calculating 
the mean velocity of flow, by injecting a tracer into the fluid stream at one point and measuring 
its maximum concentration at a given point downstream, may in some cases lead to an appreciable 
error, even in situations where the diffusion-type model is applicable. Methods are shown for 
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evaluating D from experimental measurements, examples are worked out and conditions for 
applicability of the model are discussed, 


Résumé— En écoulement des fluides, le mélange dans le sens longitudinal peut étre caractérisé 
par un paramétre unique D: ™“ coefficient de disy ersion longitudinal ” analogue au coefficient 
de diffusion moléculaire et exprimé avec les mémes unités. Les résultats d'une étude de ce 
type montre que le nombre sans dimension de Perclet peut étre utilisé comme critére de similitude 
pour les mélanges “* longitudinaux.”” De méme les auteurs signalent que la méthode classique 
et directe pour la mesure de la vitesse d’écoulement d'un fluide par introduction d'un traceur 
et mesure de.sa concentration maximum en un point donné, peut en certain cas donner des 
erreurs appréciables, méme dans les cas ob le type de diffusion modéle peut étre appliqué. Ils 
indiquent des méthodes pour calculer D a partir de mesures expérimentales, donnent des exemples 
avec les conditions d'utilisation, 


LONGITUDINAL self-mixing is of interest in the to the problem of fluid mixing in pipes, packed 
study of fluids in flow and it is only in the last 
few years that significant progress has been made _ the 


diffusion. 


beds, etc., an obvicus first approach is to consider 


analogy between mixing and molecular 


in the analysis and understanding of this pheno- However, to do this we must have a 
menon. With respect to an analysis of fluid 
self-mixing the various physical situations may 
be placed in two broad categories depending on 
whether there is or whether there is not macroscopic 
variation in flow conditions at different locations 
along the flow path. We may consider that no 
such variations occur in flow through long pipes 
and long packed beds which may be characterized 


counterpart to the molecular diffusion coefficient 
which we shall call the “ longitudinal dispersion 
coeflicient.’”” The model we propose to explore 
assumes that the longitudinal dispersion 
coefficient, D, as defined in equation (1), is 
independent of position and provides a measure 
of the degree of mixing occurring during flow. 


We use the word “ longitudinal ” above because 


as infinite in length. On the other hand significant 
variations may be expected to occur in flow 
through tanks, short fluidized _ beds, 
filter cakes as well as through and around obstruc- 
tions. This paper is primarily concerned with the 
case characterized by the infinite pipe. 

In an attempt to provide an analytical solution 


pipes, 


22 


we wish to distinguish mixing in this direction 
from mixing in the lateral direction which may be 
quite different in magnitude. For example, in 
streamline flow of fluids through 
longitudinal mixing is due mainly to the fluid 
velocity gradients while lateral mixing is a result 
diffusion. | Numerous authors 


pipes the 


of molecular 


‘ 
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[2, 3, 6, 7, 8, 9, 18, 14, 15, 16] have pointed out 
that longitudinal mixing can be treated like 
diffusion and much has been done to develop 
the consequences of such a treatment. 


THEORY 


Consider a fluid flowing v ft*/sec through a 
pipe which stretches to infinity in either direction. 
Let the fluid enter an experimental section of this 
pipe at 2 = 0 and leave at a = L, the volume 
of this section being V ft’. 

Now at about time @ = 0 rapidly inject into 
the entering fluid a quantity of tracer which at 
unit concentration would have a volume Q ft*. 
If the flow rate and degree of mixing are indepen- 
dent of radial position, the slug of tracer will 
move downstream at an average velocity 
u = vL/V ft/sec and its location at time @ will 
be z =u@. With these conditions the mixing 
of tracer with the surrounding fluid is given by 


(1) 


where D, the longitudinal dispersion coefficient, 
macroscopically determines the mixing process. 
If the tracer is injected essentially instantan- 
eously into fluid at rest at z = 0, the solution to 
equation (1) which gives the tracer distribution 
as a function of @ and z is, from Carsiaw [1], 
— -, * 
2V4/7D0 
a family of normal or Gaussian error curves 
symmetric about the origin, with @ as parameter. 
However, in our case the fluid is not at rest but 
is moving down the pipe; thus equation (2) is 
modified to yield 


ee 
~ 2V4/eD0* 


Cc (2) 


(2—uo)* 
bs ede (3) 
At position z = L the concentration of tracer as 
a function of time is given by replacing x by L 
in equation (3), which on rearranging becomes 
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Equation (4), when plotted on CV/Q vs. v0/V 
co-ordinates, results in a family of curves with 
D/uL as parameter. This is the C-curve [2] 
and it represents the concentration of tracer 
as it passes the downstream recording point. 








1. C-curves for small values of the longitudinal 
dispersion number D/ul (equation (5)). 


Figs. 1 and 2 show the C-curves for various values 
of the parameter D/uL, which is the reciprocal of 
the Peclet number. Note that the skewness 
increases with D/ul ; also for very small values 
of D/ul equation (4) approaches the normal 


Fic. 2. C-curves for large values of the longitudinal 


D 
(equations (4) and (8)). 
aL 


number 


error curve. This can be shown as follows : 
when D/uL is very small the values of v@/V 
must be taken close to 1 (otherwise the values of 
CV/Q are essentially zero); and for v@/V close 
to 1 equation (4) is approximately 
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(Ean 


(9 


1 e (7) : 
2, | D 
ul 
which is normal. 

Physically this means that when D/uL is very 
small the concentration curve, which is normal 
in shape at any instant, does not change appreci- 
ably in the time interval required for the fluid 
to pass by the recording point. However, when 
D/uL is of the order of 0-01 or greater the concen- 
tration curve flattens out rapidly and its shape 
changes significantly during the time required 
for the tracer to pass the recording point ; this 
is the cause of the skewness. 

It is often sufficiently accurate to assume 
(e.g. Danckwerts [2] and TayLor [13, 14]) that 
the maximum value of VC/Q occurs for v@/V = 1. 
This is, in fact, not quite the case, as will now be 
shown. If equation (4) is differentiated with 
respect to v@/V and the derivative equated to 
zero the solution, which represents the time at 
which the maximum concentration occurs, is 


vd D\? D 
; lh ae ; 3 
V (a) . (a (9) 


which can also be written as 


v 1 . 
V iT. so 
Ma) +3 +a 


It is evident from this latter form that the 
maximum concentration always occurs for values 
of v8, V between 0 and 1 and that for very small 
values of D/ul the maximum value of CV/Q 
is achieved when v@/ V is very close to, but less 
than 1. However, if D/ul is of appreciable size 


= 


the value of v@/V given by equation (7) is 


approximately and thus may depart 


1 
2(D/uL) 
measurably from 1. 

When the value of v@/V from equation (6) 
is substituted in equation (4) the result is 


The dotted line in Fig. 2 shows the locus of 
") 
Q / max 


Fig. 3 gives the plot of (S) 
versus D/ul and also | . 


max 
versus D/uL, 

vé 

which from equation (4) is 


Cc r 1 /(D\-3 7 
Q ve , z7z\az (9) 
vy" 
Examination of Fig. 3 shows that the two CV/Q 


values are practically identical for all but large 
values of D/ul. 














Output 


, — , D 
as functions of the longitudinal dispersion number I 
ul 
(equations (8) and (9)). 


Finally, for use in the next section, we develop 
the variance o* of the function in equation (4). 
This is found, using the formula 


b b 2 
= | sede —| [asceyae (10) 
a a 
[17], where x corresponds to v@/V, f(a) is the 
function governing its variation — equation (4) 
in the present case-and the interval from 
a to b is the domain of the variable z, namely, 
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0 to @ for v@/V. Thus, the equation for the 
variance becomes 


= IOP) . (11) 


0 


This expression can be evaluated by first letting 


00 _2+Ve+4, 
V 2 
then later substituting 
2 


Re 
(a) 


This enables one to convert the integrals in 
equation (11) to gamma functions whence it is 
a simple matter to find 


(12) 


=0(2) 
o =S\az) * 


3 
u 


or, solving for D/uL, 


i ee 
= == (V/8o#+ 1 — 1). 


aL (18) 


If D/uL is small enough so that (D/uL)* can be 
neglected, the solution of equation (12) is 


(14) 


EXPERIMENTAL DETERMINATION OF THE 
DISPERSION COEFFICIENT 

In this section methods will be suggested for 
determining the value of D from experimental 
results, following which examples will be worked 
out as illustrations of the methods. These may be 
broadly classified as methods requiring only one 
tracer concentration reading and methods involv- 
ing a series of tracer concentration readings. 

If only one reading is to be taken it may be 
either the maximum concentration or that at 


v@/V =1, whichever is more convenient to 
obtain. It uses the fact that either of these 
values is a function of D/uL; thus, knowing 
Q, u, L and the concentration value, it is possible 
to find D. More data are required if a series of 
readings is to be taken; however, it probably is 
worthwhile, for the results will be more informative 
(e.g., as a check on the suitability of the model) 
and the estimated D will be more reliable. Also, 
one may use either the C- or F-curves (2) directly, 
the F-curve being the output for a step function 
input. If the C-curve is used, the method exploits 
the relation between its variance, o*, and D/uL ; 
again, knowing u, L, and the computed variance 
but this time not requiring Q it is possible to solve 
for D. If the F-curve is used one may either 
reduce it by differentiation to the corresponding 
C-curve, or else if D/uL is small one may plot 
the F-curve on probability graph paper from 
which the variance is found directly. 
Example 1: use of one concentration reading of 
the C-curve 

Consider water flowing at 1-17 ft/sec through a 
1-124” diameter pyrex tube. A volume of 1%, 
KMnO, solution, which would fill one inch of 
tube, is rapidly injected into the flowing water 
and 9 feet downstrean from this point of injection 
the concentration of KMnO, in the water is 
measured by an emission type photoelectric cell 
hooked up to a milliammeter. A camera photo- 
graphs the milliammeter dial 7-7 seconds after 
introduction of KMnO,. This corresponds to the 
time for one test section — volume of fluid to pass 
through the test section. The photograph shows 
that a 0-00555°, KMnO, solution passes the 
recording point. 

From the above data : 


L =9 ft. 

u = 1-17 ft/sec. 

v = 117A ft" /sec. 
V =9A ft’. 

0 = 7-7 sec. 

C = 0-00555%. 


Q = < ft® of 1%, KMnQ, solution, considered 


of unit concentration ; 
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therefore v0/V =1 

and CV/Q = 0-6 

From Fig. 3 D/uL = 0-217 
hence D = 2-28 ft? ‘sec. 


Note that from figure 4 CV/Qinax = 0°64, a value 
which is only slightly higher than that found at 
v@/V = 1. However, from Fig. 2 or equation 
(6) we find that C,,, occurs at 00/V = 0-8, so 
if this method is used to find the flow rate of 
fluids and it is assumed that one volume of 
tracer has passed through the test section when 
Cirax Occurs a 20% error will result. The flow 
condition and D value for this example approx- 
imate Run 194 of Fow.er and Brown [4). 


Example 2: use of a series of readings of the C-curve 


Consider the experimental set-up of example 
1 but with one modification, readings of the millia- 
mmeter dial are taken at 2 second intervals. 
The results for the flow conditions of example 
1 are shown in Table 1. 


Table 1. 





Time KMnO, Concentration 
(sec) (arbitrary units) 


—~—— — — © 8 Me OF 








To estimate the variance of the time-concentra- 
tion curve formula (10) is modified, essentially 
by replacing the integrals by finite sums and by 
replacing the theoretical function of equation 
(4) by the observed concentration readings. 
Thus, the modified formula is 


s 227 (ae a 

% = Zf =f 
where x corresponds to @, f corresponds to or is 
proportional to CV/Q and the summation is 
taken over all the uniformly spaced readings. 
Evaluating the individual terms of equation (15), 


Zf = 386 
2af = 4252 sec 
Z2* f = 65392 sec*. 


(15) 


Hence, from equation (15) 
o} = 48 sec’, 
or in terms of reduced time units v@/V, 
o” = 0-81 


From equation (13) 


(2) oan 


u 
hence 
D = 2.28 ft? /sec. 


Note that neither the quantity of injected tracer 
nor the actual outlet concentration need be known. 
Also, it may be seen from the data in Table 1 
that the time-concentration curve has a long 
‘tail’. This is in accord with the theory which 
indicates that such a situation should occur for 
the D/uL value of this example. Skewness of 
this sort has been found experimentally [13, 14] 
and various explanations have been proposed 
to account for it. 


Example 3: use of a series of readings of the 
F-curve 

To illustrate the use of the F-curve in finding 
the dispersion coefficient consider the data of 
von RosEnBERG [11] for the displacement of 
benzene by n-butyrate in a 14” diameter packed 
column. Take the run at the lowest flow rate, 
2-19 x 10° ft/sec, which is about 2 feet/day. 
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Concentrations were determined by refractive 
index measurements and when plotted, the time- 
concentration curve was found to be S-shaped. 
However, when plotted on probability paper the 
raw data [12] for this run fall approximately on 
a straight line as shown in Fig. 4. 











10 20 S080 BOOS Oe 
n-Butyrate in outlet ——— 
Fic. 4. Data of von Rosenperc plotted on probability 


(or normal curve) paper (example 3). 


The standard deviation, o, of a norma! distribu- 
tion which plots as a straight line on probability 
paper is found by using the following property 
of a normal distribution : namely, that the interval 
between the 16th and 84th percentile points 
is two standard deviations [17]. The 84th per- 
centile line intersects the straight line at 187750 
secs and the 16th percentile line intersects it 
at 178550 secs, so the difference, 9200 seconds, 
is taken as the value of 2c. Therefore 


a, = 4600 sec 


But the time necessary for one pore-volume of 
fluid to pass through the packed bed is 183,000 
secs. Therefore, in terms of v@/V 


o = 0.0252. 
From equation (14) 
D/uL = 0-00082 


But u = 2-19 x 10° ft/sec 
and L = 4ft, 
so D=28 x 10° ft?/sec 


SMITH 


which is in the order of magnitude of the mole- 
cular diffusion coefficient. 

This method should only be used when the 
F-curve plots as a straight line on probability 
paper, or equivalently, when D/ul is small so 
that the C-curve is approximately normal. 
As D/uL < 0-01 in this example, we are justified 
in using this method. 


DISCUSSION 


In this treatment, D has been considered to 
be independent of position. This holds for either 
of two conditions, the first being that in which 
fluid velocity and mixing are uniform. In flow 
through a bare pipe this condition certainly does 
not hold at low Reynolds numbers but is more 
closely approached at high Reynolds numbers. 
In packed beds of solids it may be expected that 
these conditions hold fairly closely, although 
recent experiments [5, 10] show that there is an 
appreciable variation in flow velocity with 
lateral position. 

There is another condition for which we may 
apply this model even if the velocity is far from 
uniform, namely, when the lateral dispersion of 
material (for instance, by molecular diffusion or 
turbulence) is great enough to ensure a uniform 
tracer concentration at any given cross section. 
In determining when this holds, 
the diameter and length of the experimental 
section must be considered as well as the mixing 
of material. For example, for given flow condi- 
tions the longer the experimental section the 
smaller the longitudinal mixing effect will be 
relative to the total amount of fluid. However, 
the lateral mixing is unaffected by such changes 
in length of experimental section. As a result, 
the longer the tube the greater is the lateral 
mixing with respect to the longitudinal mixing 
and the more reliable is the model [13, 14]. 

As D is independent of position when either 
of the above two conditions is approximated 
we may expect this model to be fruitful in study- 
ing conversions in homogeneous flow reactors 
(high flow rate, hence condition 1 holds), in caleu- 
lating the degree of contamination in the succes- 
sive flow of fluids in pipe lines (large length to 


condition 
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diameter ratio, hence condition 2 holes) and 
finally in estimating the depth of the frontal 
mixing zone between fluids when an oil reservoir 
is being invaded (turbulences restricted to the 
size of the pores, large length to pore diameter 
ratio, hence both conditions hold). 

If varying flow conditions exist in different 
parts of the experimental section, an average 
D value will result. For example, consider flow 
through a reaction vessel in which the tracer is 
introduced at the inlet and the C-curve is obtained 
from the outlet stream. In this situation, end 
effects may be appreciable since flow conditions 
close to the injection and recording points may 
differ widely from those found in the rest of the 
equipment. The magnitude of these effects would 
depend on both the actual hydrodynamic situa- 
tion and D/uL for the vessel. For small D/ul 
end effects can be ignored but for large D/uL 
these effects should be considered, due to their 
distortion of the expected C- and F-curve. 
As yet no treatment of end effects is available 
and, until such time that their magnitude is 
known, they will have to be ignored. 

It should be noted that the longitudinal and 
lateral dispersion of material are intimately 
related and for flow through a tube this relation- 
ship forms the basis for a theoretical prediction 
by Tay.or [13, 14] of the longitudinal dispersion 
coefficient. Considerable data appear in the 
literature which if reanalysed should yield useful 


correlations relating the longitudinal dispersion 
coefficient to the properties of the system. Such 
experimental results should also be compared 
with Tay.or’s theoretical predictions. A study 
of this kind is being made and will soon be 
reported. 


NOTATION 


cross sccti»nal area of tube (L?) 

= concentration of tracer in the fluid (dimen- 
sionless]} 
initial concentration of tracer in fluid (dimen- 
sionless) 
the longitudinal dispersion coefficient, defined 
in equation (1) (L?/T) 


the reciprocal of the Peclet number (dimen- 
sionless) 
length of experimental section of tube (L) 
volume of tracer of unit concentration which 
would correspond to the actual amount of 
tracer introduced into the fluid (L%) 
vL/V, the average flow velocity (L/T) 
volumetric fluid flow rate (L°/T) 

= volume of tube (L*) 


reduced time (dimensionless) 


distance from entrance of the experimental 
section of tube (L) 

time, measured from time of injection of tracer 
into the flowing fluid (T) 

the variance of the C-curve (dimensionless) 


2 
) o”, the variance of the time-concentration 
v 


curve (T?). 
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Book 


©. Kupascnewski and E. Li. Evans: Metallurgical 
Thermochemistry. (Second Edition), Pergamon Press, 
1955, xiv 


. 
410 pp. 55s. 


THe appearance of a second edition of “ Metallurgical 
Thermochemistry ~ indication of both the well- 
deserved success of its predecessor and the continued 
interest in the practical applications of thermodynamics, 
especially in the field of process metallurgy. The basic 
arrangement of the book has not been altered, but it has 
been significantly revised and extended. The theoretical 
treatment in Chapter 1 (which still assumes a familiarity 
with the elements of the subject) now includes a con- 
sideration of recent work on non-regular solutions and a 
section on slag-metal equilibria. The following 
chapters, on experimental methods (with emphasis on 
calorimetry) and the estimation of thermodynamic data, 
have also been slightly enlarged. The main contribu- 
tion to the new edition’s extra 42 pages however, is the 
revision of the tabulated data which form Chapter 4. 
Together with the earlier one on experimental methods, 
this chapter is a forceful reminder that although the 
science of thermodynamics is exact, a precise measurement 
of the thermodynamic properties of a system is often a 
matter of considerable difficulty. Thus, as the authors 
observe, a revised value does not always fall within the 
limits of accuracy assigned to the original one ; continued 
caution in the setting of accuracy limits is therefore wise, 
especially since a value can often be useful despite a fairly 


is an 


two 


large probable error. 

This compilation of data is compact and well-presented, 
and its revision will serve to retain for the volume a leading 
position among works of reference in this field. It is not, 
however, the book's only asset, as those actually engaged 
in the improvement and extension of the data of chemical 
thermodynamics will discover. The chapter on experi- 
mental techniques for example, should do much to eli- 
minate inconsistencies of the type that the authors have 
had to face in preparing Chapter 4. 

The printing and binding are adequate, and the price 
is modest. 


Hi. C. Cowen 


Thermo- 
Press, 1956. 


CATTERALL : 
Pergamon 


O. Kupascnewski and J. A. 
chemical Data of Alloys. 
vi + 198 pp. 45s. 


Tuts book is a critical compilcation of data for a total of 
about 230 binary and ternary systems. It includes in 
addition to purely metallic systems the binaries formed 
by the transition metals with carbon, phosphorus, 
nitrogen, sulphur and oxygen. Some of the data are 


reviews 


included in “ Metallurgical Thermochemistry,” but in the 
present volume cach system is dealt with in a separate 
section which contains, as far as possible, the partial and 
integral values of the enthalpy, entropy and free energy 
changes occurring throughout the system together with 
volume changes and latent heats of fusion and transfor- 
mation. The data have been gathered from over 500 
sources, and include values reported up to 1955. There 
is undoubtedly a real need for such a volume, especially 
in view of the rate at which information in this field is now 
accumulating. Despite this activity however, the authors 
reveal many gaps and uncertainties in our present know- 
ledge. The book will be of value mainly to research workers, 
who will be grateful that so much experience in the handl- 
ing and assessment of such data lies behind this compila- 
tion. 
The book is satisfactorily printed and bound, and 
reasonably priced. 
z&. ¢ 


COWEN 


Combustion Processes. Edited by B. Lewis, R. N. 
Pease and H. S. Tayior. Volume 2 of the series of 12 
volumes on High Speed Aerodynamics and Jet Pro- 
pulsion, Princeton University Press, London: Cum- 
berlege. 1956. pp. xv 662, with 16 plates. 84s. 


Since the end of the war a vast amount of work, experi- 
mental and theoretical, has been done in the field of com- 
Outside the universities the chief stimulus has 
developments of new 


bustion. 
been the tremendous technical 
engines for aircraft and rockets. The series which includes 
this volume has been prepared in order to * bring together 
in convenient, accurate and organized form the basic 
information that has evolved in the last decade in the 
related fields of gas dynamics, aerodynamics, combustion 
and jet propulsion.” 

This volume sets out to treat aspects of combustion 
processes other than physical measurements in com- 
bustion (to which volume 9 of the series is devoted). There 
are 15 chapters grouped into six principal sections : 
1. Thermodynamics of Combustion, 2. Chemical Kinetics 
of Combustion, 3. Flame Propagation in Gases, 4. Com- 
bustion of Liquids and Solids, 5. Detonation Processes in 
Gases, Liquids and Solids, 6. Energy Production by 
Nuclear Reactions. Although the standards of the various 
sections differ, the wide scope attempted gives this book 
a superficial attraction and libraries already embarked on 
purchasing this series may not wish to omit one volume. 
But it is a book with serious weaknesses. 


Firstly, much of the work described in it is already 
out-of-date. The preface states that the first draft was 
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completed in 1951 and this impression is confirmed by 
the duplication of material here in a text written inde- 
pendently by two of the authors and published in 1952. 


Secondly, the value and usefulness of the various 
sections is markedly uneven. Who in 1957 would consult 
this work for information on energy production by nuclear 
reactions’ Who would use in combustion the tensor 
formulation of the so-called “thermodynamics” of 
irreversible processes? Certainly the other authors of 
the volume themselves do not. 


Thirdly, and most seriously, there are too many authors. 
There are eighteen : their separate contributions have not 
been successfully unified. Editors and authors, of whom 
we read in the general preface as “ sacrificing their spare 
time " to write the book, must share the blame. There is 
still need for a book on combustion whether by one man 
or a few which is conceived as «a whole and written as a 
whole-time task. 

Perer Gray 


Waxrer L. Bapoer and Junius T. Bancnero: Intro- 
duction To Chemical Engineering. McGraw-Hill 
Book Co. Ine., New York, 1955. ix + 753 pp. £3 11s. 6d. 


Over a quarter of a century ago Professor W. L. BapGEer 
produced a book entitled Elements of Chemical Engineering 
in association with Dr. W. L. McCase. It had an imme- 
diate success not only because of the shortage of books 
on this new subject but also because infinite pains had 


been taken with the drawings and presentation of facts 
to make the subject attractive to the young engineer. Since 
then Professor BADGER has become a consulting engineer 
and has brought the practical assessment to the academic 
approach. In this new book with Professor BancnEro of 
the University of Michigan he maintains this outlook. As 
the book jacket notes * there is much we do not yet know 
about all the unit operations and there is a difference 
between what is postulated in theory and what may 
actually be encountered in practice.” Items not covered 
before are “ dimensional analysis, orifices, rotameters, 
modern evaporators and Ponchon diagrams.” 


In dimensional analysis the treatment consists of 
selecting a number of properties and arranging these in 
certain dimensionless groups. This gives the beginner the 
wrong idea. Most of the groups used today have been 
derived from the standard mathematical equations for 
hydrodynamics, heat and mass transfer, etc. This line 
of thought leads for example to persistence in using the 
Dittus and Boelter equation (which is only a slight modi- 
fication of the original Nusselt equation of 1910), instead 
of realizing that, with a slight transposition, an equation 
can be produced which is exactly analogous with the 
empirical equation relating the friction factor with the 
Reynolds number. The difficulty of revising a book on a 
quickly expanding subject is that, sooner or later, the 
time comes when the whole concept has to be altered 
and perhaps, when this time does arrive, it is better to 
start from scratch and write a new book. 

M. B. DonaLp 





